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A 625 gallon powder-insulated Dewar has been built and appropriately instrumented 
to provide information concerning pressurization gas consumption, ullage pressure, liquid 
level (ullage volume, indirectly), horizontal and vertical temperature surveys throughout 
the test fluid, and lapsed time. Pressurization levels of 25, 75, 125, and 180 psig have been 
examined individually. Two techniques—one based on the case of one-dimensional heat 
flow in a semi-infinite solid and the other founded on an electrical network analog of the 
thermal system—are given for making reasonable predictions concerning liquid-phase hydro- 
gen temperature distribution in well-insulated, rapidly-pressurized storage vessels. A ther- 
mal analysis is presented which indicates the mechanics of heat transmission and adsorption 
in the fluid. 


1. Introduction top dished heads as shown in the drawing. The 


diameter of this cylinder was chosen to provide a 


The so-called temperature stratification phenom- 
enon, which has been observed in containers of cryo- 
genic fluids under certain conditions, is of sufficient 
interest in the field of cryogenic engineering to war- 
rant special investigation. Areas in which the 
results of such research may be applied profitably 


manway into the test container for relatively easy 
access to the internal instrumentation. Evacuated 
perlite is the insulating medium, and filtered vacuum 
pump connections are located on both the top and 
the bottom dished heads; flanges attached directly 
to the base of each evacuation port permit filling and 


include those concerning the transfer of liquefied 
gases by pressurization, the transfer of liquefied gases 
by pumping, and the sealed storage of liquefied gases. 
In the first case, the stratification of contents in a 
vessel used to supply transfer fluid permits the utili- 
zation of subcooled liquid, the properties of which are 
desirable from the aspect of single-phase flow charac- 
teristics. In the second case, stratification provides 
a means to attain, and sustain, Net Positive Suction 
Head for transfer by pumping, the absolute magni- 
tude of NPSH required here depending upon both 
fluid and pump characteristics and the geometry of 
the pumping situation. For the final case, the effect 
of stratification is undesirable, generally, due to the 
accompanying promotion of an excessive rate of 
pressure rise in the storage vessel ullage space; the 
undesirability is pronounced for sealed ‘no-loss”’ 
transport containers in which the maximum time 
permissible in transit is a direct function of the con- 
tamer design pressure. 


draining of the insulation space. A static Dewar 
vacuum of 10 to 15 u is achieved in the cold, hydro- 
gen-filled Dewar. 
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2. Experimental Equipment 





An elevation assembly drawing of the 625 gallon 
Dewar used in this study is presented in figure 1. 
The inner vessel is supported within the outer shell 
by means of a 15 in. diam stainless steel cylinder 
which is fastened between central portions of the two 
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FiaurE 1. 625 gallon test Dewar. 
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Essentially all critical test components are 
attached to the manway cover, thereby providing an 
orderly means for placing, removing, and examining 
internal equipment, when necessary, with a minimum 
of difficulty (fig. 2). Horizontal, bright stainless 
steel baffles are soldered to the fill, vent, and instru- 
ment lead tubes at regular intervals within the access 
cylinder to minimize hydrogen gas circulation and 
thermal radiation in this region. One common 
method of insulating a zone such as this, which may 
be exposed to temperatures approximating 300 °K 
on one side and something less than 100 °K on the 
other, is to build an evacuated-powder or high- 
vacuum insulated plug with a cylindrical stainless 
steel wall. However, in the case of a moderate or 
high-pressure rated container, the plug must be built 
to withstand vessel design pressure—implying rela- 
tively heavy plug walls. As a consequence of the 
thermal conduction characteristics of various ma- 
terials involved, the hydrogen gas-baffle arrange- 
ment provides a better heat barrier across the man- 
way than would a plug of adequate pressure rating. 
For this reason, and for the sake of simplicity and 
economy, the baffle construction was used. It 
may be of interest to note that the thermal conduc- 
tivity of hydrogen gas is approximately one hundred 
times that of evacuated powder and one-hundredth 
that of stainless steel. 

Located immediately beneath the bottom opening 
of the access cylinder is a gas diffuser which diverts 
the incoming pressurization gas 90° through curled 
brass shavings and introduces the stream tangentially 
to the underside of the inner shell top dished head. 
The diffuser is fed by a 114 in. tube which also serves 
as a gas vent when the test vessel is being filled or 
depressurized., 

A capacitance probe, hanging from the bottom of 
the diffuser, is used to provide the liquid contents 
(level) determination—the signal of the probe being 
furnished to a direct-reading recorder. Auxiliary 
level information may be obtained by observing the 
Dewar ullage pressure, by then noting the saturation 
temperature corresponding to this pressure on a 
hydrogen vapor pressure-temperature diagram, and 
finally by locating this temperature in data provided 
by the vertical thermocouple rake. Knowledge of 
the liquid surface location derived by the latter 
method has proved to be quite precise in the zone 
of closely-grouped thermocouples; in other regions, 
the precision is understandably lacking. 

Further examination of figure 1 reveals a Venturi 
tube which is attached to the inlet of the fluid dis- 
charge tube to permit evaluation of such a device in 
cryogenic media. Results of this work are presented 
in another paper.‘ 

In the experimental program, two thermocouple 
rakes are utilized—one installed vertically in the 
center of the vessel, and the other fastened to the 
side wall, capable of sensing horizontal temperature 
profiles from zero to 11 in. from the wall at three 
distinct levels (approximately 17, 49, and 56 in. from 

4J. R. Purcell, A. F. Schmidt, R. B. Jacobs, The Venturi tube as a liquefied 


gas flow measuring device, Advances in Cryogenic Engineering, 5, pp. 282 to 288 
(Plenum Press, New York, N.Y., 1960 


the tank bottom, respectively ). Figure 3 illustrates 
the placement of each thermocouple in the vertical 
direction. 

The system of differential thermocouples—59 in 
all—is composed of one 36 gage gold-cobalt lead and 
an appropriate number of 36 gage copper leads, each 
of which passes through a rubber compression-type 
seal at the top of the tank. <A reference junction 
for the entire series of ’couples is attached to a small 
parahydrogen-filled vapor-pressure bulb located be- 
low the inlet of the Venturi; reference temperature 
information, therefore, is derived from bulb pressure 
data and a parahydrogen vapor pressure-temperature 
curve. ‘Temperatures, corresponding to the vapor 
bulb pressure readings, are accurate to £0.05 °K. 

Hydrogen pressurization gas is stored at high pres- 
sure in a standard tube trailer of 260 ft® free volume, 
from which it passes through a preset 
recorder-actuated solenoid valve and 
Dewar during pressurization. 


pressure- 


on into the 


3. Test Procedure and Results 


In an effort to present the data accumulated from 


each pressurization run most descriptively, tests 
A-2, A-3, A-4, and B-1 have been analyzed graphi- 
cally through temperature-time and temperature- 


height curves (fig. 4 to 9 Test identification sym- 
bols listed above refer to the conditions of table 1. 
Tests A-1 through A-4 were conducted in the 
following manner: The test vessel liquid inlet con- 
nection was sealed off, and the container was main- 
tained at *% psig for a minimum of 16 hr by venting 
into a hydrogen gasholder. At the end of this period 
(actually the start of the test day), during which 
time the fluid contents had reached an equilibrium 
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FicurE 4. Temperature versus height, hydrogen test 
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state, a complete set of thermocouple readings was 
taken and recorded as the initial condition. The 
vent-to-gasholder valve was then closed, a pressure- 
recorder-actuated solenoid valve between the hydro- 
gen tube trailer and the test container 
diffuser was opened and pressurization of the Dewar 
to a preselected level was effected within 30 see. 
Additional small quantities of gas were generally 
required for 15 min after the initial pressurization in 
order to maintain the desired pressure level, after 
which time the liquid evaporation rate became more 
than sufficient to furnish the needed gas. Rather 
than disturb the temperature gradients being 
established in the test space by continuously releasing 


gas gas 
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FicurE 5. Temperature height, 


versus 


125 psig. 


hydr ogen 


A, O+ time 
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FicurE 6. Temperature versus height, f, 


180 psig. 


hydroge n test A- 


A, O+ time 
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the gas overpressure, the vent was kept closed and 
the pressure was permitted to rise naturally through- 
out the test period. Thermocouple and vapor 
pressure data were recorded every ' hr for 8 hr. 
Termination of the test after 8 hr was accomplished 
simply by deactivating the recorder-actuated solenoid 
valve and depressurizing the tank to gasholder 
conditions for overnight standby. Gas consumed in 
ach of the tests listed above was measured by 
opening enough tubes in the gas trailer manifold to 
exceed, by a small margin, a precalculated demand; 
the temperature and pressure of this gas were then 
recorded, pressurization of the test Dewar occurred, 
and the final trailer manifold temperature and 
pressure were noted. Knowledge of the fluid 
capacity of each tube in the trailer, coupled with the 
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t pressure-temperature data recorded above, then 
+ provided information concerning the amount of 

20 = t tpt tt tt +t pressurizing gas required in each test. 
r) v2 ; 11/2 2 22 3 Pressure build-up test B-1 was effected by com- 


easiand pletely sealing off the experimental chamber and 


Ficure 8. Temperature versus time, hydrogen test A-4, 180 | Tecording the necessary data three times a day over 
j uo. sue ; ‘! 4 4, o . . 
psig. an 80 hr test run. 





TABLE l. 


Test | Final | Ullage 


Initial | pressure * pressure > Pressuriza- Liquid volume 
| (zerotime | (zero time tion gas con- volume (including 
(zero time) | plus % min) | plus7% hr) | sumption | |} manway 
| 


} } volume 





- 
Test al pressure 
| 
| 


psig | psig } psig scf cf } cf 
A-l % 25 | 46 | 110 67.2 | 20.3 
A-2 % | 75 } 102 | 310 63.2 | 24.3 
A-3 % | 125 } 160 465 59.7 27.8 
A-4 4 180 | (95)e | 700 57.3 30. 2 
B-1 % % 10 oe 56.3 31.2 





* Test pressure refers to level achieved with tube trailer pressurization gas. 

b Final pressure refers to level achieved at completion of test (pressure resulting from combination 
of tube trailer gas and gas evolved during test period). 

¢ Final pressure after 34% hr. 
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4. Discussion 


Data obtained from the program described here 
were studied and used in several different ways: 

(1) Pressurization gas consumption data were 
recorded and several attempts were made to correlate 
this information; an insufficient number of vapor- 
phase thermocouples plus the additional complica- 
tion of irregular vessel geometry due to the manway 
volume discouraged extensive work in this area; 

2) Considering system geometry and certain 
properties of the test fluid, metal containing wall 
and evacuated-powder insulation, a thermal analysis 
was made to determine the mechanics of heat 
transmission and absorption in the fluid; 

(3) A method useful in predicting the change of 
temperature with time at any given point in a solid 
was applied to the experimental quantity of liquid 
hydrogen—-the theoretical results were compared 
with experimental information ; 

(4) An electrical network analog of the thermal 
system was constructed and a comparison with the 
theoretical method (3) was made. 


4.1. Pressurization Gas Consumption 


Information pertaining to the consumption of 
pressurization gas, pressure level and ullage volume 
in this series of tests is presented in table 1. One in- 
teresting observation can be made from these data. 
For a perfect gas—using Boyle’s Law (pV=p’V’) 
the amount of gas required to elevate the pressure p 
in a known volume V’ to a specified level p’ may be 
calculated (assuming no condensation, evaporation, 
etc.). Results of the four tests reported here indi- 
cate that the pressurization gas consumed in any test 
can be predicted to an accuracy of 10 percent by 
first using the perfect gas relationship to compute a 
required volume V of gas at standard conditions, 
then doubling the answer. To restate, the amount 
of gas needed to pressurize a vessel of this particular 
geometry is approximately twice that computed from 
Boyle’s Law. Gas temperature, pressurization time 
and diffuser configuration should also be considered 
in the final analysis. Although the correlation is 
crude, it may be of some value in estimating the 
requirements of subsequent work of this nature. 


4.2. Thermal Analysis 


Heat transmission into the liquid hydrogen was 
considered to occur in the following manner: 

(1) The uniform (bulk) temperature portion of 
the liquid is warmed by heat transferred from the 
containing surface and thence by convection; 

(2) The stratified liquid layer is warmed by heat 
conducted from vapor above the interface and by 
condensation at the interface; 

(3) Vaporization occurs at the outer periphery of 
the interface where liquid contacts the warmer 
metal. 

Validity of these assumptions may be tested by 
heat’ balances using internal energy and apparent 
heat transfer information. Taking experimental 





data from the last hour and a half of test A-3 (fig. 


5) as an example, calculations indicate a_ total 
internal energy change rate of 53.8 w: 
Internal Energy Change Rate 
In uniform temperature portion of liquid (below 
a 17.5 w 
In stratified portion of liquid (above 38” to 50.3"). 24.0 w 
To raise mass of liquid to be vaporized to satura- 
tion temperature ____- -- E 6.5 w 
To raise mass of liquid to be vaporized from satu- 
rated liquid to final vapor state : 22.3 w 
70. 3 w 
To cool initial mass of vapor to final state —16.5 w 
53. 8 w 


In determining apparent heat transfer through the 
evacuated-powder tank insulation and down the man- 
way, a Value of insulation conductivity (k=40 pw/em 
°K) was chosen such that the heat transmission 
computed from Dewar geometry (below) equalled 
the computed rate of change in internal energy 
(above). The reason for choosing a value of & in the 
manner described relates back to some uncertainty 
existing in the absolute value of this property, even 
with a fair knowledge of material substance, size 
and interstitial pressure. The value computed is 
within the range of values considered reasonable for 
the material. 


Heat Transfer Rate 


Through bottom head of test vessel : 10.7 w 

Through cylindrical wall of uniform temperature 
portion of liquid (below 38’’) : 14.2 w 

Through cylindrical wall of stratified portion of 
liquid (above 38’’ to 50.3’’) : 7.3W 

Through vessel wall surrounding vapor phase (above 
50.3") i 8.1 w 
Through manway (hydrogen gas and metal neck) 13.5 w 
53.8 w 


Checking premise (1), The uniform (bulk) temper- 
ature portion of the liquid is warmed by heat transferred 


from the containing surface, and thence by convection: 


Internal energy change rate ~ Heat transfer rate 


17.5 w=24.9 w 


Then with respect to the stratified layer, an as- 
sumption is made that heat transferred through the 
wall surrounding this portion of fluid provides the 
energy required to raise the temperature of the liquid, 
which is to be vaporized, from the bulk temperature 
to a saturated condition. 
Internal energy change rate ~ Heat transfer rate 
6.5 w=7.3 W 


An apparent surplus of transferred heat (7.4+ 
0.8=8.2 w), from the two zones (or regions) con- 
sidered, may be assumed to contribute to vaporiza- 
tion at the periphery interface. 

From premise (2), The stratified liquid layer is 
warmed by heat conducted from vapor above the interface 
and by condensation at the interface: 








Condensation rate ~ Internal energy change rate— 
(stratified layer) 


(interface) 
~24.0 w 
=18.5 Ww 
where the heat conduction rate is computed from 
the gas phase portion of curves F and G@ (test A-3) 
near the interface. 

To complete the account of heat distribution in 
the Dewar, from interface up, gas phase temperature 
and pressure information over the chosen hour and 
one-half test interval is used to calculate a system 
mass (gas) increase. This Am is then multiplied by 
the latent heat of vaporization of hydrogen, result- 
ing in a net vaporization energy value of 8.5 w. 
Total vaporization rate, therefore, is equal to the 
condensation rate computed previously plus the net 
vaporization rate computed above, or 


Total vaporization rate=18.5+8.5=27.0 w 


This rate, plus the 5.5 w conducted through the 
vapor to the interface, must be supplied from above 
the interface and from the apparent surplus beneath 
the interface. 


4.3. Semi-Infinite Solid 


The experimental volume of liquid hydrogen has 
been treated as a semi-infinite solid with the liquid- 
vapor interface representing a constant temperature 
bounding face. Since equilibrium temperatures are 
observed throughout the test fluid prior to pressuriza- 
tion, and since the plane surface temperature imme- 
diately after pressurization is also known (saturation 
temperature corresponding to interfacial pressure), 
the uniform temperature of the liquid and the sud- 
denly-changed surface temperature are known. With 
this information, and certain physical fluid properties, 





the case of one-dimensional heat flow in a semi- 
infinite solid may be calculated. 

A solution for the temperature history 7(z, @) 
must satisfy the partial-differential equation of 
Fourier, 

oT OT 
—a — 
o0 Ox" 


By imposing appropriate initial and boundary condi- 
tions and proceeding through the necessary mathe- 
matical sequence, the desired relationship is found 
to be 


= 9 "12/08 
Cans _. e-* dX 
i 1 y7 /70 
where 
t=temperature at distance z below liquid-vapor 
interface, 
t,=saturation temperature corresponding to pres- 


sure at interface, 
initial fluid equilibrium temperature, 
distance below interface at temperature f 














| 
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Heat conduction rate 
(vapor to liquid) 


-5.5 Ww 





é=time, 
a=thermal diffusivity, =k/pC, 
k=thermal conductivity, 
p=density, 
C,=specific heat at constant pressure, 
2 (7? Va e-g,y—error function, or probability 
vx Jo integral. 


A primary decision one faces in working out the 
mechanics of this problem concerns a some what arbi- 
trary choice in the value of thermal diffusivity. 
Clearly, it must be evaluated at some temperature 
between ¢; and ¢,. For the case of liquid hydrogen, 
the variation in @ is not considerable—the value at 
20 °K being approximately one-third higher than 
that at 31 °K; further, it is va@ that is used in this 
computation. With this knowledge then, and the 
additional stipulation that C, be used as heat cs upacity 
in the diffusivity determination, one may choose es- 
sentially any value of a between ¢; and ¢,, run through 
the simple mathematics of the problem and arrive at 
an excellent approximation of fluid temperature dis- 
tribution. In most the accuracy associated 
with such a method of choosing a is more than suf- 
ficient. If, however, it is desired to refine the caleu- 
lation, the initial approximation (above) may be 
used to provide basic temperature versus distance 
curves, from which any value, or values, of a can be 
selected for the final computation. 

This method of theoretically simulating the actual 
physical conditions described previously is not with- 
out limitation when applied to liquefied gases of small 
latent heat of vaporization and low temperature. An 
efficient thermal insulation was used to envelop the 
test container in this work, thereby permitting rela- 
tively little heat to flow into the fluid. Should the 
liquid be submitted to a progressively higher heat 
flux, degradation of the established temper rature gra- 
dients would be noted at some increased value of this 
flux due to greater boiling, and hence agitation, of 
the liquid. A poorly insulated volume of liquid, fluid 
containers utilizing high-thermal conductivity ma- 
terials of construction shell walls, liners, baffles, 
etc., and Dewars in transit, may not be adapted to 
this type of analysis. 

As implied in preceding paragraphs, the represen- 
tation of a quantity of liquefied gas by a semi-infinite 
solid interesting and apparently feasible. The 
method has been applied to the reported experimen- 
tal data with excellent results. Test data were not 
completely analogous in that the bulk fluid tempera- 
ture shifted slightly with time, due to the combined 
effects of heat transmitted into the liquid from the 
containing wall and compression of the liquid, and 
also because of the fact that the interface temperature 
shifted somewhat during the test period as a result of 
the pressure increase in the nonvented Dewar. 


cases, 
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However, these deviations had only a secondary 
effect on the problem. A comparison of theoretical 
and sexperimental values for the 125 psig test (A-3) 
reveals the following information (see fig. 10): Using 
an interface temperature of 30.8 °K and a liquid 
bulk temperature of 20.3 °K, a temperature difference 
of 10.5 °K has been established. For distances 
greater than 3 in. below the surface of the liquid, and 
for any time between zero and 8 hr, the maximum 
deviation between theoretical and experimental work 

0.3 °K, or 3 percent of the temperature span, 
while the average deviation is less than 0.1 °K. The 
first 3 in. beneath the surface are affected more by 
an increase in interface (saturation) temperature 
with time and by the slight increase in volume of the 
liquid, so that a direct comparison here is likely to 
show the greatest error. Even so, the maximum 
error is seen to be less than 1.5 °K. 


Is 


1.5 


4.4. Electrical Network Analog 


In order to evaluate more directly the effect of 
bulk and surface temperature shift previously men- 
tioned, some time was devoted to the use of an 
electrical analog as a means of most conveniently 
examining this shift. Also, the ability to physically 
observe experimental phenomena requiring 8 hr 
being simulated electrically in a few seconds is of 
considerable interest. 
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The electrical analog is based on the equation 
MT (compare with the Four quati 
== - compare with the Fourier equation 
00 RC ioe Bs : 
oT oT , , cae 
3074 32 previously mentioned), if distributed 


constants are used in the circuitry. Using lumped 
constants—that is, conventional resistors and capaci- 
tors—the error is slight if reasonable consideration 
is given to circuit construction. It is seen from the 
above equations that voltage is analogous to tempera- 
ture in the thermal circuit and that the electrical 
diffusivity 1/RC corresponds to the thermal diffusivity 
a. In the overall circuit, it is not necessary to define 
resistance and capacitance individually—only the 
product which determines the time constant and 
thus the ratio a/1/RC, or aRC. This ratio is gener- 
ally quite large (approximately 20,000 in the present 
case) due to practical limitations of the circuit. 
However, the advantage of such a circuit is under- 
standable in that the curves are plotted very rapidly. 

Results obtained by this technique were exactly 
comparable with those achieved by the method of 
the semi-infinite solid. Unfortunately, circumstances 
did not permit further refinement of the circuit to 
provide for the gradual shift of bulk and surface 
temperatures with time; it is hoped that this may 
be accomplished in the near future. 


5. Conclusions 


It is evident from the preceding paragraphs that 
reasonable predictions can be made concerning liquid- 
phase hydrogen temperature distribution in well- 
insulated, rapidly -pressurized storage vessels. 

By tre ating the liquid as a semi- -infinite solid, one 
may foretell the temperature history, temperature 
gradient and heating rate of the stratified fluid 
portion beneath the “liquid- vapor interface; where 
numerous succeeding initial and boundary condi- 
tions are to be imposed in the manner of trial and 
error, where these conditions are to vary with time 
for any one given situation, or where a variation of 
thermal properties must be taken into account, the 
electrical analog method may be modified to per- 
mit more rapid solutions. 

Additionally, if heat flux through the containing 
vessel wall is known, bulk liquid temperature rise 
may be predicted in a very rough manner. 

Finally, it will be noted by the curves presented 
here, and by theory, that an increase in ullage pres- 
sure—hence an increase of interfacial temperature— 
results in the transfer of more heat to the stratified 
layer, the quantity of which is calculable using the 
method of the semi-infinite solid. Knowledge of 
this may serve as one criterion in the optimization 
of pressure level within a liquefied gas storage or 
transfer system. 


(Paper 65C2-58) 
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Temperature Dependence of Elastic Constants of 
Some Cermet Specimens’ 


Sam Spinner 


(June 16, 1960) 


The temperature dependence of both Young’s and shear modulus of four types of cermet 


specimens, known as ‘‘nickel-bonded titanium carbide,” 


determined by a dynamic resonance method 


as well as Ni and TiC have been 


The Young’s modulus temperature curves of the cermets are characterized by a linear 


decrease from room temperature until about 700 to 1,000 


‘C. In this upper temperature 


region, an inflection in the relation develops, accompanied by an increase in internal friction. 
Both these effects are attributed to viscous grain boundary slip. 
For shear modulus, only the linear portion of the modulus temperature relation was 


obtained. 


In this linear region, the relative decrease in shear modulus for the cermet speci- 


mens is greater than the relative decrease in Young’s modulus for the same type of specimen. 
Chis means that Poisson’s ratio rises with temperature for the cermets. 


1. Introduction 


Among the more noteworthy properties of cermets 
in general are their high values of elastic moduli. 
Although the literature dealing with the properties 
of cermets devotes considerable attention to their 
elastic constants [1,2,3],? comparatively little atten- 
tion has been paid to the values of these constants at 
elevated temperatures.* 

Inasmuch as these cermets are intended to be 
used as structural members at elevated temperatures, 
and inasmuch as these studies [4, 5] have shown that 
their elastic moduli may decrease by as much as 
15 percent at around 1,000 °C from their room 
temperature values, and may decrease by another 
20 percent in the next 200 °C or so, it would appear 
that further data concerning the variation of these 
elastic constants with temperature would be most 
useful. The main purpose of this paper is to present 
such data. Furthermore, since the four types of 
cermets selected for obtaining these data are com- 
posed basically of the same constituents, varying 
from each other mainly in the relative proportions 
of metal and ceramic, comparison of the results 
among these types presents interesting possibilities 
for analysis. 


2. Experimental Procedures 


2.1. Specimens 
The cermet specimens used in this investigation 
r'} t s| ised th vestigatior 
are\known as ‘‘nickel-bonded titanium carbide’’. 
They were taken from among those used by Lang 


[1] for his room temperature measurements. These 





1 This investigation was carried on under the sponsorship of the Division of 
Research, U.S. Atomic Energy Commission. 

3 Figures in brackets indicate the literature references at the end of this paper. 

* A paper by Watchman and Lam [4] containing one curve of Young’s modulus 
as a function of temperature for one of the types of cermets used in this study, and 
a WADC report [5] showing Young’s modulus temperature curves on several 
cermet materials are the only data along these lines known to the author. 
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specimens were all in the shape of rectangular bars 


about 6 in. by 4 in. by } in. Table 1, which gives 
the number of sets of determinations for both 


Young’s and shear modulus performed on all of the 
specimens, retains Lang’s scheme for identifying the 
cermet specimens well the manufacturer’s 
° 4 rT’ bd : 

designation. The table also includes the nominal 
metal content by weight of each type of cermet; 
the remainder (of each type) consisting of a solid 
solution composed of TiC and 8 percent of specimen 
weight of (TiC+NbC+TaC). A description of 
the composition and microstructure of many types 
of cermets including some used in this investigation 
is contained in a report by Shevlin et al. [6]. The 
a a : . ° ‘ 

TiC and nickel specimens, which are included for 
comparison, are polycrystalline materials of com- 
mercial purity. The TiC specimen was cylindrical 


as as 


in shape, about 6 in. long and 3 in. in diameter. 
The nickel specimen was about 7 in. long and 3 in. 


square in cross section. 











TABLE 1. Designation and number of determinations made 
on all specimens 
ee 

Specimens Number of | Number of 

__| Nominal sets of sets of shear 

| metal | Young’s | modulus 

Lang’s desig-| content modulus | determins- 

M’ frs. designation nation =| | determina- | tions 
| tions 
| | | — 
. = —| oe MOE ena - 
: f | 6 2 
K-150 B } 31-4 | 710% Ni | 3 3 
ae 
K-151 B. = 1 |}20% Ni 1 = 
| } 
mas f 33-1 |). . 1 iemandas 
K-152 B..... ee 33-10 |;30% Ni 9 3 
' 
ar , | (25% :— 3 | 1 
K-162 B_.-. 34-1 15% Mo | eocee tale s E E 
bikie — | =. 


Titanium Carbide 
Nickel 





* The number preceding the dash indicates the type of cermet and corresponds 
to the manufacturer’s designation. The number following the dash indicates 
the particular specimen. 








2.2. Method and Calculations 

The elastic moduli measurements at elevated tem- 
peratures were performed by a dynamic resonance 
method as described by Wachtman and Lam [4]. 
Essentially, the method consists of suspending the 
specimen in a furnace from two fibers, one fiber 
being attached to an electromagnetic driver and the 
other to a crystal pickup. The driver and pickup 
are outside the furnace and connected to appropriate 
electronic equipment for arousing and detecting the 
mechanical resonance frequencies of the specimen. 
The evidence of resonance is indicated by a Lissajou 
pattern on a cathode ray oscilloscope. The most 
satisfactory suspension elements were found to be 
molybdenum wires coated with nickel and chromium 
to prevent oxidation [7]. 

If the suspension wires are tied at opposite sides 
of the specimen, as shown in the sketch (fig. 1) from 
a previous paper [S], then, in addition to the usual 
flatwise flexural resonance vibrations, the edgewise 
flexural resonance vibrations and, even more impor- 
tant, the torsional resonance vibrations of the spec- 
imens can be detected. The response of one of the 
flexural resonance frequencies, usually the edgewise, 
usually died out at a lower temperature than the 
other. Also, because of the geometry of the spec- 
imens (the fairly low ratio of width to depth of the 
cross section) the response of the torsional resonance 
frequency was weaker than the flexural ones and 
died out at the lowest temperature. 

Occasionally, the torsional response would disap- 
pear at one temperature and reappear at a higher 
one before dying out completely. The larger num- 
ber of sets of Young’s modulus-temperature deter- 
minations in table 1, as compared with the number 
of sets of shear modulus-temperature determinations, 
is mainly due to the fact that the experimental phase 
of the investigation of the cermet specimens had been 
about half completed before it was noticed that the 
torsional resonance vibrations could be detected at 
all. Nevertheless, at least one set of torsional res- 
onance vibrations was obtained for each type of 
cermet. The flexural vibrations are related to 
Young’s modulus while the torsional vibrations are 
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related to the shear modulus. Because of the square 
and circular cross section of the Ni and TiC spec- 
imens, a satisfactory torsional mode could not be 
obtained at elevated temperatures, so that for these 
specimens only Young’s modulus as a function of 
temperature could be determined. 

The ratio of a particular elastic modulus at ele- 
vated temperature to the modulus at 


Same rool 


temperature, /,//, was calculated from the well- 
known equation, 

E, (i.¥ _ Al 

(7) /G+7) (1) 


where f,/f, is the ratio of the corresponding resonance 
frequencies, 

The factor, {1 A///)| is the correction for thermal 
expansion of the specimen from room temperature 
to the temperature of measurement. Values of A///, 
the linear thermal expansion over the required tem- 
perature range, were obtained for the cermet spec- 
imens from dilatometer measurements. These data 


are shown in figure 2. Corresponding data for TiC 
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FIGURE 2. 
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Thermal linear expansions for materials used in thi 


nvestigation., 


and nickel were taken from the literature [9, 
They are also included in figure 2 for 
with the cermets. 

Equation (1) makes no restriction as to the type 
of vibration frequencies involved in obtaining /,//) 
{12}. The conditions that must be fulfilled are that 
the specimens do not change in shape on heating 
and that f, and f,, in a partie » ratio, both rep- 
resent the same type of vibration frequency, whether 
torsional, flexural flatwise flexural edgewise. 
Thus, if f,/f) are flatwise or edgewise flexural 
onance frequencies, then /2,//) becomes the ratio of 
Young’s moduli, Y,/)>, and if tidlfo are torsional 
resonance frequencies, then /,//) becomes the ratio 
of 7 shear moduli, @,/G,. Equation (1) also holds 
if f,/fo represent overtones of the fundamental of any 
aed ‘of vibration. However, in this study, the 
fundamental resonance frequencies were used almost 
entirely. 

If more than one resonance frequene: v for caleu- 
lating a particular elastic modulus is used, then the 
agreement in /,//) from these resonance frequencies 
provides an internal check on the consistency of the 
data. The agreement in Y,/¥> from the two reso- 


10, 11). 


comparison 


or 


res- 


nance flexural frequencies at any particular temper- 
ature at which both could be detected was on the 
order of 0.1 percent justifying the use of the average 
of these two values for that temperature. For the 


G’,/G) determinations, usually only one ratio of reso- 
nance frequencies was available for a particular 
point (under some unusually favorable conditions, 


the first overtone of torsion could be obtained and 
similarly averaged with the fumndamental), so that 
the agreement of repeated sets of determinations 
served as a measure of the precision of the data. Re- 
peated sets of determinations were available in all 
but (table 1) for G,/G,. For Y,/Yo, a 


one r 
larger number 


case 
of determinations was available in all 
From these repeated measurements, if 
estimated that the precision over the linear range of 
Y’, is on the order of 0.3 percent and that of G, is 
about 0.4 percent. These repeated measurements, 
and consequently the above figures of precision, 
include the data for both specimens of each type 
(when two specimens were used), and no distinction 
is indicated in the notation in figures 3 and 4 between 
the two specimens of each type. 

If G, as well as Y, are determined, then the speci- 


causes. Is 
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FIGURE 3 Young's modulus—temperature relations for the materials of this investigation. 
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Ficure 4. Shear modulus—temperature relations for the four 


types of cermets used in this study. 


mens are completely specified elastically over the 
temperature range over which these two moduli 
can be measured, since, for isotropic materials * all 
elastic constants are related by well-known equa- 
tions in such a way that if any two are known the 
others can be determined. In this study Poisson’s 
ratio, u, was computed from these two moduli. 

In order to plot Y, and G, directly as a function of 
temperature, a room temperature or base value of 
both Young’s and shear modulus, Y) and G) must be 
obtained. For the cermets, these values (of Y) and 
G)) were obtained by averaging the values obtained 
by Lang[1] for the particular specimens of each type 
used. For TiC, Y, and G, were determined from 
the flexural and torsional resonance frequencies, 
respectively, the mass and the dimensions based on 
Pickett’s relations [13]. The value for nickel was 
taken from the literature [14]. The values so ob- 
tained for all the types of specimens used are given in 
table 2. Since Y, and G, are available for each 
type of specimen, Poisson’s ratio, u, may also be 
calculated from the familiar relation, 


Y 
u= 


2G 


This is given in the final column of the table. 

It should be noted that the choice of the base 
values, Yo and Gp, is not critical. Had different 
values of these moduli been chosen on the basis of 
some other reasonable selection, it would not have 
significantly changed the shape of the temperature 
modulus curves shown in figures 3 and 4, but would 
only have caused a slight upward or downward 
displacement of these curves. 


* The materials are considered to be isotropic because, although the individual 
grains of which they are composed are antisotropic, nevertheless their orienta- 
tion and distribution are random, so that macroscopically this condition of iso- 
tropy is substantially fulfilled. Actually, some of Lang’s{1] results indicated 
that some preferred crystal] orientation or inhomogeneity might exist in the speci- 
mens. However, these departures from isotropy, if they exist at all, are small 
and do not significantly affect the relations based on the assumption of isotropy 
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TABLE Room temperature base values of 
specimens used in this investigation 


Specimen | Young’s mod-} Shear modu- Poisson's 


type ulus ¥i lus Go ratio uw 
kilobars kilobars 
3l 4074 1669 0. 220 
32 4077 1681 213 
33 3791 1550 223 
34 3952 1701 162 
ric 1774 748.9 184 
Ni 2147 


2.3. Oxidation 


The specimens oxidized slightly at elevated tem- 
peratures. This oxidation was noted by a decrease 
in the resonance frequencies, as well as the usual 
change in surface appearance and increase in mass, 
after the specimens had been returned to room tem- 
erature, as compared with the original values. 
his irreversible decrease in resonance frequencies is 
superimposed on the reversible decrease associated 
with increased temperature in the course of a 
modulus-temperature determination. Table 3 gives 
typical data of such changes in mass and one of the 
resonance frequencies for some specimens before and 
after heating. 


TaBLE 3. Typical changes in mass and a resonance frequency 


before and ofter heating 





Ma Resonance 
frequency 
Specimen n ( a 
° m Ie 
Before Aft Before After 
m Neo fe fe 
g g 8 c/s 
31-1 " 71. 296 71. 561 0. 9963 2534 2515 0. 9851 
32-1 77. 205 7 979 2410 2404 9950 
33-1 75. 489 Yy¥ubY | 2312 2302 9914 
34-1 . 73. 960 999s } 2445 2444 9992 
| 


Ratio of mass before to m 
2 Square of ratios of fundar 
and after heating 


iss after heating 


nental flatwise flexural resonance frequencies before 


No significance can be given to the different 
changes for the different specimens inasmuch as they 
were not subjected to the same time-temperature 


schedules. It is seen that the largest increases in 
mass are associated with the largest decreases in 
resonance frequency. However, the decrease in 


resonance frequency cannot be accounted for com- 
pletely by the increase in mass. In order for this 
condition to obtain the following relation would have 


to hold: 
fa\? 
=( . ) 


where m is the mass; f is a particular resonance 
frequency, and the subscripts 6 and a stand for 
“before” and “after” heating. Comparison of 
columns 4 and 7 of table 3 shows that the change in 
mass accounts for only a part of the change in 
resonance frequency. If the oxide only forms a 


Ms 


(3) 


Ma 


coating over the main body of the specimen, as seems 
likely, this would account for the fact that (fa/f»)? ¥ 
Myp/Mq.* 

It was first hoped that if a specimen was heated 
once or, if necessary, several times, an oxide coating 
would eventually form which would inhibit further 
oxidation, so that subsequent determinations on such 
specimens would require no correction for this 
phenomenon. Unfortunately, this did not occur. 
Successive heatings, especially if the times involved 
were not short, always resulted in some further 
oxidation. 

Clearly, the most effective way to avoid oxidation 
would be to have the specimen in a vacuum or inert 
atmosphere during the measurements at elevated 
temperatures. An apparatus for making such 
measurements in vacuo is available and was used 
to obtain the data for the TiC specimen, in which 
the oxidation was more extensive than in the cermets. 
Wachtman and Lam [4] also used this same apparatus 
for their determination on the K-150B_ cermet 
(No. 31). However, it was found more convenient 
to use the larger furnace in open atmosphere for the 
rest of the specimens, mainly because this larger 
furnace had smaller thermal gradients, and also 
because the torsional response could be obtained 
much more readily in the larger (air) furnace. 

After a number of sets of determinations were made 
in this furnace in an air environment, the amount 
of further oxidation was considerably reduced by 
having the furnace flushed with nitrogen gas during 
a determination. The sets of determinations for 
the cermet specimens listed in table 1 include those 
made in ordinary and nitrogen flushed atmosphere. 

The temperature-frequency relations were adjusted 


for the decrease due to oxidation by raising the 
recorded frequency at the —— temperature 


reached by multiplying by the factor, fo/fa, and 
linearly interpolating this increase from the highest 
temperature down to 800 °C. Below this tempera- 
ture, it was established experimentally that no 
significant oxidation took place ia the times ordinar- 
ily involved. 

Also, since the oxide coating which a specimen took 
on on heating did not apparently affect its intrinsic 
elastic modulus nor its temperature dependence, 
the new lower values of fp after a particular heating 
could be used with the original values of Yo and G 
in subsequent sets of determinations. This procedure 
appeared justified by the reproducibility of successive 
sets of determinations when it was adopted. 

The nickel specimen did not oxidize sufficiently to 
require correction. 


3. Results 
3.1. Grain Boundary Slip 


We consider first the Young’s modulus temperature 
curves in figure 3. The general pattern is seen to 
be an approximately linear decrease until a tempera- 
ture around 800 to 1,000 °C in which region a sharper 

’ The changes inthe resonance frequencies of the specimens on oxidation suggest 
the possibility of their utilization as a means of studying oxidation as a function 


of time and temperature in addition to the conventional methods (of changes in 
mass or gas pressure), 
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decrease takes place. This inflection in the modulus- 
temperature characteristic is attributed to grain 
boundary slip (henceforth designated by the initials 
g.b.s.). This phenomenon has been discussed by 
Zener [15] using data from polycrystalline metals 
and extended by Wachtman and Lam [4] to poly- 
crystalline ceramics. 

A qualitative explanation of what is believed to 
occur may be presented in the following way: 

Until the temperature region of g.b.s. is reached, 
the grain boundaries do not contribute significantly 
to the elastic modulus, the main contribution arising 
from the elastic properties of the grains themselves. 
The “normal” linear decrease in modulus with 
temperature is due to the odd and even powered 
terms of the anhormonic expansion as described by 
Stern [16{[. (The thermal expansion is due only to 
the odd powered terms of the same expansion.) 

At some elevated temperature (which marks the 
onset of g.b.s.), viscous flow between grains adds a 
nonrecoverable, relaxed component to the purely 
elastic (unrelaxed) strains which were predominant 
up to this region. Consequently, there is a rapid 
decrease in the measured modulus, which is super- 
imposed on the linear decrease that has taken place 
up to this temperature region. This process con- 
tinues at a faster rate as the temperature rises and 
the effective viscosity decreases. Since this viscous 
movement between grains inv olves energy losses, the 
region of g.b.s. is characterized not only by a more 
rapid decrease in the measured modulus but also by 
an increase in internal friction (Q The decrease 
in amplitude and sharpness of the resonance re- 
sponses, indicative of this rise in Q', did indeed take 
place and was chiefly responsible for the loss in pre- 
cision in this region. In fact, this increase in Q™ i 
as much a feature of the region of g.b.s. as the in- 
flection in the temperature-modulus characteristic. 
Had these inflections taken place without the con- 
comitant increase in Q~', then some other mechanism 
would be required to explain them. 

As the temperature continues to rise, the effective 
intergranular viscosity also continues to decrease so 
that the energy losses associated with g.b.s. should 
also decrease. This should eventually result in a 
decrease in Q~', though not to the level existing be- 
fore the region of g.b.s. This should be evidenced 
by a partial restoration of the amplitude and sharp- 
ness of the resonance responses and also by a second 


inflection in the modulus temperature curve. This 
marks the end of the region of g.b.s. Zener [15] 


cites experimental evidence of this second inflection 
with metals. 

However, no reliable evidence of having traversed 
the region of g.b.s. was found for any of the materials 
of this investigation, nor, for that matter, with any 
ceramic materials measured by the author, even when 
the temperature was raised more than 200 °C higher 
than that for which a final reliable resonance response 
could be detected. Wachtman and Lam’s data also 
gives no evidence of having traversed the region of 
g.b.s. Two possible explanations which suggest 
themselves for this are that the region of g.b.s. 
covers a much wider temperature region than antic- 








ipated, or that, if the region is traversed, the im- 
proved response is not sufficient to be detected above 
the higher general noise level associated with the 
various mechanisms which are more likely at ele- 
vated temperatures. 

The upper temperature limit of recorded elastic 
moduli was reached, then, when the responses be- 
came so diffuse and weak that they could not be 
distinguished from the general background noise and 
spurious responses which inevitably develop when it 
becomes necessary to drive the specimens harder 
and harder as the true responses become progressively 
weaker. For the nickel specimen this limit occurred 
when the modulus temperature relationship was still 
linear, as it did mainly for the torsional responses 
of the cermet specimens, shown in figure 4. 

The general characteristics of the modulus tem- 
perature curves as well as the recovery of sharpness 
and amplitude of the resonance responses are found 
to be reversible on cooling. Also, in the absence of 
oxidation, the specimens generally return to the 
same room temperature values of resonance fre- 
quency that they had originally, indicating no per- 
manent average structural change due to g.b.s. 

It is also noted from figure 2 that no noticeable 
corresponding inflection occurs in the thermal ex- 
pansion curves around the temperature region in 
which the inflection associated with g.b.s. occurs in 
figure 3. 

We now compare the types of specimens with 
respect to g.b.s. It is noted that for TiC the inflec- 
tion indicating the start of this region comes at a 
higher temperature than for the cermets; and, once 
started, the modulus drops more sharply for TiC 
than for the cermets. The presence of the metal in 
the cermet has the effect of hastening the start of, 
and broadening the measurable region of g.b.s. The 
contribution of the nickel on this region in the 
cermets is difficult to evaluate since the responses 
for the nickel itself died out before g.b.s. developed. 
However, it seems reasonable that the presence of 
more than one type of grain boundary should result 
in a broadening of the g.b.s. region. The presence 
of two types of grains such as are present in the 
cermets results in three types of grain boundaries, 
Ni-Ni, carbide-carbide, and Ni-carbide. 

In general, if there are n types of grains, and it is 
assumed that each grain forms a boundary either 
with another grain of the same type or with a differ- 
ent type of grain, then the number of types of grain 
boundaries will be given by n(n+1)/2. 


3.2. Temperature Coefficient of Elastic Moduli 


The relative decrease with temperature in Young’s 
and shear modulus in the linear region, —1/Y)(AY/At) 
and —1/G,(AG/At), and their standard deviations are 
given in table 4. These values were computed from 
a least square solution of the numerical data on 
which figures 3 and 4 are based. The range of Af 
selected (and also given in the table) was as wide as 
possible to give representative values without enter- 
ing into the inflection region marking the onset of 


g.b.s. 


—_ ————— 





PaBLe 4. Relative decrease in Young’s and shear moduli over 
the linear ranqe 
l ( \} ) l ( AG ) 
Specimer : F 
—— ) \ ‘ ( G M At bin © ¢ 
iT q ( 10-5 per © ¢ 
ric 10. 4+0.3 On) 
31 12.4 l 900 13. 240.3 900 
32 13. 2 SOO 14. 94 S 690 
$3 14.4 2 700 15. 94 3 675 
34 14 ] SUK) 15. 64 1 SO) 
1 35. 74 t 600 
Relative decrease 1 Your modulu 
lemperature range from room temperature valu iven 
Relative decrease in shear modulu 
Standard deviation 


It isseen from table 4 that both —1/¥,)(AY/At) and 

1/G)(AG/At) increase in absolute value with in- 
creasing metal content. TiC and Ni (for which only 

1/¥)(AY/At) are available) have the lowest and 
highest values of relative slope respectively, forming 
the “end members” in the foregoing pattern. 

Cermets No. 33 and 34, which have the same per- 
cent of metal content and vary from each other only 
in that one (No. 33) contains only nickel as the metal 
while the other (No. 34) contains 25 percent nickel 
and 5 percent molybdenum, show very little differ- 
ence in 1/¥)(AY/At) and 1/G)(AG/At). From 
figure 2. il may also be seen that the specimens 
arrange themselves in the same order with respect 
to positive coefficients of thermal expansion, 

/,(A//At) as they do for the negative temperature 
coefficients of elastic moduli; 1.e., 
for both these coefficients regardless of sign ure 
associated with nickel and decrease in the same order 
for the remainder of the specimens. 


the st eepest slopes 


3.3. Titanium Carbide and Nickel 


The value of }Y, obtained for the Tic specimen 
(given in table 3) is not definitive for TiC in general. 
First, the specimen is of commercial purity and small 
impurities present probably affect. the measured 
value of Y,) to some extent. Second, in the process 
of manufacture, TiC is often not produced in the 
proper stochiometric ratio [17]. This can also affect 
the value of Yy. The third and most important 
factor affecting the value of Y) of the TiC specimen 
was that it was quite porous. Its density, as derived 
from the mass and volume, computed from its dimen- 
sions was 3.558 gm/cm Using the theoretical 
density of 4.938 gm/cm® of Schwartzkopf [17] one 
calculates the porosity of the specimen to be 28 
percent. 

Mackenzie [18] has developed a theoretical rela- 
tion showing how the elastic modulus increases as 
the porosity Coble and Kingery [19] 
have obtained experimental results in good agreement 
with Mackenzie’s relation using Al,O; specimens of 
different porosite 8. Using Mackenzie’s relation, the 
Y, of the TiC spec imen was estimated to be 3411 
kilobars at theoretical density. One should also be 
cautioned here that the ceramic phase of the cermets 
is not “pure” TiC but a solid solution consisting of 
TiC and other carbides (sec. 2, 2.1, a) 


decreases. 


Coble and Kingery also showed in the paper just 
referred to that the relative slones of the Young’s 
modulus temperature curves, —1/¥,(AY/Ad#), in the 
linear region were essentially independent of the poros- 
itv of the specimens. If it is assumed that this also 
holds for TiC, then the value of —1/¥,(A¥Y/At) given 
in table 4 will not be affected. Incidentally, the 
relative changes in the elastic moduli may be com- 
puted directly from the primary data, Y,/¥> or 
G,/G,, and are independent of the particular values 


selected for ), or G). 
The value of 1/Y,(AY/At) obtained here for 
nickel (3.6 107* per °C) was in fairly good agree- 


ment, as far as could be determined from the graphs, 
with that obtained by Koster [19] for a cold worked 
unmagnetized specimen (3.8 107* per °C). Koster’s 
range of Af was from 180 to 400 °C, but At selected 
for the value just given was from 20 to 400 °C 
Further data taken with this specimen were also in 
veneral agreement with those obtained by Koster for 
annealed specimens. 


3.4. Comparison of Relative Slopes of Young's and 
Shear Modulus 


Inspection ol table 4 also shows for each type of 
specimen (W here both values exist) that 


1 (AG l fAY » 
a. \ ae ) AGT ) 9 
This means that Poisson’s ratio, M, increases with 


increasing temperature. 

Individual values of uw calculated, from the nu- 
merical values of Y, and G@, in all cases where both 
values existed, from eqn (2) confirmed this small but 
significant increase in w with increasing temperature 

Coble and Kingery [19] also found Poisson’s ratio 
to increase with temperature for Al,O; as did Lang 
and Knudsen [20] using ThO,+-0.5 percent CaO by 
weight. 

Zener [15] has deduced from theoretical considera- 
tions that the relaxation associated with g.b.s. should 
cause the relative decrease in the measured value of 
G to be three times that of Y. Unfortunately this 
result cannot be checked experimentally, not only 
because this region could not be observed for G, but 
also because the upper temperature region of g.b.s 
could not be determined for Y. 


4. Summary 


l. The temperature dependence of Young’s mod- 
ulus and shear modulus on four types of nickel- 
bonded TiC specimens have been determined by a 
dynamic resonance method. Similar measurements 


have also been made on nickel and TiC for 
comparison 
») 


2. The relative decrease in both Young’s and shear 
modulus with temperature was in the following order 
with the material having the smallest relative de- 
crease listed first: TiC, 10 percent nickel content 
cermet, 20 percent nickel, 30 percent nickel, and 25 
percent nickel +-5 percent molybdenum, and nickel. 





| 


| 
| 
| 


3. The inflection in the Young’s modulus tempera- 
ture curves was attributed to grain boundary slip. 

4. The relative decrease in shear modulus with 
temperature was greater than that for Young’s 
modulus, where both values were known, for the 
sume type of specimen. ‘Therefore Poisson’s ratio 
rises slightly with temperature. 


The author is indebted to J. Pastine for perform- 
ing a portion of the calculations on which figures 3 
and 4 were based, and to H. S. Parker and S. M. 
Lang for making the thermal expansion measurements 
on the cermet specimens. 
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Analog Simulation of Zone Melting 
H. L. Mason 


(January 13, 1961) 


An analog computer was used to simulate zone melting in small-diameter tungsten rods 


of 0.2 to 1 centimeter. 


for the chosen spatial subdivision 


Heat balance equations take account of electron bombardment, heat 
capacity, radial conduction, axial conduction, and surface radiation. 
and for 


Diagrams are given 


the computer mechanization. Predicted 


temperature-time curves are shown for a rod of 0.4-centimeter diameter with optimum 


power input of 650 watts. 
1-centimeter rod. 


* As part of the NBS program for attempting to 
supply metal samples so pure that no spurious lines 
can be detected by spectroscopy, G. A. Moore of 
the Chemical Metallurgy Section has developed 
an apparatus for zone melting of small-diameter 
tungsten rods. Because of the high cost of the semi- 
purified material used in the process, it became 
desirable to simulate the heat transfer operation as 
a means of determining optimum relations between 
rod diameter, heating rate, and progress of the molten 
zone. 

The technique is applicable 
and the problem is interesting from the simulation 
viewpoint for several reasons. The conventional 
diffusion equation frequently used for analog com- 
putation is not suitable because it requires precise 
knowledge of the temperature distribution at the 
surface. A model must be set up in terms of heat 
flows; these too are unknown, but can more readily 
be estimated from the physical situation. The 
model must include radiation losses, which at these 
high temperatures require fourth-power nonlinear 
terms. ‘The change of state at the melting point 
results in a local interruption of the temperature 
rise in the melting zone. Conductivity, thermal 
capacity, and emissivity will change with temperature 
and with state. Time rates of temperature change 
will be treated as infinitesimals, spatial rates as 
discrete increments. The incremental dimensions 
for radial and axial subdivision must be carefully 
chosen if optimum use is to be made of the limited 
number of integrating amplifiers available. 

The radii 7, of the 5 rods to be simulated range from 
0.1 to 0.5 em; lengths (actually 3 or 4 ft) are assumed 
infinite. Heating is by continuous electron bom- 
bardment in a vacuum from a heated wire of 1.2 mm 
diameter, shaped to an annulus of 4.0-cm diameter 
surrounding the rod as shown in figure 1. The 
heater is initially fixed near one end, but after 
melting begins it is moved very slowly up (or down) 
the rod. ‘The impurities are retained in the molten 
zone and are carried to one end, leaving the purer 
metal to solidify again to rod form. This motion 
neglected in the simulation. Circumferential 
symmetry and an isomorphic rod are assumed, and 
the heat input from the wire is taken to be a step 
function of time. 


to other materials, 


is 
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Unsatisfactory operation is predicted for zone melting of a 


Because of symmetry on either side of the heater 
midplane, it is enough to consider only half the length 
of the rod. Preliminary estimates indicated that 
temperatures would rise only slightly at points 20 mm 
distant from that plane, so that heat flow beyond this 
was neglected. Steep temperature gradients may be 
expected near the heater, and it is necessary to assure 
that metal will not drip off the outer surface before 
the core is melted. One purpose of the simulation, 
therefore, was to find whether the molten zone ad- 
vanced radially at least as fast as it advanced axially. 
Subdivisions into annuli or cylinders with lengths!Ay, 
outer radii 7,, and inner radii r; were decided upon. 
Their generalized proportions relative to surface radii 
r, are given in table 1. 





POWER 


INPUT — 





HEATER WIRE 





FicurRE 1. Assumed subdivision of rod. 








The heat stored per unit time in zone 1, which is 
assumed to have a uniform temperature of 7, degrees 
Kelvin, at time @ seconds after heating begins, is the 
thermal capacity of the zone volume multiplied by 
the rate of temperature rise. With density yg/em* 
and temperature-dependent specific heat e=c(7;) 
cal/g-deg K, this storage rate is yer(r?—7r?) Ay dT;/dé. 
Similar expressions hold for the other zones, with r, 
replacing 7, and 7;=0 for the inner ones. 


TABLE 1. 


Vim 


t 


Heat inflow to zone 1 will be 239 P eal/see, where P 
is the steady kilowatt power input to the half-rod. 
Heat outflow by reason of conductivity k=k(7)) cal 


(sec) (em)? (deg K/em) will be radially inward to 
zone 2 across the area 2ar,;Ay with the gradient (7; 
T,)/(m,—1r2), where 7; and rz are mean radii. Also, 


there will be conduction axi: ally to zone 3, at the rate 


r?) (T,—T. 


rk(r? Yi— Ys) where y, and y; are Mean 
lengths. No heat is transferred bevond the heater 
midplane, since there is no gradient across that 


boundary. However, zone 1 loses heat by radiation 
to ambient temperature 7, at the rate 1.37¢ (2a7r,Ay) 
[(7,/1000)*—(T,/1000)*| where e=e(7)) is the tem- 
perature-dependent emissivity. The heat flow bal- 
ance equation for zone 1, at any time up to the 


instant of local melting, is then 
T, 
V2 


T; 
239 P—2rr kay | 
r7; 
z 
) 
Yi 
3 i 
—1.37¢€(27r,Ay) | ( ) —( 


Similar heat flow balance equations can be written 
for the — zones except that 239P is replaced by 
radial or axial conductive inputs, and there is no 
scalianines for the interior zones where rp replaces 7, and 
r,=—U. 

For zone 3, 


dT, 


ray 10 


VCR Ts: 


—ark(y 


T; 


1000 1000 


ry, I 


for example, the equation is 


—2rr,kAy E 2 ] 
r3—T"4 
teary [BE ater [BEE 
Yi—Y Ys 
a al rl f - 


T 
—1.37€( eres) | (7535) ~ 1000 


dT, 
r;)Ay dé 
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Numerical values of the properties y, ¢, &, and 
were known only imprecisely for tungsten, but it ap- 
peared that they change only slowly in the range 
from T,=300 to 7;,,=3700, the melting point. Ac- 
cordingly, the values of table 2 were used in the com- 
putation, based on preliminary estimates of temper- 
atures with zone 1 melted and zones 2 and just 
beginning to melt. The heat of fusion in zone 1, 
about 44 cal/g, was neglected in this approximation. 
For each zone, the initial temperature was assumed 
to be 300 °K, and regions bevond zone 9 were as- 
sumed to remain at this temperature, 

Substitution of numerical values in the 9 expres- 
sions for heat flow gave for the rod of l-em diameter 
the following set of simultaneous first-order differ- 
ence-differential equations for solution by the analog 
computer: 


-0.27,=.015P—.0394 T-T; 569( 7\—T,) 
— .0498T4/10' 
oT. .118(7;-T;) + .363(7,-T,) 
0.27. .369(7,—T;) +-.0388(7;-T,) +.173(T7,-T;) 
0498 74/10! 
0.27, 406(7.-T,) + .1302(T,-T,) —.23(T,-T,) 
0.27; 39(7T;-T. .3195(T;-T,) —.261(7;-T;) 
134 T2/10' 
27, .115(7T.-T, .2116(7;-T,) + .067(7,-7;) 
T; .129(7;-T;) —.0504(7,—-T;) +-.167(T,-T, 
T 167(7;-T.) +.167(T; Ts) 
a 167(7,;-T») —.167 7» 
TABLE 2 
Zon I € 
3700 0. 84 0. 31 0. 36 
3000 73 30 ‘a 
11K a4 13 
729 sD 0 ro 0. 


Figure 2 shows the computer patchboard diagram 
for the first two and last two of these equations. 
The time scale was slowed in the ratio of 1:20, and 
1 machine volt represented 100° Kelvin. Heater 
input P was adjusted by trial with the aim of estab- 
lishing an optimum value which would use a mini- 
mum of energy without causing zone 3 to melt 
ahead zone 2. For the 1l-em rod, such a value 
could not be found, since zone 3 always reached 
3,700° before zone 2, while the temperature of zone 
| went up to about 4,600°. This represents an 
undesirable physical situation which in practice may 
prevent satisfactory processing of a rod of this large 
diameter. In view of the uncertainty as to actual 
values of 7, c, #, and €, it was considered unnecessary 
to carry through the computation for a more refined 
approximation. 
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A set of temperature-time curves for one of the 
smatler with 650-watt input shown in 
figure 3. It will be noted that the melting of zone 
3 lags that zones 1 and and hence that the 
power can be adjusted so that zone 3 just fails to 
melt. This assures formation of a thin liquid laver 
which ts phy sically stable. The rod is nearly melted 
through in about 3 sec. This finding required modi- 
fication of the power control circuits of the experi- 


rods ul 
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Zone 


FIGURE 3. melling 2-mm-radius tungsten rod. 


mental apparatus by addition of rapid-response 
components to supplement the original  servo- 
mechanism, which had required about 380 see to 
establish an equilibrium setting. 

Assistance on this problem came from NBS 


colleagues G. A. Moore, who supplied the physical 
constants and offered much constructive advice; P. 
K. Wong, who operated the computer; and R. 
Shives, who plotted the curves. 
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Residual Losses in a Guard-Ring Micrometer-Electrode 
Holder for Solid-Disk Dielectric Specimens 


Arnold H. Scott and William P. Harris 


(January 19, 1961) 


A guard-ring micrometer-electrode holder especially designed for use in making dielectric 
measurements on solid-disk specimens of low-loss materials has been constructed. A 


method of determining the separation of the electrodes to 


gages, is described. 


1 micron, using ball reference 


Residual losses in the guard-ring holder, and also in the bridge standard 


capacitor used in the Schering bridge employed for the measurements, caused by losses in 
surface films on the electrodes and by series resistance in the leads, were accurately measured 
It is shown that at frequencies from 100 cycles per second to 100 kiloeycles per second 


these residual losses are not more than a few microradians. 


Using corrections thus obtained, 


the holder-bridge combination can be used to measure the losses of low-loss materials with 


greater accuracy than heretofore. 


1. Introduction 


This paper describes a guard-ring micrometer- 
electrode bolder and methods by which the separa- 
tion of the electrodes, the capacitance, and the 
residual losses of the holder can be precisely deter- 
mined. The holder is one which was designed and 
constructed to meet the need for a very accurate 
method of measuring the dielectric properties of 
some standard reference specimens of dielectrics. 

Such standard reference specimens [1]! are useful 
in checking the various techniques and equipment 
used for measurements of the dielectric constant 
and dissipation factor of insulating materials. Two 
general classes of dielectric measurements are used: 
two terminal and three terminal. The three ter- 
minal method is usually more accurate, but may 
require more time, especially if a three-terminal 
bridge with an earthing circuit is used. For routine 
measurements in quantity, the two terminal method 
is often used. This method generally requires cor- 
rections for the edge capacitance and ground capaci- 
tance of the specimen under test. Standard dielectric 
reference materials permit a check on the accuracy 
of such measurements. 

Recently there have come into general use two- 
terminal holders which eliminate most of the edge 
capacitance and all of the ground capacitance from 
the measurement, and permit relatively rapid meas- 
urements to be made. However, the loss measure- 
ments include the loss over the edge of the specimen 
and the loss through the insulation of the high- 
potential electrode of the holder. Accurate deter- 
mination of the electrode spacing is also required. 
Here again, standard dielectric reference specimens 
are useful in checking the techniques and measure- 
ments involved. 

In order for standard reference specimens to be 
useful, it is necessary that their properties be known 


1 Figures in brackets indicate tho literature references at the end of}this paper. 
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accurately, apart from the methods and equipment 
used to measure these properties. This requires 
that residual losses and capacitances in the specimen 
holder be reduced to negligible proportions, or if this 
is not feasible, that they be accurately determined 
and separated from the measurements made on the 
specimens. <A further complication arises from the 
fact that most materials that might be used for 
standard reference specimens exhibit properties that 
are affected by the ambient humidity, and must 
therefore be studied under different conditions of 
humidity in order that this effect might be evaluated, 
and that the dielectric properties under a given set 
of conditions might be known. ‘This involves the 
determination of the variation of the residuals 
as a function of humidity. Of course, temperature 
dependence must also be considered. 

The foregoing considerations prompted the con- 
struction of the specimen holder to be described, and 
led to the measurements made to determine the 
residual losses under various conditions. Methods 
are described which permit the determination of the 
losses in the holder, arising from series resistance in 
the connecting leads, from losses inherent in any film 
that may exist on the surfaces of the electrodes, or 
from any other source. From this, corrections can 
be applied to measurements made on specimens 
placed in the holder. 

A useful bonus of the method used was that the 
losses in the standard bridge capacitor used were 
also determined, and a correction for this source of 
error was obtained. 


2. Description of the Specimen Holder 


The holder consists of a heavy cylindrical case 
with the electrodes mounted coaxially inside. The 
central guide rod of the movable electrode extends 
through the top of the case and its position is con- 
trolled by a thin nut located just above the top of 
the case. The position of the movable electrode is 








sensed by a micrometer which comes in contact with 
the top of the guide rod. Figure 1 is a cross sectional 
diagram of the holder, 
installed in a controlled-humidity chamber. 

The guarded electrode and guard-ring are mounted 
on a ceramic disk. This disk was ground to nearly 
optical flatness and fits into a machined recess in the 
base of the case structure. A hole was ground in 
the center of this ceramic disk so that the centering 
rod of the guarded electrode would fit snugly into 
this hole thus holding the guarded electrode ac- 
curately centered on the ceramic disk. The guarded 
electrode is held snugly against the ceramic disk by 
means of screws in recessed holes in the ceramic disk. 
The guard-ring was recessed to fit snugly over the 
outside edge of the ceramic disk so that the guard- 
ring is also accurately centered on the ceramic disk. 
This permits the maintenance of a uniform gap 
between the guarded electrode and the guard-ring. 
Screws through the base plate and the ceramic disk 
into the guard-ring hold this electrode assembly 
rigidly against the base plate and make the necessary 
electrical connection between the guard-ring and the 
surrounding assembly. The surface of 
electrode assembly was ground optically flat. 

The upper movable electrode is also mounted on 
a ceramic disk whose surfaces have been accurately 
ground. This ceramic disk is in turn mounted on 
the carefully machined surface of a metal disk 
attached to a guide rod which extends through the 
top of the barre! assembly. This guide rod passes 
through two guide bearings about 7 em apart. A 
metal bellows is attached to the mounting disk and 
to a collar which fits over the guide bearing 
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Figure 1. Schematic diagram of the gquard-ring micrometer- 
electrode holder for solid-disk dielectric specimens. 
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and figure 2 shows the holder 
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FIGURE 2. 


in place in controlled- 


atmos phe 


A set screw in the collar permits the location of the 
collar so that the bellows can be made to apply a 
force tending to move the movable electrode toward 
the guard-ring assembly. The collar and bellows 
also prevent the rotation of the movable electrode 
assembly. 

The upper end of the guide rod is threaded, and a 
large diameter knurled nut is used on this to set the 
spacing between the This nut works 
against the top of the guide bearing assembly. Fine 
threads (32 threads per inch) are used so that with 
the spacing between the 
within 1 x. 

The spacing between the electrodes is determined 
by means of a micrometer which is mounted on a 
frame so that it is held directly above the end of the 
guide rod. A steel ball is mounted on the end of the 
guide rod and the position of the guide rod is deter- 
mined by serewing the micrometer down until it 
just touches this ball 

Contact with the ball can be accurately determined 
by a very light touch. The micrometer is equipped 
with a large graduated wheel so that the smallest divi- 
sions are equivale nt to 2 yp. With a little practice 
and skill the micrometer measurements can be made 
repeatedly within a micron, by different observers. 

The top and bottom of the cylindrical case can be 
removed to expose the electrode surfaces. The ends 
of the evlindrical ease are carefully machined to be 
accurately perpendicular to the axis of the cylinder. 


electrodes. 


eare electrodes can be set 


The end plates and electrode assemblies are carefully 
ground and machined so that when they are placed 
on the end of the eylindrical case the electrode sur- 
faces are also perpendicular to the axis of the ceyl- 
inder. Only slight adjustments are necessary to 
make the electrode surfaces parallel to the accuracy 
which is required in this work. 

The electrodes were made of A.I.S.I. Type 304 
stainless steel. No plating or coating was applied. 
The surfaces were ground and lapped smooth. 

Connection is made- to the gui arded electrode by 
means of a coaxial connector whose outer case is 
attached to the base plate of the holder. Connec- 
tion is made to the upper electrode by a flexible wire 
which is clamped under one of the mounting screws 
which hold the electrode to the ceramic disk. The 
other end is connected to the core of the coaxial cable 
which passes through the wall of the evlindrical case. 

Two windows about 7-em wide and 3-em high are 
cut in opposite sides of the cylindrical case at the 
electrode level to permit the introduction and _ re- 
moval of samples from between the electrodes. A 
sliding band around the 
to cover these 
required, 


evlindrical case can be used 
openings if complete shielding is 


3. Determination of the Area of the Guarded 
Electrode 


To determine the area of the guarded electrode 
the base plate and cuarded electrode assembly were 
removed from the evlindrical case. The diameter of 
the guarded electrode, and the gap between the 
guarded electrode and the guard-ring were measured 
by means of a traveling microscope. The effective 
area of the guarded electrode was then computed 
using the gap corrections published in ASTM D-150 
[2]. The effective area was found to be 20.366 em?. 


4. Determination of the True Spacing of the 
Electrodes 


In order to adjust the upper and lower electrodes 

be parallel and to measure the actual separation 
of the electrodes, a set of ball reference gages was 
prepared and accurately measured, Steel balls with 
nominal diameters of 39, 's, °49, 316, 72, and 4 in. were 
rigidly mounted in strips of plastic. These strips 
served as handles for placing the gages at the desired 
points between the electrodes, and the rigid mounting 
insured that a definite diameter was always presented. 
These reference gages were measured by the Engi- 
neering Metrology Section of NBS, and their diam- 
eters determined under no-load condition with an 
accuracy of + 0.00005 em. 

Parallelism of the electrodes was established by 
using a ball mounted similarly to those described, 
but not necessarily accurately measured, to complete 
the electrical circuit of an ohmmeter connected to 
the two electrodes. The electrodes were brought 
together by the thumb nut to such a distance that 
the ball just made contact between the electrode 
surfaces, but vet was free enough to move about. 





| different 


If there was any lack of parallelism, this was made 
readily apparent by the existence of areas on one 
side of the electrodes where contact was not made 
while at the opposite side contact was made, or 
possibly the ball would not slide easily between the 
electrodes. Adjustments were made by tightening 
or loosening the appropriate screws located around 
the circumference of the case. 

It was found that after such adjustments had been 
made, their full effect was not immediately seen. It 
took several hours for the dimensions to stabilize. 
This is attributed to two effects. First, the act of 
making the adjustment requires some handling of the 
cell, with some disturbance of the temperatures of 
the various parts of the cell case. The other factor, 
and perhaps the major contributor to the long time 
required for dimensional stability to be reestablished, 
is that the metal probably creeps for a time after the 
stresses are altered. The result of these effects is 
that the final adjustments must be made by making 
very small changes and then waiting several hours 
or overnight for the ensuing changes in the metal 
take place. 

This same ball-gage technique can be 
determine whether the guard ring and the 
electrode are in the same plane, 
electrode as a reference surface. 

After the electrodes had been made pare allel, a 
comp: irison was made between the measured capac i- 
tance and that computed for the seperation of the 
electrodes as determined by ball reference 
diameters. Typical results are Bieri in 
table 1. It was coneluded, from the fact that 


used to 
cuarded 
using the other 


rages ot 


PABLE 1. Comparison of measured capacitance with the capac- 


ilance computed from the electrode separation 


S pac y Computed C) Measured C AC 
‘ 
1 ] pf 
0. 31773 5. 6754 5. 6755 0). OO 
17638 $3. 7853 3. 7856 O003 
558 3. 2457 3. 2456 On 


that the reference 
now be used for occasional check meas- 
urements, the measured capacitance used to determine 
the electrode separation with the holder empty, and 
these measurements used to find the zero correction 
to be applied to the micrometer readings. It is 
desirable to have a rapid means of determining this 
micrometer zero correction, because any adjustments 
of the micrometer mount, or of the electrodes, such 
as resetting them to be parallel, changes the microm- 
eter zero correction.2 It was found that the zero 
correction can be determined readily and repeatably 
from capacitance measurements. Some results of 
such measurements are shown in table 2. 

The agreement of the micrometer zero correction 
determined for the various spacings is acceptable, 


excellent agreement was obtained, 
gages could 


2 A significant change of temperature ers ilso change the correction, For 
this work, temperature is held constant to +1 degree C. The effect of this small 
temperature change was found to be negligib x 
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TABLE 2 


2. Determination of micrometer zero correction from 
measured capacitance 
Micrometer | Capacitance | Computed Difference 
reading spacing 
em pi cm cm 
0. 2499 7.177 0. 2513 0. 0014 
1194 14. 926 . 1208 . 0014 
0793 22. 324 . 0808 OO15 
0590 29. 847 . 0604 0014 
. 0471 37. 159 0485 0014 
0390 44. 566 0405 0015 
0332 52. 077 0346 . 0014 


and following these measurements a zero correction 
of 14 w was used until further adjustments were 
made to the electrodes. 

It is estimated that the separation of the electrodes 
can be determined to about +1 yu, using the microm- 
eter in this fashion. 


5. Losses in the Specimen Holder 
Standard Bridge Capacitor 


5.1. 


and 


General Considerations 


It was intended that this holder should be used 
during accurate measurements of the losses in 
materials that have rather small losses. Therefore it 
was necessary to evaluate the holder for the effects 
of inherent losses. It is evident that any series 
resistance in the leads connecting the holder into the 
measuring circuit will appear in the measurements 
as a loss. Also, losses in any surface films on the 
electrodes must be evaluated. It is known that 
some electrode materials form oxide films which 
might be expected to exhibit losses [3, 4, 5, 6, 7, 8, 9, 
19, 13]. 


5.2. The Bridge 


A low-voltage Conjugate Schering Bridge was 
used to make the measurements. The substitution 
method was employed. A diagram of the bridge 
circuit is shown in figure 3. Because the residual 
losses in the variable bridge standard capacitor C, 
are involved in the measurements, it is necessary to 
evaluate these This requirement was no 
appreciable handicap, because the method employed 
yielded simultaneously the losses in both the speci- 
men holder and the standard bridge capacitor. 

Measurements by the substitution method are 
made by setting Cy to some value of capacitance by 
adjusting the separation of the electrodes, and setting 
C,, the standard bridge capacitor, to some initial 
value, usually 20.000 pf. The capacitor C; is also 
set to some convenient initial value. After setting 
Sw 1 to position A, the bridge is balanced by varying 
C1, Cy, Cs, and R;, while alternately setting switch 2 
to the A and B positions, until the detector shows a 
null in both A and B positions. The values of (, 
and ©; are noted, and will be designated C} and C3, 
respectively. Sw 1 is now moved to position B. 
Leaving ©, and C, unchanged, the bridge is rebalanced 
by varying ©,, C;, C;, and R;. The new values of 
C, and C; are noted, and will be designated as C’; and 


losses. 















C, Se 
—Te eee ree, See : 
== = == so ae Por sae <4 
ol * ijt 4 aa! 

' if! rot yis™! ae 
' ' si - ' 1 
(|! Saepera. 7, 
DETECTOR - ‘fa CH A! 
| ‘he \ -i+, (SPECIMEN HOLDER) __ |! 
_— ~ ' 
| — rB ' 
Ra ~ 
° ee ee, eee | 
ee 2 Sin «eo ar > = ; 
ele ae « offs, it 
4 i= 
' Ce 

















Figure 3. Schematic diagram of the conjugate Schering bridge 


used for this investigation. 


C";. The difference between the second and first 
values of (©,, AC), is the capacitance of the holder,’ 
Cy. The difference between the second and first 
values of C;, AC; is used to compute the value of the 
dissipation factor of the holder, Dy, according to the 
equation * 


wR CAC, , GI—G, 
Deed, (5.1) 
AC, wAC’, 
where 
D,y,z=dissipation factor (tan 6) of the specimen 


holder ; 
R;,=value of FR, figure 3, in ohms; 
C"J=capacitance of C, with Sw 1 in position B, 
which is the total capacitance in the meas- 
uring arm of the bridge; 


AC,=C", —C;=difference between final and initial 
values of (C;; 
AC,= C7 —C)= difference between final and initial 
values of C); 
G,=equivalent parallel conductance of C, cor- 
responding to the initial value C}; and 
G’; =equivalent parallel conductance of C, cor- 


responding to O’;. 


The first term on the right in eq (5.1), wh3C,’’AC;/ 
AC, is the dissipation factor computed without re- 


> This statement is true for the case at hand, where the correction terms in the 
exact expression are negligible. The exact form is given inthe appendix, section 7. 





* This equation is derived in the appendix, section 7, and appears there as 
eq (7.7) 
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gard to the second term, which is a correction term 
taking into account the effect of non-zero and non- 
constant dissipation factor, 2, in the bridge standard 
capacitor, C,. Let us then write 


Dy=Dy+ Dr, (5.2) 
where 
Dy=woR,Cy AC;/AC,=the measured or appar- 
ent dissipation factor 
as given by the uncor- 
rected bridge measure- 
ment. 
and 
D,.=(G) —G))/wAC,=the correction term to be 


applied to the apparent 
dissipation factor to 
take into account the 
dissipation factor of the 
bridge standard capac- 
itor, C). 


In order to consider the effect of elements that be- 
have as series resistances, it will be more convenient 
to have eq (5.2) written in its series equivalent form, 
neglecting for the moment elements that behave as 
parallel conductances. If we separate the series 
resistance associated with the bridge standard capaci- 
tor, C;, into two parts, one part corresponding to the 
series ohmic resistance of the leads to C,, the other 
part corresponding to the losses caused by a surface 
film on the electrodes of C,, eq (5.2) becomes ° 


Dy Du T QwhhosrC > T wos Al p T wO Rog’ (5.3) 
where 
Ros, =the equivalent series resistance of the 
ohmic resistance of the leads, connec- 
tions, ete., of C,, 
and 
Rosp=the equivalent series resistance of the 


film on the surfaces of the electrodes 
or C;. 


5.3. Consequences of Equation 5.3 


Examination of eq (5.3) reveals that there are two 
correction terms involving /2s,;, one dependent on 
C,', the initial value of the bridge standard capacitor, 
the other dependent on AC, There is one term 
involving Page. 

Let us now consider the behavior of these correc- 
tion terms as a function of the capacitance of C2, 
and examine their effect on the value of Dy, the un- 
corrected or measured value of Dy. In the following 
arguments, it will be assumed that all measurements 
are made at a fixed frequency, except as stated. 
Also, for any single set of curves, it will be assumed 
that C) has a single value, usually, in the actual 
measurements to be presented later, 20.000 pf. 


§ This relation is derived in full in the appendix, section 7. 
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Taking first the simplest case, involving only the 
term due to losses in a surface film, assuming for the 
moment that Rs, is negligible, we see that the 
correction term is 

wRosp’ Ci= Der (5.4) 
This term is a constant, even for different values of 
©}, because R.sp’ is inversely proportional to the 
capacitance,® in this case C). Therefore we may 
plot this term versus capacitance as a straight line 
of zero slope, as shown in figure 4. 

If we now consider the effect of a lossy surface 
film on the holder electrodes, where the area is 
constant and the capacitance is a function of the 
separation of the electrodes, we see that a dissipation 
factor proportional to the capacitance results.’ Let 
us designate this component of dissipation factor, 
due to a lossy surface film in the holder, as 

Dyp=oRysrC y= oh ysrAC2, (5.5) 
where Ryspy is the equivalent series resistance cor- 
responding to the surface film losses in the holder. 
A plot of Dy versus capacitance would be a straight 
line with positive slope, passing through D=0, 
AC,=0 (or nearly so, depending on the ground 
capacitances [11, 12]). This is shown in figure 4, 

From eq (5.2), we see than Dy=D,—D,, or in 
this case Dy=Dyr— Dy. Dy is plotted in figure 4, 
where it is seen that it is a straight line with the 
same slope as the plot of Dyr, and with a negative 
intercept equal in magnitude to the constant value 
of D.y, if the previously assumed conditions exist. 




















Due 

Dor 

D pe 

0 —_ - —— 

1) 20 40 60 80 100 Co 
0 20 40 60 80 AC, 
Ficure 4. A plot of assumed, arbitrary values of dissipation 
factor due to a lossy surface film on the electrodes of the 
specimen holder, Dyry, and those of the bridge standard 


capacitor, Dor. 


The dissipation factor that would be indicated by the bridge under these con- 
ditions, Du= Dar—Dzr, is also plotted. These are plotted against AC) (see 
text), assuming an initial value of C2 of 20.000 pf. 


6 It seems reasonable to assume that a surface film would be uniformly dis- 
tributed over the surface. Thus, the conductance of the film would be propor- 
tional to the active electrode area. The capacitance is also proportional to the 
active electrode area. Therefore the resistance Rosr is inversely proportional to 
the capacitance. 

Near the edges of the capacitor places, where fringing occurs, this proportion- 
ality would not hold exactly. However, results show that the nonlinearity is 
negligible in the present case, as might be expected from the geometry of the 
capacitor used. 

? This assumes that the equivalent series resistance corresponding to the surface 
film losses is independent of the electric field. Tests at various voltages applied 
to the bridge yielded the same results, thus making this assumption appear 
reasonable for this case. 








Let us now consider the effect of the series 
resistance of the leads, switches, or other circuit 
components used to connect the specimen holder 
into the measuring circuit. A component of dissipa- 
tion factor, which we shall designate Dy,, would 
result. This dissipation factor would behave in 
the same manner as Dy,y, because it is of the same 
form: 

Daur ORyspAC>. (5.6) 
Figure 5 shows a plot of Dy,, Dur, and the total 
dissipation factor of the holder Dyr. 

Dyr=Darit+ Daur. (5.7) 
Assuming the same dissipation factor for the bridge 
standard capacitor as in the previous case, D,,, 
and plotting Dy—Dyr—Dy, we have the results 
plotted in figure 5. We see that the slope of the 
Dy, plot is the same as that of Dyr, with negative 
intercept equal in magnitude to Dy. 

Let us now consider the effect of series resistance 
in the leads connecting the bridge standard capacitor 
C, to the measuring circuit. From eq 
see that there are two terms involving this resistance : 
QwRospCs and wlhyszAC, Let these be designated 
Dozq and Dy respectively. These terms plot as 


shown in figure 6. Also plotted is D,;7, the total 
dissipation factor of C). 
Der Dep T Dyra T Doz». (5.8) 


If we plot also the dissipation factor of the specimen 
holder, Dyr, and Dy, the results are as shown in figure 
6. Once again we see that Dy plots as a straight 
line, with a negative intercept equal to the value of 
D7 at AC=0. The slope of the plot of Dy is now 
different from that of the plot of Dy, because of the 
fact that D.,, is not constant. 

Because the purpose of all these hypotheses is to 
provide a means of evaluating the various compo- 
nents of dissipation factor and to separate the dissi- 
pation factor of the specimen holder from that of the 
bridge standard capacitor, let us suppose that we 
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Figure 5. A plot similar to figure 4, but including now a 
component of dissipation factor due to series resistance 
the leads to the specimen holder, Dy. 


in 
The total] dissipation factor of the holder, Dar 


Dur+Duz,is plotted, and now 
Du=Dur—Dsy 


(5.3), we | 
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FicurE 6. A plot similar to figure 5, including components of 
dissipation factor due to series resistance in the leads to Co, 
the bridge capacitor. 

As discussed in the text, there are two such components, one constant with 


capacitance, the other proportional to AC. Der is the sum of all components of 
the dissipation factor of C>. 


have just made a set of measurements that yielded 
the plot * of Dy seen in figure 6. The slope of this 
curve will be the same as that of the plot of Dy only 
if the slope of Dy, is zero. To determine the slope of 
Dor, we can make another set of measurements 
starting from a different C,’. This would change 
the value of D.,,, and would permit an evaluation of 
this component of dissipation factor. Another re- 
vealing test would be a set of measurements at 
different frequencies. It can be seen that such 
measurements would permit the separation § of 
dissipation factor components due to series resistance, 
in either the bridge standard capacitor leads or the 
holder leads. The behavior of the slope and inter- 
cept of the plot of jy, would indicate the presence or 
absence of appreciable resistance in the leads of the 
bridge standard capacitor. 

We have now dealt with losses in a surface film, 
and losses due to series resistance. Let us turn our 
attention to losses due to a fixed parallel conductance. 

If a fixed parallel conductance occurs in the branch 
containing the bridge standard capacitor, it is readily 
apparent from eq (5.1) that it will not make any con- 
tribution to the apparent dissipation factor of the 
unknown being measured, for if G@,’=G,’’, the cor- 
rection term is zero. 

If, however, there is a fixed conductance in paral- 
lel with the specimen holder, it will contribute a 
component of dissipation factor that varies with 
capacitance in the manner shown by the plot labelled 
Dyp in figure 7. When this is added to the other 
components of dissipation factor of the holder, the 
sum is Dyr, plotted in figure 7. Subtracting from 
this D2, the total dissipation factor of C2, we have 
Dy, also plotted in figure 7. It is evident that this 
plot behaves in a manner quite different from those 
previously considered, and a plot of measured values 








* Combining equations 5.2, 5.3, 5.4, 5.5, 5.6 and 5.7, we can write the collected 
expression for Dm» 
Dy =wRus tACr+eRusrAC,—wRogr’ Cz —2wRy81,Cy! —wRs AC 
in Dy 
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Figure 7. A plot similar to figure 6, including a dissipation 


factor component due to a fixed parallel conductance in the 
holder, Dyp. 


A similar conductance in C,; would produce no effect in the measurements 
see text 


of Dy would readily reveal such a parallel conduct- 
ance. In practice, such a conductance is not to be 
expected in the present set of circumstances, and 
actual measurements confirmed the absence of paral- 
lel conductance. 


5.4. Results of Measurements 


From the foregoing discussion, it is evident that 
we have available a means of evaluating the dissipa- 
tion factor of the empty specimen holder, and also 
that of the bridge standard capacitor, Cy. It is only 
necessary to use the Schering bridge to measure the 
apparent dissipation factor, 1, at several values of 
Cy, the capacitance of the specimen holder, recalling 
that (footnote 3) CymACs. 

Measurements were made of the losses of the new 
specimen holder, known as Holder No. 4, using Gen- 
eral Radio Co. Type 722—-DS829, Serial 7128 as the 
standard capacitor. 

The measurements were made during the winter 
months, with room relative humidity ranging be- 
tween 10 percent and 40 percent, mostly at about 
20 percent. 

The bridge standard capacitor, C,, was maintained 
in an enclosure dried by a stream of dry nitrogen. 
This nitrogen stream was dried by bubbling it 
through concentrated sulfuric acid, thence passing 
it through a drying column containing calcium chlo- 
ride or calcium sulfate, and finally through a dry 
cotton filter. 

The effect of relative humidity on the dissipation 
factor of the holder was of considerable importance, 
because the holder was constructed to measure the 
losses in low-loss materials under different conditions 
of humidity, principally very dry and 52 percent 
relative humidity, in order to compare the results. 
Therefore, measurements were made on the empty 
holder under two conditions of relative humidity. 
The holder was enclosed in a chamber which was 
very thoroughly dried by exposing therein a gener- 
ous quantity of P,O;, which held the relative hu- 
midity well below 1!4 percent, which was the lower 





limit of the relative humidity measuring apparatus 
used. By exposing within this chamber a saturated 
solution of Mg(NQs)., the relative humidity could 
be maintained very close to 52 percent. 

Note that the controlled atmosphere applies to the 
specimen holder only, and not to the bridge standard 
capacitor, which for all measurements was dried by 
a stream of dry nitrogen. 

Two typical sets of measurements are presented in 
figure 8, for the purpose of showing that a straight 
line provides a reasonable fit for all the data. There 
is no pronounced tendency to exhibit a curve, such as 
the plot of Dy in figure 7, thus indicating that there 
is no appreciable parallel conductance in the holder. 
Note that the dissipation factors are all very small, 
and the slope is also small. It is not profitable to try 
to place any interpretation on the small digressions 
from the straight line, because the resolution of the 
apparatus is not great enough to permit more accu- 
rate measurements of such small dissipation factors. 
The resolution is much better at high capacitances 
than at the low end. Equation (5.1) gives, neglecting 
the correction term, 


wR,O3AC, 
) = : 2 7 
I H AC ', | 


from which it is seen that as the measured dissipation 
factor becomes very small, AC; becomes very small. 
Since it is very difficult to determine very small differ- 
ences in relatively large quantities, the method gives 
rather poor resolution at the lower values of dissipa- 
tion factor, and this effect is heightened by the fact 
that the smallest dissipation factors occur when AC; 
is smallest. For these reasons, no measurements at 
AC, less than 20 pf are included. 

Most of the dielectric specimens for which this new 
holder was designed will have capacitances of 5 to 50 
pf. Referring again to figure 8, notice that in this 
range the dissipation factor does not vary more than 
about 4 x 107°, or in terms of the loss angle,*® about 4 
microradians (yradians). 
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Figure 8. Typical results of actual measurement of dissipation 
factor of the holder under two widely different conditions, 
showing that straight-line plots adequately represent the data, 
and that the dissipation foctors are small. 


§ The dissipation factor, D, is the tangent of the loss angle. D=tan 6. For 
small angles 6&tan 4, and often the term loss angle, expressed in w radians is used. 


107 








The magnitudes of the intercepts of the two sets of 
measurements are less than 2 wradians. In the light 
of the preceding analysis, this can be interpreted to 
show that the sum of the components of loss angle 
caused by the surface film on the electrodes of the 
bridge standard capacitor, C,, and that caused by 
series resistance in the leads to C, is less than 2 
pradians at AC,=0. This sets a limit to the slope of 
the plot of the total dissipation factor of C, versus 


AC,, as is evident from examination of eq (5.3). For 
convenience this expression is rewritten here: 
Dy =Dy+ 2whos,( 2 + whos, AC, + wl hosp’. (5.3) 


Note that the third term on the right is the same as 
the second term on the right when AC,=2C,’. Thus 
if the loss angle of the bridge standard capacitor is 
less than 2 uwradians at AC,=0, it must be less than 4 
pradians at AC,—2C,’=40.000 pf. Note also that 
this is true at the highest frequency used, 100 ke/s. 
Even if the third term on the right were at its maxi- 
mum, almost 2 wradians, at this frequency, it would 
be only one-tenth of this or 0.2 ywradian at 10 ke/s, 
and would be no more than 0.02 yradian at 1 ke/s. 

To sum up this line of reasoning, we have proven 
that at frequencies up to 100 ke/s, and at specimen 
holder capacitances up to 40 pf (CyAC,), the maxi- 
mum error in the determination of the loss angle of 
the empty specimen holder due to neglected series 
resistance in the leads to the bridge standard capaci- 
tor is no greater than 2 wradians. Other tests show 
that the actual error due to this cause is much smaller 
than this maximum. Such tests include measure- 
ments made at various frequencies, from other initial 
values of C,, and on actual specimens of very low 
loss materials, such as polyethylene. 

Our early measurements showed a much greater 
slope for Dy at 100 ke/s. Examination showed this 
loss to be attributable to the series resistance of cer- 
tain leads used at that time. The offending cable was 
removed and replaced with a type with much lower 
resistance, and a great improvement was observed. 
The effect of the residual series resistance in the leads 
to the specimen holder are still observable in figure 9, 
where the results of measurements at frequencies of 
100 c/s, 1 ke/s, 10 ke/s, and 100 ke/s are shown to- 
gether. Note that the slope of the 100 ke/s plot is 
greater than the other slopes. The results shown in 
figure 9 show that at all measured frequencies except 
100 c/s, the intercept is about 1 wradian or less. This 
confirms and emphasizes that the series resistance in 
in the leads to C, is low. 

The data presented in figure 9 were obtained 
with the specimen holder maintained in a very dry 
atmosphere. In figure 10 we see the results obtained 
when this atmosphere had a relative humidity of 
about 52 percent. Note that in both these cases, 
the bridge standard capacitor, C,, was maintained 
as dry as could be attained with a stream of dry 
nitrogen flowing through the enclosure containing 
it. This means that losses in (,, indicated prin- 
cipally by the intercepts, should remain nearly 
unchanged for the two sets of data. The higher 
relative humidity might be expected to increase 











6xI0% J 

NOLDER rh < 1.5% 

100 kcA 

4r 4 

2F J 

D * 4 

os 10 kc 

1@) 





nD 


100 c/s | 
a ae 4 


-4- 4 
4 


12) 20 








1 4. 


80 too 





A. A. 

40 60 
CAPACITANCE (ACo), pf 
Figure 9. Results obtained at 100 c/s, 1 ke/s, 10 ke/s, and 


100kc/s, with the holder maintained in a dry (relative humidity 
less than 1.5 percent) atmosphere. 
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FicurE 10. Results obtained at 100 c/s, 1 ke/s, 10kce/s, and 100 
ke/s, with the holder maintained in an atmosphere of 52 


percent relative humidity. 


the loss in the surface film on the electrodes of the 
holder. This would be evidenced by an increased 
slope for Dy plots at the higher relative humidity. 

Comparison of figures 9 and 10 shows that actual 
measurements agree with the behavior anticipated. 
The intercepts, within experimental error, are the 
same for both cases, while the slopes are noticeably 
greater for the higher humidity case. 

The agreement of the intercepts is even more 
evident in figure 11, where values of D, obtained 
under the two different conditions of relative 
humidity surrounding the holder are plotted to- 
gether. The agreement is satisfactory, and the 
curve drawn through the average values shows the 
variation of D, as a function of frequency. The 
higher dissipation factor at the lower frequencies is 
evident here, as well as a very slight increase at the 
high end of the frequency scale. This latter is 
probably due to the residual series resistance of the 
leads to the bridge standard capacitor, C). 

Early in the course of this investigation it was 
discovered that the apparent loss of the empty 
specimen holder was profoundly affected by the 
presence of small amounts of foreign materials on 
the electrode surfaces. It was necessary to keep 
the electrodes scrupulously clean, by occasionally 


cleaning them with distilled water followed by 
petroleum ether. It was found that electrodes 
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Figure 11. The dissipation factor, D2, of the bridge standard 
capacitor, C2, plotted as a function of frequency. 
Each value is obtained from the intercept of one of the preceding plots, one set 
from the dry holder plots of figure 9, the other set from the 52 percent relative 


humidity holder plots of figure 10. Note that C) was maintained in an en- 
closure dried by a stream of dry nitrogen for all measurements (see text). 


applied to dielectric specimens must be very clean, 
to avoid carrying dust particles, etc. to the holder 
electrodes. The care previously used in the handling 
of specimens was found to be inadequate, and new 
precautions were initiated, especially in the case 
of low loss materials, where a loss contribution due 
to foreign matter would be a greater percentage 
of the whole loss. 


5.5. Correction for Holder Losses 
Measurements 


in Dielectric 


After the values of the dissipation factors of the 
bridge standard capacitor and the empty specimen 
holder have been determined, the next problem is 
how to use them properly to correct the measured 
dissipation factor of a dielectric specimen. 

In all cases, it is clear that D,, the dissipation 
factor of the bridge standard capacitor, should be 
added to the measured value of dissipation factor 
of the specimen, Dy. 

If the dielectric specimen has no contact electrodes 
applied to it, and is larger than the guarded elec- 
trode, it is evident that the measured dissipation 
factor will consist of two components, one due to 
losses in the dielectric material, the other due to 
losses in the holder. The latter component, which 
can be determined by the method described herein, 
can be subtracted from the measured value to yield 
the dissipation factor of the dielectric material 
itself, taking into account the fact that the dielectric 
constant of the specimen multiplies the effect of 


Dy. The correct relation to use in this case is: 
D,=Dy—e&DatDz, (5.9) 
where 
D,=dissipation factor of the specimen 
and 


e, dielectric constant of the specimen. 

In the case of a specimen with contact electrodes 
applied, this question must be answered: does the 
surface film remain intact, or is electrical contact 
made directly from the conducting contact electrode 
to the conducting substrate of the holder electrode 








material? If the latter is true, then the dielectric 
film is in effect short-circuited, and does not affect 
the measured dissipation factor. The problem is 
then shifted to that of determining what losses occur 
in the contact electrode material, in the interface of 
the contact electrode and the dielectric specimen, 
and in the material serving as adhesive, if any. 
This is a separate matter, and will not be further 
treated here. 

If the specimen is smaller than the guarded 
electrode, so that some air capacitance is included 
in the measurement, the problem can be resolved 
into two parts, one dealing with the portion involving 
air capacitance in the space surrounding the speci- 
men, the other dealing with the portion of the be a 
between the electrodes that is occupied by the dielec- 
tric specimen, taking into account its contact elec- 
trodes, if any. 

Admittedly, some uncertainties arise here with 
respect to what happens in the case where contact 
electrodes are used. However, the fact that the 
surface film contributes (in this case) such a small 
component of loss, saves one from commiting a 
grave error, and it is probably best to make reason- 
able assumptions and proceed, secure in the knowl- 
edge that at worst the error will not exceed a few 
microradians. 


5.6. Implications for Further Investigation 


A second specimen holder is under construction, 
quite similar in basic design to the one described 
herein, except that provision is being made to make 
the electrodes interchangeable, and several sets of 
electrodes are in preparation. These electrodes are 
made of various materials, and will be used to deter- 
mine which material yields lowest losses in surface 
layers. In addition to measuring several solid elee- 
trodes, it is intended that platings of various kinds 
will be applied to some of the electrode sets, and the 
losses in the surfaces measured. 


6. Summary and Conclusions 


A guard-ring micrometer-electrode specimen holder, 
especially designed for use for precision measure- 
ments of the dielectric properties of dielectric refer- 
ence standard materials, has been constructed and 
evaluated. A technique has been devised for deter- 
mining the spacing between the electrodes to +-1 
u, using the micrometer. 

It was found that resistance in cables often used 
to connect the components of a bridge can produce 
losses which ave appreciable, and unless a correction 
is made for such losses, they can produce errors in 
the values of the dissipation factor of the holder, 
and also that of the bridge standard capacitor. A 
preferable method is to keep these resistances as low 
as possible by proper choice of cable type. 

It was found also that minute quantities of foreign 
matter (such as dust or grease) on the holder elec- 
trodes will contribute an appreciable increment to 
the dissipation factor. It is necessary to be very 
careful to keep the electrodes scrupulously clean 
when low-loss materials are being measured. 
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The dissipation factor of the carefully cleaned 
holder has been measured, and it has been concluded 
that at frequencies in the range 100 c/s to 100 ke/s, 
almost all of the loss is in the surface film on the 
Type 304 stainless steel electrodes. The value of 
this dissipation factor was found to be directly 
proportional to the capacitance of the holder, and 


therefore inversely proportional to the separ ation of 


the electrodes. The relative humidity of the at- 
mosphere surrounding the holder had an effect on 
the surface losses. The dissipation factor varied 
from about 1>10~* per picofarad at 1 ke/s in a very 
dry atmosphere, to about 710~-® per picofarad at 
100 ke/s with an atmosphere of about 52 percent 
relative humidity. 

Knowledge of the values of these empty specimen 
holder losses is very useful for making corrections 
to the measured loss of dielectric materials with low 
losses, where the loss contribution of the electrode 
surface may be an appreciable fraction of the total 
loss. 

Measurements made by the substitution method 
with a Schering bridge should also be corrected for 
the error due to the dissipation factor of the bridge 
standard capacitor. The method used vielded the 
value of this correction. The loss in a General 
Radio Type 722—DS29 variable air capacitor was 
concluded to be almost entirely in the electrode 
surfaces, and it was confirmed that the dissipation 
factor at any single frequency is almost constant 
over the entire capacitance range. When the 
capacitor was maintained in an enclosure dried by a 
stream of dry nitrogen, the dissipation factor was 
found to range from less than 0.5 10-° at 10 ke/s, 
to about 3.5107 at 100 e/s. 

With this new specimen holder, it possible to 
measure the spacing of the electrodes with about 
the same accuracy as that with which the thickness 
of the specimens can be determined. With the 
application of corre ‘ctions determined by the methods 
described in this paper, the holder can be used to 
measure the loss angle of low-loss materials to 
within a few microradians. 


7. Appendix 


The exact expression for Cy as a function of the 
Schering bridge parameters 1s 


‘Ad ld yr yr 
v1 G3 ( 37 a i 3 


Cy= C2 See 04+ G' Se as (7.1) 
o * G.*-+a'C;* dil. G, ?7+w°C;? ee 
The expression for Dy is 
wC,’ GY 43400303 
y) —s ? Cs ( — > aC -)-C | 
Di " "Cr " \ G3 *+w*Cs 
in G3 a C: (“3 oa) = G; ‘ (7 2) 
“Oh H I G, . G ir oy wf ar wl H 








where 
Dyz=dissipation factor of the specimen holder 
and, in general, any specimen contained 
therein ; 
C3=initial value of C,, with switch 1, figure 3, 
in position A; 
C:’=final value of C,, with switch 1 in position 


B; 
AC,=C,'’—C,’: 


(’3=initial value of C3, with switch 1 in position 
a: 

(’;’=final value of Cy, with switch 1 in position 
B; 

AC’, Obs $A 

Gi=initial value of equivalent parallel con- 
ductance associated with ¢ 

(,'=final value of equivalent parallel conduc- 
tance associated with Cs; 

(7;=initial value of equivalent parallel conduc- 


tance associated with Cy; and 
final value of equivalent parallel condue- 
tance associated with C4. 
Since G; is principally the conductance of a fixed 
precision resistor of more than 10° ohms, it 
safe to assume that 


Is 


not 


We can now write 


: 
wl’, 


ve a G4 ang A 
2 IC ( Bi Sad ad )—C 


G:C 
, Gy [ ey = 
+—? ap | (7.3 
wl'y ° wy 


For the Schering bridge used, w/G,; was held constant 
at 24107, C2’ never exceeded 1100 pf, and CY —CI 
never exceeded 1000 pf. Under conditions 
we nay write 


Dy 


ONC! 
G? 


aT, 


these 


wl,’ [_,,,(G3+w Gy’ oie 
Du=77-| C3'( +- E2, (7.4) 
G3Uy Gr wf wl Ho Wy 
with no appreciable error in Dy. And 
cc 7, —G, 
D = (CU; — (7.5) 
H G ( a \ + w Cw ‘ y 
will give no more than 0.5 percent error in Dy. For 
less extreme values of C2’, the error is less. The 
present work involved very low Dy, and therefore 


very small C3’—C3 resulting in wholly 


negligible 
errors due to terms omitted from eq (7.5). 


The first term on the right in eq (7.5) may be 
expressed as 
whCy , , 
— (C3'—C)» 
Cx 
yielding 
R ( oA AC A “Ed ol = 
Dy Y + “2 3, (7.0) 
Cn wx 
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In eqs (7.2) through (7.6), the exact expression 
for Cy given by eq (7.1) may be used. For the 


present work, however, C3; was about 100 pf and 
C3’ was less than 150 pf, yielding an error in Cy less 
than 0.003 percent if the correction terms in eq (7.1) 
are neglected. Therefore, we may with negligible 
further error write AC, for Cy in the preceding equa- 
tions. Thus our expression for Dy becomes 


wh ,( pad AC / Gs’ GQ: 


) ; 
Du AC; wAC, 


(7.7) 
The last term in eq (7.7) is a correction term 
necessary if Gi’#G). If this correction term is ig- 
nored, or is negligible, Dy is the term wl,C,’’AC,;/AC, 
which we will call )),,, the measured or uncorrected 


value of Dy. Then 
G',' —G; 
Dy=Dy a; (¢.3S) 
wA 2 
Kor present purposes it will be useful to express 


this in terms of equivalent series resistances. 
In parallel notation D=Gp/wl'p, 
where 
Gp—equivalent parallel conductance 
(’p—equivalent parallel capacitance. 


Therefore, let us write 


G’—G, = CG C345 CS es, 
2 a =F ee. Di —-—* D). (7.9) 
wAfl 2 AC owl 2 AC ow 2 AC 2 2 AC 2 * 
In series notation, D=wh.( 
where 
?s=equivalent series resistance 
Cs=equivalent series capacitance. 


Therefore if we let 

Chs=series equivalent of C; 

Czg=series equivalent of C%’ 

Rjs;=equivalent series resistance corresponding to @ 
R35=equivalent series resistance corresponding to @}’ 
then 


"6 qn ¢ eee 44 ea) ee 
= EA _— = Do= = @ togl s~ = 
AC.” * st.” * pidoadaee oA 


Because the dissipation factors involved in this work 
are very small, the series equivalent 
very nearly 
That is 


capacitance is 


the same as the parallel capacitance. 


Crs@Co, and Cys@C%’. 
Therefore 
$44 ore * ore Ce eRe Rss 
—— wh? ( = F wh sof 2S y si me 11) 
AC AC; AC AC, 








Let us now consider the series resistance Pos to 
be composed of two components: (1) a series resist- 
ance term /?,s; representing series ohmic resistance 
in leads, connections ete. and (2) a term [sr repre- 
senting the series equivalent resistance of a surface 
film on the electrodes. That is 


Res Rost T Rise 
Rys= Rysi+ Ryser 


(7.12) 


Because Ps; represents series ohmic resistance, this 


term is inde ~pendent of C;. That is 
Rosr Resz Rosr 
and 
) "a(R 2SL T Ras .) Cs*w(Rest4 | ) 
Du me : AG a — = 
\é 13) 
(03° —CP)wRos 
) Pee —2SL 
Dut 
le. Cyw tos 
nd Bn F. (7.14 
+ AC, AC, 7.14) 
Because C3/=C}+AC, 
(¢ ga -Cy \oRea: wl 2 Ry spC 
)u=D -_ 
Dy ut AC, r AC) 
CY Riehl, wR mei 
—=- (7.15) 
AC { AC . ‘ / 
The conductance of the surface film, 1/Rossr, is 


proportional to the active area of the capacitor 
plates, which may be expressed 


Ross aA. 


From this Rosp=1/aA. 
The capacitance C, is also proportional to the active 
area of the capacitor plates, or 

G =BA. 


The product P.-C, is now seen to be a constant 
because 


pA 8B 


aA Q@ 


R, spl 


? 


the ratio of two proportionality constants. 


From 
this we see that 


Riel! = Roop 
Therefore 
(¢ yer8_¢ 2 )oRasr ws Mis 2sPL 
AC; AC; 
_ Wl oR sar i. 
AC, 


Dy=Du + 


24 "R's sp 


on (1%) Roe " ‘ ates 
Da ; AC PSE + oof 2 tosr* 7.17) 








Attacking now the second term on the right, we can 
write 
( phere 


C)(Cy’ +¢ "oR osr 


Du=Du+* ai 


- +O DR ose 
(7.18) 
Since C,’—C,;=AC, 


Dy=Dy+ (Cy + Cp) oRos14 (7.19) 


y »” 
wO Rasp, 


and because C;’=AC,+ C3, 


Da = Da T 2wO>Rosr + wAC oftosr + wl oltosr. (7.20) 
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A Bolometer Mount Efficiency Measurement Technique’ 
G. F. Engen 


(December 14, 1960) 


In the measurement of microwave power by means of the bolometric technique, the 
efficiency of the bolometer mount must be measured and applied as a correction in order 
to meet many of the accuracy requirements of today’s technology. 

An impedance technique of determining this efficiency was proposed over ten years 
ago, but has found but little use to date because of the rather severe performance require- 
ments imposed on the attendant instrumentation. This paper describes an improved 
method of implementing the technique which is based on the reflectometer concept. 

A particularly attractive feature of the new method is its substantial independence of 
connector discontinuity which has been a particularly troublesome source of error in coaxial 
systems. In further contrast with some of the earlier proposals for implementing this 
impedance technique, the new method is readily applicable to both matched and unmatched 
mounts and does not require mathematical approximations (in the first order theory). 

\ comprehensive error analysis indicates that an accuracy of 0.5 percent is possible in 
the existing state of the art. 


1. Background troublesome source of error in coaxial systems. This 
technique constitutes an improved method of making 
The “impedance”? method of measuring bolometer the measurements implicit in the “impedance”’ 
mount efficiency devised by Kerns [1]! is one of the method formulated by Kerns, and the subsequent 
few basic techniques developed thus far for deter- discussion and error agalysis will be limited to the 
mining this parameter. As originally outlined, how- procedures to be described. For a discussion of the 
ever, the accuracy which could be achieved was | More basic postulates upon which the procedure is 
rather severely limited by the state of the impedance based the reader is referred to the original paper [1]. 
measuring art, and this led in turn to the develop- | In particular it should be noted that the technique is, 
ment of a number of modifications of the technique | © Its present form, applicable to barretter but not 
with the objective of reducing the overall error. | thermistor type bolometers.? 
Beatty [2], for example, proposed a modification 


based on certain mathematical approximations and 2. Introduction 

restrictions in generality which provided improved 

accuracy, but the associated operating procedures It was shown by Kerns [1] that thé efficiency of a 
proved to be nonetheless time consuming and exacting, bolometer mount. may, under sulable conditions. 
and the accessory instrumentation was never devel- | be determined from three impedyice measurements 


oped or refined to the point where one was, on a | at the bolometer mount input tefminals correspond- 
routine basis, able to place a great deal of confidence ing to three different values of Dolometer resistance. 
in the results. Additional refinements or modifica- | Tf one of these resistance valu is chosen to coincide 
tions of the technique have also been suggested by | with the value of resistancefor which the efficiency 
Weinschel [3], Ginzton [4], Lane [5], and merase | is deaived. this reanti weoe expressed in terms of 
others, but it is probably safe to say that none of | the input reflection coefftients thus [2]: 
these proposals has, as vet, come into widespread use. 
At the Boulder Laboratories of the National . 
Bureau of Standards another variation of the im- n=K 
pedance method has been deve loped which provides 
improved accuracy and simplified operational pro- 
cedures but, unlike the earlier modifications, requires 
neither mathematical approximations (in the first 
order theory), nor restrictions in the generality of the 






where the T,, To, W3 are the reflection coefficients 
corresponding to bdlometer resistances R,, Rs, Rs, 
respectively, 7 is thé efficiency when the bolometer 


method. <A particularly attractive feature of this | — : 
new version is its substantial independence of con- 2 In order to apply the technique to a given bolometer mount, it is a sufficient 
‘ ; A : ; (but not necessary) condition that it be possible to determine the microwave 
nector diseont muity, W hich has been an especially impedance of the bolometer element from measurements of its d-c resistance. 
A somewhat different sufficient condition is that the resistive component of the 
se microwave impedance be related to the d-c value by means of a real proportion- 
ility factor (whose magnitude need not be known), and that the reactive compo- 
*A preliminary report of this technique was given at a joint meeting of the nent be constant or independent of the d-c resistance. Necessary conditions have 
International Scientific Radio Union and the Institute of Radio Engineers, April not been determined, but it has not been possible, thus far, to relate the microwave 
26, 1958, Washington, D.C impedance to the d-c resistance in an adequate fashion in order to utilize this 

Figures in brackets inc Li ite the literature references at the end of this paper. technique with thermistors, 
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resistance has the value /?, (typically 200 ohms), and 


2R(R;—R, ) 


Kk=|— - For the special case of a 
\(R.—R))(R3—Rz2)| ; ’ 
“matched”? mount (T,=0) the expression becomes 
-| Ts 
7=K\= — |e (2) 
r;—T, 


The determination of the factor AK, while present- 
ing a number of practical problems, does not pose 
the difficulties encountered in measuring the factor 
containing the I’s. It will be noted that this latter 
factor is characterized by differences or changes in 


the reflection coefficient values as the bolometer 
resistance is varied. Because the resistance varia- 
tion which can be achieved with the available 


barretter elements is at best rather limited, the 
corresponding changes in the reflection coefficient 
are not large, and if the individual [’s are measured 
(by means of a slotted line for example) a rather 
severe requirement on the accessory instrumentation 
results. Beatty [2] has shown that an error of +1 
percent in the VSWR measurements may lead to an 


error as large as +6 percent in the determined 
efficiency. 
More recently, the reflectometer technique has 


found increased usage in the measurement of reflec- 
tion coefficient m: agnitudes. This device, shown in 
figure 1, ideally viel Is a response of the form: 


=f iT (3) 


where 6; and b, are the wave amplitudes of the signals 
at the respective detectors, T is the reflection coeffi- 
cient of the bolometer mount or other load terminat- 
ing arm 2, and & is a (real) constant whose value 
may be determined by observing the detector outputs 
with arm 2 terminated by a load of known reflection 

a fixed short for example. 

While the usual reflectometer fails to produce 
this type of response because of imperfect directivity 
and other deviations from ideal behavior of the 
directional couplers employed, it is possible to com- 
pensate for these imperfections at a single 
by means of auxiliary tuners [6], 
9 


frequency 
as shown in figure 
Procedures for the adjustment of these tuners, 
such that the response of eq (3) may be realized, 
were ‘given in the cited reference [6]. In essence 
the procedure is to adjust the tuning transformer 
T, such that the directivity of the associated direc- 


tional coupler becomes infinite, while 7, is adjusted 
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Ficure 1. Ideal reflectometer 
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FiGuRE 2. Reflectometer with auxiliary tuners. 

such that the reflection coefficient of the equivalent 
generator, at the terminals where the 
mount is connected, vanishes. 

The reflectometer does not, in this form, provide 
a determination of the argument or phase angle of 
the reflection coefficient, however Beatty [2] has 
noted that for matched mounts of high efficiency 
the reflection coefficient vectors T, and Ty; are nearly 
colinear so that to a good approximation [';—T, 

Tr, + T,. (The proof of a more general counter- 
part of this statement will be given later in this 
paper.) To the extent of the validity of this approx- 
imation, the cited reflectometer technique may thus 
be emploved to determine the efficiency of matched 
bolometer mounts by means of eq (2). The appli- 
cation of the reflectometer method to the more 
general problem of unmatched mounts is, however, 
by no means obvious. 

Within the limitations discussed, 
method is practical in rectangular 
tems. In coaxial systems, however, 
leaves a great deal to be desired 
gular waveguide, as will be subsequently demon- 
strated, the procedure, as given, is unnecessarily 
complex for the application. 

The problem in coaxial measurements centers 
around the impedance discontinuity which is part of 
the coaxial connector. In order to provide a useful 
result, an overall accuracy in the efficiency deter- 
mination of 1 percent would appear to be a reasonable 
goal, and since this figure will include contributions 
from a number of different sources it is desirable to 
reduce the individual contributions to the order of 
0.1 percent. Applying this criterion to the amount 
of mismatch permissible in the equivalent generator 
impedance implies a generator match with a VSWR 
of 1.002 or less. While this value may be realized in 
rectangular waveguide systems by means of suffi- 
ciently refined tec hniques, the prospect of achieving 
such a result in a coaxial system, where > ep 
discontinuities with VSWR’s of the order of are 
typical, is certainly far from encouraging. 

It may be noted, in passing, that this problem has 
appeared so regularly in attempts to refine this 
technique for use with coaxial type bolometer mounts 
that one is led to wonder if any significant improve- 
ment in the technique is possible with the existing 
coaxial connectors. For example in Beatty’s modi- 
fication, which was based on a variation of the slotted 
line technique, the measured efficiency included the 


bolometer 


the described 
waveguide sys- 
the procedure 
, and even in rectan- 


114 








losses in that portion of the slotted section between 
the probe position and bolometer mount, leaving one 
with the awkward additional problem of measuring 
and applying as a correction the losses in this portion 
of the slotted line. 


General Theory 


The contributions of the new method to be de- 
scribed in the following paragraphs include an 
“exact”? procedure for determining the vector differ- 
ence between two reflection coefficient values, and a 
tuning procedure which virtually eliminates the 
potential error from connector discontinuity. 

The magnitude of the vector difference between | 
two reflection coefficient values may be determined | 
with the help of the cireuit arrangement shown in | 
figure 3. For convenience an ideal reflectometer | 
will be assumed. Terms of the form [;—T,| may 
be determined as follows: The barretter resistance 


is first adjusted such that the bolometer mount has | 


the input reflection coefficient T). 
detector 3 is then ‘‘nulled”’ 
of the auxiliary signal 


The signal at the 
or balanced out by means 


channel. Under these cir- 


cumstances the signal provided by the auxiliary 
channel is proportional to —T, the barretter 
resistance Is now adjusted to produce I's, and if 
suitable precautions (to be described later) have 


been taken to isolate the two channels, the detector 3 
signal will now be proportional to ['3—T) as required. 
In order 
dure it will prove helpful to consider a special case. 
If a perfectly matched bolometer mount is assumed, 
the tuning of the reflectometer may be effected as 
follows. First, with the bolometer mount connected | 
and the bolometer resistance adjusted to the value 
corresponding to the perfeet match, the transformer 
T, is adjusted for a null in arm 3, corresponding to 
infinite directivity for the associated coupler. 
the bolometer mount is replaced by a sliding short 
and the transformer 7’, adjusted such that the ratio 
2 is constant for all positions of the short. As was 
) 
explained in the reference [6], this adjustment pro- 
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duces an equivalent generator match (provided that 
the previous adjustment has been correctly made), 
and having completed these two adjustments, the 
measurement of the mount efficiency may be carried 
out using eq (2) as already described. (In practice 
additional tuning elements, whose use will be de- 
scribed later, are required to isolate the auxiliary 
channel.) 

It should be noted that to the extent that the 
(coaxial) connectors may be considered dissipation 
free, the attendant impedance discontinuity has no 
effect upon the sliding short adjustment since an 
ideal sliding short preceeded by a lossless discon- 
tinuity still presents a reflection coefficient of unit 
magnitude and (except for trivial cases) of variable 
phase angle. The place where the connector dis- 
continuity is potentially important is in the initially 
assumed impedance match for the bolometer mount. 
But, as will be proved, the assumed impedance match 
is entirely unnecessary. To be sure, if the mount 
is matched, the measurement may be carried out as 
described; while if the mount is not matched, the 
identical of measurements and operations 
described above will yield the factor 


series 


(T;—T.)(T:—T 2) 

Psa U lr, |?) 
as required in the more general expression for 
mount efficiency. 


In other words, the reflectometer adjustment and 
measurement proceeds on the basis of an assumed 
perfect impedance match for the bolometer mount 
(but employing only the specific procedures men- 
tioned) It then develops that the method gives 
the correct value of efficiency regardless of whether 
the mount is matched or not! The proof for these 
assertions will be found in the following paragraphs. 

It is a general property of a four-arm junction 


(of which the reflectometer of figure 2 may be 
regarded as a special case) that the ratio of the 


emergent wave amplitudes in arms 3 


and 4 may be 
written in the form: 


“<p (4) 
where the A, B, C, D are functions of the parameters 
of the four-arm junction and the detectors termin- 
ating arms 3 and 4, and T is the reflection coefficient 
of the load terminating arm 2. 

If the reflection coefficient of the botometer mount 
is designated by T, when the bolometer at its 
nominal operating value F, (the value at which the 
mount efficiency is to be determined) and the 
junction has been adjusted (by means of 7,) such 
that the signal 6; vanishes when arm 2 is terminated 


is 


by this load (T,), then B=—AT,. A solution of 
eq (4), subject to the condition that the ratio _ 
M4 
remain constant (as provided by adjustment of 7,) 
; : : : BfG 
while the phase of IT varies, yields ( ) r|? 
Lii¢ 1¢ | lase O ( 5 re dD 


where the asterisk (*) denotes the complex conjugate. 








(A second solution, A/C=B/D is trivial since it 
gives for 7 a value which is independent of TI.) 
For a result 
a=(5): Substituting the first and last of these 


results in eq (4) gives: 


4 
sliding short the above becomes 


bs 
b, 


A(t—T;) (5) 
~DQ—Tr*) is 

Substituting I, and YT; for T and taking the absolute 
value of the ratio of the product to the difference 
gives: 
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DI W—rsPs] |I—PiPs|_ [Al |Ps—Tal|P— x 
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The term on the right is just the one required in 
the general expression for mount efficiency, while the 


factor may be determined by observing the sys- 


i 


tem response with the bolometer mount replaced by 


a short or other termination for which |[}=1. This 
may be demonstrated by letting [=e. Then: 

|bs|_|A) | e*—T2|_|Al yy | 1-Tre"* |_|A] 
b| |D| |ime"rs| (DI'* “lide "| (D| 8 








as required. 
Comparison of eqs (5) and (6) with (3) and (2) re- 
spectively, indicates that the measurement may be 
| 
: ee 
effected as outlined with D| 
— 


TT? 


taking the part of k, and 


- being the more general counterpart of T. 


4. Circuit Arrangement and Tuning 
Procedure 


A practical waveguide system utilizing these con- 
cepts is shown in figure 4. For convenience the 
directional couplers will be designated by P, Q, R, 
and S. The auxiliary arm between couplers F and 
S contains an on-off switch, phase shifter, and atten- 
uator as required to balance out the signals from 
coupler P. As a matter of operating convenience 
it is generally helpful (but not essential) to adjust 
the system in such a way that the signal at detector 
4 is nominally constant or independent of T. This 
is achieved by an adjustment of 7’, to be described. 

In order to obtain the vector difference between 
two reflection coefficients by the procedure described 
earlier, it is important to investigate the possible 
effect of variations in T upon the signal delivered to 
the detector 3 via the auxiliary arm. It is intui- 
tively evident (and proof will be omitted) that the 
desired operation will be achieved if the following 
criteria are satisfied. 








DETECTOR (3) 


TO 
SIGNAL 
SOURCE 











FriGure 4, 


Waveguide system for measuring bolometer mount 
efficiency. 


The signal coupled into the sidearm of coupler 
FR is independent of I’, or more specifically since it is 
the ratio of the signal in detector 3 to detector 4 
which is observed, this ratio must be independent 
of T. 

The impedance “looking into” 
sidearm is independent of I. 

As a consequence of criterion 2, it will also be true 
that the impedance looking into this coupler and 
detector from the main arm will be independent of 
the operation of the phase shifter, attenuator, or 
switch in the auxiliary arm. 

The adjustment of the tuning transformers to 
achieve these conditions may be done by the follow- 
ing steps: 

(1) Directional coupler P and tuner T, are tem- 
porarily removed from the system. In their place a 
sliding short is connected to coupler Q, and an 
auxiliary detector connected at the point X. The 
microwave signal is then delivered to the detector 3 
via the auxiliary arm, and tuner 7, is adjusted for : 
null in the auxiliary detector connected at point X. 
This satisfies criterion 2 above. 

Tuner 7’, is now adjusted such that the ratio 


coupler S at its 


- 


\bs| 
|b, 
and TJ, is similarly adjusted such that |b,| remains 
constant as the position of the short is varied. If, 
as is often the case, the type of instrumentation 
2) is of such a nature 
bs 
b,| 
b,| 

presence of large variations in |b,|, 7, may be first 
adjusted such that |);) remains constant as the 
sliding short is varied, and 7, then adjusted for a 
constant value of |4,|. Ideally, the two procedures 
will yield identical results, while in practice, depend- 
ing upon the degree of precision toward which one is 
working, it may prove desirable to check the result 
for conformity to the criteria stated earlier. 

If attention is centered on the four-arm junction 
comprised of couplers R, Q, and tuner T,, it is of 
interest to note that the adjustment of 7, to make 


is constant for all positions of the sliding short, 


employed to measure |),| and 


of 


| 


that the measurement is awkward in the 


b.| 

\b,! 
solution A/C= B/D rejected in the earlier discussion 
of areflectometer. Although a sliding short has been 
specified, other types of variable loads, for example, 
a fixed short preceded by a variable attenuator, 
might also be employed in this particular adjust- 
ment. For practical purposes, the method of con- 
necting the directional couplers will determine which 
of the solutions A/C=h/D or B/A=(C/D)* 
achieved by the sliding short adjustment, but should 
the question ever arise as to which of the two condi- 
tions has been realized by the given procedure it is 
sufficient to note that the former solution makes 


the ratio a constant corresponds to the “trivial” 


is 


1h. | 

completely independent of T while this is not true 
i“4 
of the latter one. 

(3) The coupler P and tuner 7, are next replaced, 
the bolometer mount connected and sof adjusted for 
a null at detector 3 with the bolometer resistance at 
its nominal value. The switch in the auxiliary arm 
is in the “‘off’’ position for this and the subsequent 
adjustments. 

(4) The bolometer mount is then replaced by the 
sliding short and 7’, adjusted for a constant ratio 


as the position of the short is varied. This com- 


D3 
1b, 
pletes the alinement procedure. 
. . Ig ° 4 9 
If the ratios ; | ure designated by subscripts 1, 3, 
M4 
and S, corresponding to values for the bolometer 
resistance of Ff, 3, and the response when the 
bolometer mount is replaced by a short, the desired 
result may be expressed in terms of the measured 
values: 


b b,| 

bali [Ogi (I's P2)(Py =) (g) 
(32) (%) | . | ITs—Ti}(1 I’, |?) 
| b/s by/;| \bals 


where the first factor in the denominator of the left 
side is determined in the manner described at the 
beginning of section 3. 


5. Practical Considerations 


By of the techniques described in the 
preceding sections, the microwave impedance meas- 
urements required for a determination of mount 
efficiency have been reduced to a set of tuning ad- 
justments and the measurement or determination 

! | ! 
of the ratios ic ’ |Ps| ’ 
\b, 1 lbs 3 


techniques are available. 


means 


ete., for which a variety of 
Depending upon the type 


of instrumentation employed, the ratio i may be 
1%4 
indicated directly, or the terms |b; and {b4! may be 
obtained individually. If the adjustment of 7, has 
been made with sufficient care, the signal {d,! will 
depend only upon the output level of the signal 
source; and if the latter is sufficiently stable, |b,| 








will be constant and thus may be cancelled out of eq 
(7), leaving only the values of |b3| to be determined. 
Indeed, as discussed in the reference [6], if a suf- 
ficiently stable signal source is used the directional 
coupler Q, detector 4 and tuner 7’, may be omitted 
from the system. Under these conditions, step (2) 
of the tuning procedure described above consists 
of adjusting 7, such that |b3| is constant, while the 


| 

other steps follow as given with | 

iM4 

The detection techniques which may be employed 

include audio modulation and detection, heterodyne 

detection, and bolometric power detection. Although 

it is not within the scope of this paper to examine 

each of these methods in detail, it may prove useful 

to call attention to some of their more prominent 
features. 


replaced by |bs|. 


Perhaps the most convenient, but one of the least 
accurate, in the existing state of the art, of the cited 
methods is that of conventional audio modulation 
and detection. The detection equipment may take 
the form of a ratio-meter with barretter detectors 


such that 

tors are specified in order to ensure that the detector 
impedance will be at least nominally constant as re- 
quired above. Alternatively, a standing-wave type 
amplifier may be used, and the signal |b,! applied to 
the automatic gain control channel (if available). 
It should be noted that such an instrument is ac- 
tually a ratio-meter although the permissible excur- 
sion of the input signal applied to the AGC channel 
is usually quite limited. In practice these excursions 
may be held to a small value by adjustment of 7,, 
and once again, if the signal source is sufficiently 
stable and the adjustment of 7, is made with suffi- 
cient care, it is only necessary to measure {b3|, and 
the AGC channel is not required. The chief prob- 
lem with conventional audio techniques appears to 
be in the difficulty of accurately measuring changes 
in the level of an audio signal. In addition, because 
the measurements must be made over a nominal 20 
db dynamic range, the deviations of the barretter 
detector from true square law response would be 
ultimately a source for concern. 

Another objection to the conventional audio 
technique is that the barretter resistance (in the 
mount to be measured) tends to follow the modula- 
tion envelope, in the same manner as occurs in the 
detectors. In order to avoid this difficulty it is 
necessary to hold the microwave power dissipated 
in the barretter to a small fraction, nominally 1 
percent or so, of the total bias power. In addition, 
the resistance excursions which may be realized 
with the commercially available barretters are such 
that the reflected signal is down by another nominal 
20 db, and thus even if the modification (to be 
described below) which permits tight coupling to 
this reflected signal is employed, the signal power 
available is only of the order of a few microwatts, 
which taxes the ability of the best audio detection 
equipment. In spite of these objections, however, 


ee ; 
b,| is indicated directly. Barretter detec- 
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the audio technique does provide a convenient 
method of developing a familiarity with the overall 
system behavior, and as this is bei ing written, some 
of these objections are being overcome [7]. When 
used in conjunction with some of the additional 
refinements to be discussed below, the conventional 
audio technique would probably permit an efficiency 
determination to an accuracy of a few percent 
or better. 

The IF substitution technique of attenuation 
measurement immediately applicable to the 
measurement of relative levels of |b; , while a possible 


for of b, 
figure 5. In figure 5, the 30 Me/s signal derived from 
b, is adjusted by means of the piston attenuator to 
equal in amplitude that obtained from 63, while 
phase balance is obtained from the phase shifter as 
shown. The changes in \b;| with respect to |b.) may 
thus be read directly from the piston attenuator. 
The measurement procedure thus consists of a series 
of nulling adjustments and is independent of fluctua- 
tions in the amplitude of the signal source, or the 
gain of the 30 Me/s null detector. For these 
reasons the scheme appears attractive (on paper 
at least) for the proposed application. 

A preliminary investigation of this technique has 
shown, however, the existence of a number of at- 
tendant practical problems. First, the 30 Me/s cir- 
cuits are frequency sensitive, requiring high stability 
of the intermediate frequency. Second, for proper 
operation the conversion efficiency of the crystal 
mixer should be independent of the changes in level 
of the local oscillator signal which inevitably ac- 
company the operation of the phase shifter, but the 
conversion gain proved to be more sensitive to the 
local oscillator amplitude than anticipated. Third, 
the fact that the null is phase sensitive complicates 
the alignment procedure somewhat, particularly the 
adjustment of 7,. Undoubtedly, a more careful 
examination of the technique will uncover additional 
problems; this technique has not been investigated 
further because the method described below proved 
to be more convenient with the auxiliary apparatus 
on hand. 

The method referred to above consists of power 
detection by the bolometric technique. These tech- 
niques have been refined to the point [8] where 


is 


scheme the measurement is shown in 
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Ficure 5. Possible heterodyne detection system. 


| detector 3. 


accuracies of 0.1 to 0.2 percent may be realized at 
the required power levels. (It will be noted that 
questions as to bolometer mount efficiency of the 
detectors, etc., are unimportant here, since only 
power ratios are required.) It will prove useful to 
consider certain aspects of this procedure in detail 
since the y are also applicable to the techniques 
described earlier. 

The techniques by which the power is measured 
were discussed in detail in the cited reference [8] 
and need not be described here. In order to increase 
the signal at the detectors to the level where power 
detection feasible, the configuration shown in 
figure 6 is emploved. With the exception of the re- 
versal of the connections to the coupler P, the sys- 
tem and tuning procedure are the same as that given 
for figure 4. <A “erg nt choice of coupling ratio 
for the couplers P, Q, R, S, is 10 db. Under these 
conditions, nent esto d 80 percent-of the power 
reflected from the termination is delivered to the 
A further increase in this figure may 
be realized by other choices of coupling ratios, for 
example a figure of 90 percent may be realized by 
changing P and S to 3 and 20 db respectively, while 
a value of 98 percent would obtain with the values 
20, 20, 3, and 20 db for P, Q, R, and S, respectively. 
These two alternatives would require an additional 
3 and 10 db of power from the generator, however, 
and thus the first set of values appears to be a 
reasonable compromise when the available power is 
a potential problem. 

The power delivered to the detector 4 is stabilized 
or held constant by means of techniques developed 
in this laboratory, [9] and thus only the 
power at the detector 3 are required, 

In order that the detector signal may be as large 
as possible, it is desirable to operate the barretter, 
in the mount to calibrated, over as wide a re- 
sistance excursion as possible. The commercially 
available barretters are not all equally suitable in 


is 


values of 


be 
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this regard, but the type selected for this application 
may be operated over the range 140 to 260 ohms 
(nominal operating resistance 200 ohms). In op- 
eration the barretter forms one arm of a self-bal- 
ancing bridge [8], another arm of which is varied 
to produce the required changes in resistance. 

The bolometer mount under measurement is pro- 
vided with a nominal 10 mw of microwave power, 
resulting in reflected signals of approximately 310 
and 170 ww at the resistance values of 140 and 260 
ohms. <A practical problem exists with reference to 
the measurement at 140 ohms, since the 10 mw of 
microwave sing alone, without any additional d-c 
power, will bias the element at perhaps 180 ohms. 
This difficulty is avoided by an indirect measure- 
ment technique. In order to obtain the measure- 
ment corresponding to the 140-olhm resistance value 
the microwave power is first reduced to such a value 
that the element may be operated at 140 ohms. 
The signal in the detector 3 is then nulled or balanced 
out by means of the auxiliary signal channel. As- 
suming linearity of the system, this null will be inde- 
pendent of power level. The bolometer resistance 
is then returned 200 ohms, and the microwave 
power increased to its nominal 10 mw value. Under 
these conditions, the only signal delivered to detector 
3 will come from the auxiliary channel and will be 
equal in amplitude (but opposite in phase) to the 
signal which would be obtained from the bolometer 
mount were it possible to reduce its resistance to 
the 140 ohm value in the the 10 mw 
signal, This signal is then measured and used in 
the equations where the 140 ohm value is required. 


to 


presence of 


With the bolometer mount replaced by a short, 
the reflected signal has, of course, power level 
equal to that in the incident wave, nominally 10 
mw. The procedure thus requires measurements 
over the nominal 20 db range between the 10mw 
and 100 pw power levels. In addition to the meas- 
urement of signal levels, it is also important 
able to recognize the presence of signal levels 
in the region 60 to 80 db below these values during 
the course of the tuning operations. For this pur- 
pose it has proved convenient to split off a portion 
of the signal from arm 3 by means of a directional 
coupler and employ an susan: heterodyne detector 
A waveguide switch in the secondary arm of this 
coupler is kept in the nontransmitting position when 
the adjustment of tuner 7, is made and when the 
signal amplitudes are being measured by the bolo- 
metric detector. In other words, the switeh is in 
the transmitting position only when the high sensi- 
tivity for tuning adjustments is required. This 
isolates the remainder of the system from possible 
variations in the impedance of this heterodyne 
detector. 

Although little has been said concerning the de- 
tails of the power measurement, this aspect of the 
procedure not without its share of practical 
problems. In order to make power measurements 
to an accuracy few tenths of a percent at the 
100uw level, temperature control of the bolometer 
mount to about 0.001 ° C is required. A second 


a 


these 
to be 


Is 


of a 








problem concerns the barretter element itself. In 
practice it is frequently observed that even though 
the ambient temperature is held constant, the bias 
power required to maintain a given operating re- 
sistance may still be subject to ¢ onsiderable variation 
or drift. Certain types of elements have been ob- 
served to be much worse than others in this respect, 
but even among the better ones it is usually necessary 
to apply the bias power continuously for several days 
or even a week before stable operation is achieved. 

Although the foregoing discussion has been con- 


cerned primarily with the barretter element, the 
measurement also imposes stringent stability re- 
quirements on the associated apparatus as_ well. 
This will be easily recognized when attention is 


called to the fact that the bolometric measurement 
is in reality a differential power measurement. 
Thus if it is required to measure a one percent change 
to an accuracy of a part in 10°*, a stability or resolu- 
tion in the total d-c bias power of a part in 10° 
implied, 


IS 
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6. Approximate Methods 


The foregoing procedure is ‘fexact”’ in the sense 
that the only mathematical approximations which 
have been introduced thus far have been in assuming 
that eqs (1) and (4) adequately describe the bolom- 
eter mount and reflectometer behavior. In a large 
percentage of the cases of practical interest a con- 
siderable simplification in the technique may be 
effected with little loss in accuracy by the introduc- 
tion of suitable approximations. As a matter of 
fact, the accuracy achieved in practice is actually 
improved in certain instances by one of these approxi- 
mations which will be described. 


The first of these approximations relates to the 
substitution of ||T;3)+ 1,|| for |f;—T,| in eq (2), or 
more generally, the substitution of 

r3—Ts by Ee bs lr, Cy Fs 
a age 208 7 +x 
1—T,F,*|~ |1—T,r.*| r,r.* 1—Irilz 


ineq (6). (The minus sign is used if both R,; and R; 
are either greater than or less than Ry», while the 
positive sign is used when, as in the following ex- 
ample, one resistance value is above and the other 
below R..) 

In order to examine the validity of this approxi- 
mation it is convenient to obtain a gener ale xpression 
for this error in terms of the magnitude and argu- 
ments of the complex numbers T,, 3, or their more 


i—T; 
general counterparts — ete. Let A=|A| and 
iP. 
B= \Bie'*~®* represent two complex numbers. Then 


for small values of ¢, the error in replacing the vector 
difference by the sum of the magnitudes is given by 
the expression: 


Error= 


. (l—cos 9), 
(\A )2 % 
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The input reflection coefficient of the bolometer 
mount is given by the equation [2]: 


th. 
T=Sy i_S.T, (10) 
where I, is the load reflection coefficient, and the 
Sn.» are the scattering coefficients of the bolometer 
mount (see fig. 7). 

A substantial simplification in the analysis results 
from the fact that it is permissible to choose the 
impedance normalizing parameters in such a way 
that T,—0 when the bolometer is at its nominal 
operating value R,. For any other value of R, TY, 
has the value 


R—R, 


"ER, 


Utilizing these results, the term 


r.—T, 
i—T.r? 
becomes: 
Sbra 
n= T= SuF¥ gy 
1—T,r% ST Siz oo 


2 4 (Sat; —1Si rn 


Let V,. denote the argument of S,, and let 


F=( Sat aise (12) 


then, since the efficiency 7 of the bolometer mount 
when terminated in 2; is given by the expression [2]: 
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Ficure 7. Scattering coefficient representation of bolometer 


mount, 





It is shown in the appendix that |F| satisfies the 


inequality 
(16) 


In a large percentage of the cases of practical interest 
the efficiency will be in excess of 90 percent, while 
r,, is ordinarily limited to the range 0.1 to 0.2. 
Substituting these results in eq (9) yields the approxi- 
mate result: 


‘ 
Error «max 


(17) 


This is the approximate maximum limit of error 


which obtains when arg F 5 For other values 
of arg F the error will be smaller, becoming zero 
for arg F=0 or z. 

In a typical case where the bolometer is operated 
at the resistance values 140-200-260 ohms and for 
a mount efficiency of 90 percent the error is approxi- 
mately <0.012 percent. The uncertainty intro- 
duced by replacing the vector difference by the sum 
of the magnitudes in this example is thus wholly 
negligible. Under these conditions it will thus be 
recognized that the auxiliary waveguide channel, 
which was introduced to permit a determination of 
this vector difference, contributes virtually nothing 
to the accuracy of the final result as far as this 
aspect of the operation is concerned. 

The auxiliary loop is still a practical necessity, 
however, where power detection is employed because 
of the problem with the 140-ohm measurement as 
outlined above. In addition, this auxiliary loop may 
contribute to the accuracy of the result as follows: 

Let the resistance excursion of the bolometer ele- 
ment be chosen in such a way that Ty -T'3; then 
it can be shown that: 


Ll] T,—Ts rs) 6 (T'3—T2)(T,—T 2) 


(T,.—T,)(d 


—— (1S) 
T'.|*)| 


with an error 


Error < I'7,(1—n)? (19) 
which is approximately twice the value obtained in 
eq (17). 

One of the substantial sources of error attending 
this technique occurs in practice in the measurement 
of the signal amplitudes (corresponding to the various 
resistance values) at the detector 3. As outlined 
above, a total of four measurements are required, 
while use of the first approximation described reduces 
this number to three. Use of the approximation (18) 
however, reduces the number of required measure- 
ments to two, the vector difference and the measure- 
ment with the bolometer mount replaced by the 


The proof may be effected by substituting eq (10), (12), and (16) in (18) and 


simplifying the resultant expressions, 
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chort. In addition the power level which obtains at 
detector 3 in measuring this vector difference is four 
times that obtained from I » alone which reduces the 
dynamic range and improves the accuracy with 
which these measurements can be made. In the 
existing state of the art this latter technique is the 
most accurate, for mounts of high efficiency, of those 
described. 


7. Measurement of Attenuation 


The foregoing presentation has been directed pri- 
marily at the problem of measuring bolometer mount 
efficiencies. The efficiency, however, determines the 
dissipative component of attenuation, and a modifi- 
cation of the techniques described should also prove 
useful in other applications, such as measuring the 
efficiency of waveguide to coax adaptors, short sec- 
tions of waveguide, ete. A complete discussion of 
these potential applications and the requisite modifi- 
cations in the method is not, however, within the 
scope of this paper. 


8. Error Analysis 


In any technique or method with potential applica- 
tions in the field of standards work, a thorough in- 
vestigation of the attendant sources of error is an 
important part of the description. The sources of 
error to be discussed are the following: (1) Error due 
to misadjustment of tuning transformers, (2) error 
due to dissipation in reference short, (3) error due to 
uncertainty in the values of bolometer resistance, and 
(4) error in the measurement of the signal amplitudes. 
These will be treated in the order listed. 

The error analysis or study which has been made in 
connection with this technique, is, however, too 
lengthy to be presented in its entirety. For this rea- 
son the discussion to follow will be largely of a 
summary nature, and devoted primarily to the wave- 
guide system of figure 6. The extension of 
results to other variations of the method, however, 
should prove straightforward. The procedure to be 
followed in the analysis is that of assuming ideal 
operation of the system except for the source of error 
being considered. This technique yields the first 
order correction terms to the ideal theory for individ- 
ual sources of error. 


these 


8.1. 


Error Due to Misadjustment of Tuning 
Transformers 


a. Misadjustment of Tuner T, 


The adjustment of 7, and 7, is important in ob- 
taining the vector difference (T,—T,) ete. As a con- 
sequence of the adjustment of tuner 7, the signals 
delivered to detector 3 via arms 1 and 3 of directional 
coupler S (see fig. 6) may be considered as originating 
from two independent sources (of constant amplitude 
and frequency). If tuner 7’, is not properly adjusted, 
the signal delivered via arm 3 will vary with changes 
in T (as produced by changes in bolometer resist- 
ance). If the “null” observed at XY in step 1 of the 





tuning procedure is down from the signal at detector 
3 by at least 60 db, then the voltage amplitude of the 
signal emerging from arm 1 will be down by a factor 
of 10°. During the course of the measurement, T 
varies between the approximate limits of +1/7, and 
the maximum antic ‘ipated variation in the signal de- 
livered to detector 3 via the auxiliary channel due to 
this variation is approximately 3 parts in 10*. Finally, 
this is approximately one-half the total signal being 
measured, so the maximum anticipated error is one- 
half this value and may be neglected. 


b. Misadjustment of Tuner 7’, 


The problem once again is one of attempting to 
insure that the signal delivered to detector 3 via 
the auxiliary loop is independent of T. This time 
the possible interaction is through coupler R# instead 
of coupler S. In practice the variation in reflection 
coefficient observed “looking in” arm 3 of coupler 
P during the course of the measurement will be 
approximately +0.015. If with the sliding short 
connected to the output of coupler Q (tuning opera- 
tion 2), 7, is adjusted so that the maximum varia- 

| | 
tions in 7] are reduced to 1 percent (0.1 percent may 
[D4 
be realized in practice without too much difficulty), 
the maximum variation in the auxiliary channel 
signal will be of the order of a part or two in 10* 
and again may be neglected. 


c. Misadjustment of Tuner 7’; 


The function of tuner 7, is to adjust the four arm 
junction in such a way that it satisfies the condition 


B=—ATP,. The failure to exactly realize this 
adjustment may be accounted for by letting 
B=—A(T.+4). Substituting this result in the 
equation for determining the bolometer mount 


efficiency yields two approximate upper limits for 
the error due to this source, 


" ] 1iTn T Tis 7 l )| ° 
E< —}-_ 4 - 2\r.|+4(-—1 
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(20) 
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The choice of which expression to use obviously 
depends upon the relative magnitudes of T, and 


(!-1) 
S/] 

reduced by choosing the resistance variation such 
that Ty -T'3, and if the VSWR (c) of the mount 
is less than 1.02, the coefficient of |6; has the approxi- 
mate limit of 0.2 for an efficiency greater than 90 
percent. On the other hand, for a resistance excur- 
tion of 140 to 200 to 260 ohms and o=1.5 as in an 
earlier example, this coefficient will have the approxi- 
mate limit 2.8. The desirability of choosing resist- 
ance values such that [',——TI', is readily apparent. 


It will be noted that the error can be 


~ 
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If the approximate method implied by eq (18) is 
employed, the approximate limit of error is given 
by the expression 


’ rr’, » ‘. 
E<2| ++ r ta(l—a)yrA |. 6. 


In this expression the coefficient of |6, has the ap- 
proximate values of 0.06 and 0.5 for VSWR’s of 
1.02 and 1.5 respectively. The use of this approxi- 
mation thus reduces the error. 

If the tuner 7, has been adjusted such that the 
signal level at detector 3 is down by 60 db from the 
level which exists with the shorting plate connected 
in place of the bolometer element, |6| will have the 
nominal value of 107%, while an 80 db “null” will 
vield a value for 6) of 10-*. The error due to this 
source can thus usually be held to 0.1 percent or less, 


(22) 


d. Misadjustment of Tuner 7’, 


The function of tuner 7, is to adjust the junction 
such that it satisfies the condition D/C=(A/B)*. 
The failure to exactly realize this condition may be 
accounted for by letting D/C=(A/B)*+.. Sub- 
stituting this result into the equation for bolometer 
mount efficiency yields the approximate expression 
E= ¢ for the error due to this source. 
arises primarily in the determination of 
means of the fixed short. 

The failure to realize the proper adjustment of 7, 
can result from either failure to entirely eliminate 


This error 


A/D by 


b, ait 

, as the sliding short is moved, 
da 

or from deviations of the sliding short from the 
assumed ideal behavior. An analytic treatment of 
this problem yields for |e the approximate limit: 


e/<(1+2/r, | Fee 4 (i —/P.))ir.j+i—e | 


the variations in 


ener Fo... Pain, and P., 


represent the maximum, 


b,|? 

b,, * 
observed during the course of the sliding short 
adjustment, T,, is the reflection coefficient of the 
sliding “‘short”’ 
guide joint which connects the sliding short to the 
remainder of the waveguide system. In the above 
expression the waveguide in which the short slides 
has been assumed lossless, while in practice, when the 
correct adjustment has been realized, there will be 


minimum, and average values of the ratio 


¥ , . by 
a gradual and uniform decrease in the ratio js the 
sliding short is withdrawn. 

The error due to 7; may be eliminated by including 
a precision waveguide section in arm 2, as shown in 
figure 8, into which the sliding short may be inserted 
for this tuning operation. Other implications of this 
arrangement will be discussed in subsequent para- 


, and 7; is the efficiency of the wave- | 


PRECISION WAVEGUIDE SECTION 
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Figure 8. ZJllustration of use of precision waveguide section. 


graphs. It will be noted that the error due to the 

reflection coefficient of the sliding short (T,,) differing 

from unit magnitude decreases as the impedance 

match of the bolometer mount is improved. In 
° ° ° de” 

practice the variations in 7” can usually be reduced 

da 

to the order of several parts in 10° such that the 

error from these sources can usually be held to 0.1 

to 0.2 percent, 


e. Error Due to Imperfect Nulling When Using Auxiliary Channel 


If this ‘“‘null’’ is measured in terms of the signal 
level which exists in the detector prior to the nulling 
procedure it will be readily recognized that a 60 db 
minimum corresponds to a possible difference be- 
tween the two signals of a part in 10°, ete. Since 
this signal is next combined with another one of 
nominally the same amplitude and phase the error 
contribution to the final measurement is approxi- 
mately one-half this value. 


8.2. Error Due to Dissipation in Reference Short 


The error due to this 
with section (d) and is 
expression: 


source is closely connected 
given by the approximate 


(24) 


where [, is the reflection coefficient of the reference 
or standard short. In the present state of the art 
the highest reliable reflection coefficients are pro- 
vided by ‘‘quarter wave” shorts, that is, a quarter 
wavelength section of waveguide terminated by a 
soldered shorting plate, or the entire structure may 
be electroformed in a single piece. The advantage 
of this construction is in the elimination of the 
longitudinal component of current flow at the plane 


of connection, thus tending to minimize flange 
losses. Measurements and calculations of the re- 


flection coefficient magnitudes of such shorts [10] 
have vielded values in excess of 0.999. 

If such a device is used in conjunction with the 
uniform waveguide section as in figure 8, the meas- 
urement will include, part of the bolometer 
mount losses, the dissipation occurring at the bo- 
lometer mount input flange. If, on the other hand, 


as 
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bh. 
b, 
nated, there is some advantage in taking the average 
of the results with a quarter wave and conventional 
“flat”? short. This tends to cancel out the error due 
to imperfect tuning but introduces other uncertan- 
ties if there is an appreciable amount of flange 
loss [11]. The general features of this problem are 
discussed in greater detail in the cited reference. 


the variations in have not been completely elimi- 


8.3. Error Due to Uncertainty in the Values of 
Bolometer Resistance 





The bolometer resistances enter the expression for 
bolometer mount efficiency through the factor 


2h, R 
(RP, 


K . i 
Ri) (Rs R.) 

In practice the bolometer is biased by means of a 
self-balancing d-c bridge [8]. This device is capable 
of maintaining the bolometer resistance within 0.01 
percent of the value called for by the bridge param- 
eters. Use of NBS type resistors makes it possible to 
maintain the value of the reference arm within 0.01 
percent, while an error in the upper or ratio arms 
cancels out since the F’s may be multiplied by an 
arbitrary factor without changing the value of K. 
The F’s may thus be effectively maintained within 
the limits +0.02 percent. However, because the 
expression for A involves their differences, the total 
error due to this source is of the order of 0.1 percent. 
In addition the measurement will be in error by the 
ratio of the bolometer lead resistance to Ro. 


8.4, Error Due to Measurement of Signal Amplitudes 


An approximate limit to the error introduced in 
the measurement of the signal amplitudes is the sum 
of the errors in the individual measurements. For 
the approximation method of eq (18), this error can 
be kept within 0.15 percent if power detection is em- 
ployed, while for other variations of the procedure 
and/or other detection procedures the error will in 
general be larger. 

The error due to instability in the frequency of 
the signal generator has not been considered explicitly 
but in general will manifest itself in an inability to 
realize the tuning conditions to the required degree 
of precision. 

In summary, it will be recognized that the ultimate 
accuracy attainable will depend primarily upon which 
of the several alternate procedures is employed and 
upon the degree of refinement employed in the 


attendant adjustments and measurements. The 
foregoing results are somewhat arbitrarily sum- 
marized in the following table. The minimum 


values are intended to represent the approximate 
best values in the existing state of the art, while the 
larger values are intended to be typical of what may 
be achieved if the requirements are relaxed some- 
what. 
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Nominal error limits in percent 
Source of error 


o< 1.02 o=1.5 
Adjustment of 7 a 0. 00-0. 015 0. 00-0. 015 
Adjustment of 7, 2 . 00- . 02 .00— . 02 
Adjustment of 7, . 0O- . 02 . 05 3 
Adjustment of 7, : .05- .3 a; 5 
Imperfect null in use of 
auxiliary arm . 02— .05 02- . 05 
Reference short .05—- .3 . 08 3 
Bolometer resistance cn Ae aa 15 
Measurement of signal am- 
plitudes 15-1. 0 . 15-1. 0 
Total 0. 4 1. 8% 0. 5 -2. 3‘ 


9. Experimental Results 


An X-band waveguide setup employing the con- 
cepts developed in this paper is shown in figure 9. 
The efficiency values obtained utilizing this method 
have consistently shown agreement with microcal- 
orimeteric determinations [12] to within one-half of 
one percent and usually within a few tenths of a 
percent. This is of considerable interest since the 
two methods should, in principle, differ by the 
amount of the substitution error, and suggests That 
the substitution error is typically at least an order 
of magnitude smaller than the previously estimated 
limits. 


10. Appendix 
The inequality 
IFl<1—ny (16) 
can be proved in the following manner. 
Let the terminals of the 2 arm junction represented 


in figure 7 be reversed such that the load is con- 
nected to arm 1. Then the input reflection coeffi- 





FIGuRE 9. 
mentation employed in the measurement of bolometer mount 
efficiency. 


Picture of waveguide system and accessory instru- 








cient at arm 2 will be given by: (ef. eq™10) 


S2T: 
T'.=Sr 
, as a 


This can be written in the form, 


r al, +8 


= aT 4é (25) 
where a= S8},—S1S2 
B= So 
y=Su 
é=1, (26) 


Equation (25) indicates that Ty, is related to T, 
by means of a linear fractional transformation. This 
transformation has the well known property‘ of 
mapping circles into circles with straight lines as 
limiting cases. 

Let arm 1 be terminated by a sliding short. Then 
the locus of T, is the unit circle centered at the 
origin. From the above property of the transfor- 
mation it will be recognized that the locus of Tin 1S 
also a circle, and since a reflection coefficient mag- 
nitude can never exceed unity it is evident that this 
circle must lie completely within or on the unit circle 
as shown in figure 10. The distance to the center 
of this circle, r,, and the radius, r, thus satisfy the 
inequality: 


re<1—r. 


~ 


In terms of the constants of the 


~ transformation, 
r, and r are given by 


66*—avy* 
r ae (28) 
6/?—Iy/? 
and 
ai— By| 
r= 7". (29) 
0;|° =e - ; 


Substituting eqs (26) into (28) and (29) and these 
in turn into (27) and recalling the definition of F, 
yields the result: 


i. | ee 
4 So l<1— 
if : j i S|? | ” 
Sie 2 
where n=- -- - See ec (13 3) 
aoe Sul? | 


ee 


The author thanks R. W. Beatty and D. M. Kerns 
for their helpful suggestions in reviewing the manu- 
scripts, and W. E. Case, M. E. Harvey, and J. E. 
Gilbert who provided experimental de smonstrations 
of techniques described. 





4 See a text on complex variable theory 
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Telescope for Measurement of Optic Angle of Mica’ 
Stanley Ruthberg 
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’ 


The described instrument allows rapid measurement of the apparent optic angle to an 


accuracy of 5’ of arc for samples as large as 
pertinent to the quality of mica. 


») 


inches in diameter. 
Instrumentation is quite simple but dependent upon the 


This angle is a property 


complex phenomena of interference figures produced by biaxial crystals in polarized light. 
Magnification is great, dispersion can be determined, and the uniformity of samples can be 


observed 


l. Introduction 


The optic angle and the uniformity of mica split- 
tings are of interest, for these contribute to the 
quality of. the material, and are properties covered 
in specifications on visual quality.' 

The uniformity of the crystals can be determined 
with parallel, plane polarized light and a polarizing 
evepiece of polarization plane at right angles to that 
of the incident light. The specimen behaves as a 
wave plate with a resultant pattern dependent upon 
the variation of thickness and index. The orthoscope 
arrangement of the petrographic microscope utilizes 
such procedure, 

Direct measurement of the optic angle can be 
made with a number of instruments. For example, 
the conoscope arrangement of the petrographic 
microscope utilizes strongly divergent plane polarized 
light to produce an interference figure by which the 
angle can be determined with or without a universal 
stage, dependent upon the precision desired.2 The 
universal stage can be used to a tenth of a degree. 
A basic spectrometer, ie., collimator, goniometer, 
telescope modified with polarizing attachments, and 
rotary crystal holder can be utilized to greater 
precision. An instrument employing these com- 
ponents was developed for the precise location of the 
optic axis of polished quartz plates. Although the 
instrument and the procedure * were for the uniaxial 
quartz, the instrument itself could) be used for 
biaxial crystals. Such instruments as these are 
intricate and relatively expensive apparatus, which, 
further, examine only a very small region of the 
subject crystal. 

Rough measurements, at the other extreme, can 
be obtained with small apparatus of the type sug- 
gested in the ASTM. specifieations.° 

The described instrument, which has been used 
for a detailed study on mica behavior,® is relatively 
simple and readily assembled. It is used to deter- 


*Work supported by Defense Materials Service, General Services Administra 
tion 
American Society for Testing Materials, Tentative specifications for natural 
muscovite mica sed on visual quality, D351-53'1 
E.g., Manual of petrographic methods A. Johannsen, 2nd Ed., 1918 
(MeGraw-H Optical crystallography, E. E. Wahlstrom, 2nd Ed., 1957, 
(John Wiley and Sons 
E. Brodhun and O. Schénrock, Z. Instrumentenk, 22, 353 (1902 
‘E. Gumlich, Wiss. Abhand. der Physik.-tech. Reichanstalt 2, 202 (1895 
Appendix I, ASTM D 351-5311 
Results to be pul hed 


| 


mine both the uniformity of the crystals and the 
optic angle. In principle it is of the collimator- 
goniometer-telescope form except that the collimator 
and telescope are replaced with a primitive lens and 
the eve; yet, it has greater precision than the cono- 
scope and universal stage. It differs from the pre- 
vious instruments in that the entire area of a 2 in. 
diam specimen can be examined simultaneously to 
give crystal behavior in the large. Magnification 
of the interference pattern much greater than 
that of the conoscope; apparent optic angle can be 
measured to an accuracy of 5’ of are, and dispersion 
can be determined. 


is 


2. Theory and Structure 


In general the velocity of propagation of light in 
a crystal is a function of the vibration direction of 
the wave, or of the direction of propagation and the 
orientation of the plane of vibration. The optic 
axis of such a crystal is that direction for which the 
velocity of propagation is independent of the orien- 
tation of the plane of vibration. Then, in any other 
direction a plane polarized wave may be resolved 
into two components vibrating at right angles to 
each other with different indices of refraction. Many 
crystals, and mica is one, are biaxial. The angle 
subtended by the axes is the optic angle, labeled 
2V in figure 1. Waves propagated along the axes, 
A and B, are refracted upon emergence. The appar- 











Figure 1. 


Optic angle of a mica crystal. 








ent optic angle, labeled 2, is then that subtended | 
by the refracted rays. 

When divergent, plane polarized, monochromatic 
light of wavelength, A, passes through a_ biaxial 
crystal and then a polarizing plate, or analyzer, an 
intereference figure is produced which is composed 
of a family of contour lines of light of equal phase 
(isochromes) and a superimposed pattern of equal 
vibration direction (isogyres). When the direction 
of polarization of the analyzer is at right angles to 
the plane of polarization, when the optic plane of the | 
crystal is parallel to the direction of polarization of 
the polarizer or analyzer, and if the angle is symmet- 
rical about the surface normal, the resultant pattern 
is as in figure 2. Here, the dark maltese cross, 
A,-A,’ and P,’—P,, is the isogyre, formed by extine- 
tion of all light exit from the crystal and vibrating 
in the plane of polarization of analyzer or polarizer. 
The bright “circular” patterns are the isochromes, 
or contours of light of the equal phase difference of 
(2n+-1) 4/2 at the point of emergence between those 
components whose vibration directions are at right 
angles. 
of the family. 
section of the optic axes with the crystal surface. 
If the crystal is now rotated 45°, the isogyre breaks 
off into two hyperbolic brushes as in figure 3, with 
the optic axes centered in the vertices. The optic 
angle is then the are spanning the vertices, A and 
B, and is the property to be measured. 

White light 
various colors. In addition, if the optic angle and/or 
the orientation of the optic plane are frequency 
dependent, a dispersion pattern results for the 45° 
orientation. In this, the isogyres are separated into 
colored bands. In mica, the line of demarcation of 
the red and blue bands is sharp so that the inter- 
section of this line with the optic plane serves as a 
fairly precise reference point. 


An 


ISOGYRES 





(2n4i)472 
ISOCHROME 


Figure 2. Inte rference figure for optic plane parallel to one 
of directions of polarization. 


A,— A,’ for analyzer. P, A, B for optic axes 


P,' for polarizer, 


Points A and B represent the inter- | 


Each integral value of n describes a member | 


produces a system of isochromes of | 
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The structure of the “telescope” is as shown in 


cTOSS 


figure 4. A hair stretched across a large 
ground glass screen, illuminated from beyond by 
an ordinary incandescent light bulb, serves as an 
object and extended source. The crystal plate is 


examen ANALYZER 


: CRYSTAL PLATE 





POLARIZER 





OBJECT SPACE 


\ | 
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LIGHT SOURCE 
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FIGURE 4 Telescope 
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placed on a rotatable table equipped with protractor 
and vernier. The polarizer and analyzer are located 
before and after the table as closely as possible. A 
simple objective lens of large diameter is positioned 
approximately at its focal length from the object 
plane to produce normally a virtual image far beyond 
the object space. This can be done suitably in the 
following way. The crystal plate is rotated in its 
own plane and around a vertical axis in its plane 
until the 45° position is reached and a vertical trace 
of an isogyre is centered on the image of the cross 
hair. For this the eye is positioned above the 
analyzer in figure 4 and focused on the far off virtual 
image of the cross hair. The lens is then moved 
along the telescope axis until the image of the cross 
hair and the isogyre are parallax free. This proce- 
dure achieves a source of small divergence, which in 
turn produces an optic figure equivalent to great 
magnification. The line of demarcation of the red 
and blue fringes of the isogyre can now be fairly 
accurately centered; however, a vertical blade sight 
is also incorporated for greater precision, Aline- 
ment of the sight, isogyre, and cross hair insures 
that only axial rays are used and, hence, that the 
image is precisly centered. This also compensates 
for distortion by the large, simple lens and by crystal 
irregularity. The spanning of the specimen with 
polarizer and analy zer sheets of suitable diameter 
with close spacing and the use of a lens of large 
diameter create a large entrance aperture which is 
necessary to satisfy the demand for a beam of plane 
polarized light of large cross section needed to 
determine uniformity of fairly large splittings of mica. 

The image secured with the telescope is so large 
that only a very small portion of the optic figure is 
seen. That portion about the optic axis fills the 
view, and the vertex of an isogyre is so extended as 
to appear as a vertical trace across the entire crystal. 
This facilitates accurate reading, and in addition 
shows any variation of the angle of the exit ray across 
the crystal, which variation does occur in mica. 
However, to some, the image of the isogvre so found 
may seem reversed in color sense from that as seen 
in the petrographic microscope. Actual measure- 
ment of the apparent optic angle for the red and blue 
fringes shows that this is not the case. The appear- 
ance of reversal is due to the difference in optics. 
That is, in the conoscope arrangement of the petro- 
graphic microscope the strongly diverging light 
passes up through the crystal as in figure 5 and is 
viewed from above. Let A, and B, represent the 
optic axes for red light, with A, and B, for blue light, 
and let the red optic angle be greater than the blue. 
Thus as one looks toward the surface of the crystal 
away from the light source, he sees the blue fringe 
outside the red fringe with respect to that surface 
normal, for red light propagating along a red optic 
axis is absorbed in the analyzer to leave a blue 
excess, and vice versa for the blue light at a blue 
optic axis. The whole figure may be seen at once. 
On the other hand, the telescope with its weakly 
divergent light, extended source, and limited pattern 
has the behavior of figure 6. Here, as one looks 
down the light from the axes, B, and B, in this 


BLUE 














FiGguRE 5. Dispersion of wsogyre with conoscope. 


1,, B, for red optic axes 1,, By for blue optic axes. 2V is optic angle. 2E is 
ipparent optic angle e>b 


BLUE 











FiGURE 6. Dispersion of tsogyre with telescope. 


R, is red optic axis, By is blue optic axis. r>b 








example, toward the surface away from the light 
source, he sees the blue fringe inside the red fringe 
with respect to that surface normal, and he has only 
the one isogyre in view. This may seem reversed 
at first, but as can be seen the effect is due to the 
extension of the rays from the image plane to the 
eye through the blue and red optic axes. Thus, as 
the crystal is rotated in either direction from the 
position of symmetry, the red fringe, in this example, 
appears at the cross hair first. 

The angle so measured is not the true optic angle 
but the apparent optic angle because of refraction 
of emergent rays. 

When the crystal plate is normal to the telescope 
axis and the eye is focused on the specimen, the 
apparatus serves as a large sample orthoscope utiliz- 
ing parallel rays. The sample acts as a wave plate 
to give a color pattern dependent on_ thickness, 
birefringence, and uniformity of the surface. 


3. Operation 


An apparatus is shown in figure 7. The objective 
lens is an ordinary 4-in. diam reading glass of 25 em 


focal length. The analyzer and _ polarizer are 
Polaroid sheets of about 4.5 em diam in rotary 
holders. The crystal holder and goniometer, which 


is the only critical part of the apparatus, is pictured 
in figure 8. Here, the specimen can be rotated in 
its own plane as well as around a vertical axis. 

In use, a crystal is placed in the holder. Analyzer 
and polarizer are crossed at vertical and horizontal 
attitude. One’s eye is positioned to focus on the 
crystal itself, which is rotated for maximum and 
then minimum transmission. This eye position 
utilizes the instrument for observation of uniformity. 





FIGURE 7. 


Ope rational tele scope. 
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FIGURE 8 


Sam ple holde rand goniometer. 


The eve is then moved in to the analy zer and focused 
on the image of the cross hair. This converts the 
instrument for measurement of apparent optic angle. 
For this the specimen is rotated until the optic plane 
is horizontal. The horizontal trace of the isogvre 
may be slightly displaced from the horizontal cross 
hair for mica. The table is rotated to bring one of 
the optic axes into view. The Polaroids are then 
rotated 45 until the vertices of the isogvres are 
vertical. The eve is positioned so that the cross 
hair, sight, and isogvre are centered. The sight also 
serves to restrict change in vertical position of the 
eve. The apparent optic angle can so be deter- 
mined. 

Flat reproducibility to better 
than the accuracy of the goniometer itself, which is 


~/ 


» Ol 


specimens allow 


: are. In rippled specimens the measure is 
affected by local curvature through consequent 
variation of surface normal and hence refracted angle. 

With the components used, the Polaroid sheets 
limited the field in the examination of uniformity, 
while the lens limited the field in the determination 
and nspection of apparent optic angle. 


( Paper 65C 2-63) 


128 


An Automatic Fringe Counting Interferometer for use 


the Calibration of Line Scales 


Herbert D. Cook and Louis A. Marzetta 


(January 31, 1961 


A reversible fringe 


counting interferometer is described in which mechanical, optical, 


and electronic adjustments are maintained stable by servomechanism control or by bal- 


ancing. 
correcting by 


means of barium titanate 


rates are possible over shorter ranges. 


Design factors and details are discussed. 


Mirror parallelism is achieved by detecting the angular error electronically and 
actuators. 
recording of the court in digital form to 0.01 fringe without ambiguity. 

,200 fringes per second has been achieved over a range of 14 centimeters. 


An electronic interpolator permits 
A rate of more than 
Higher counting 
A correc- 


tion factor is derived for the error introduced by finite collimation of the interferometer beam. 


1. Introduction 


As part of a program to take advantage of the new 
radiation standard of length, an experimental revers- 
ible fringe counting interferometer was developed 
for use in the calibration of line scales at NBS. At 
the present time, line scales are calibrated by com- 
parison with a material standard, Meter bar No. 27, 
the length of which was derived from the Inter- 
national Prototype Meter maintained by the Bureau 
International des Poids et Mesures at Sevres in 
France. . 

Calibration with respect to the wavelength of a 
suitable radiation such as the lines of Kr 86 or Hg 
198 has several advantages over the use of a material 
standard. The wavelength standard is stable, easily 
reproducible in any laboratory , and can be accurately 
compared with the international standard. By 
means of fringe counting and interpolating tech- 
niques, any arbitrary length within the range of 
interference can be determined with great precision. 

Since specially trained personnel are required for 
calibrating line scales and the work is tedious, the 
number of lines that can be observed in a working 
day is limited. As a result, some scales with many 
graduations can only be partially calibrated. The 
of a mechanized instrument incorporating an 
interferometer and a photoelectric Iicroscope for 
detecting the position of the lines 
expected to be of help in making 
calibrations. 


use 


of the seales is 
more complete 


In the use of an instrument, the operation of which 
requires Many exacting procedures, the accuracy of 
measurement depends to some extent on the skill of 
the operator. Mechanization of such an instrument 
reduces the dependence on the operator and in some 
cases permits remotely controlled operation, allowing 
closer control of the environmental conditions. 
Carrying the mechanization one step further, auto- 
matic control is important if full advantage is to be 
taken of mechanization of the instrument, and the 
operator freed from full time attention to the details 
of the calibration process. 








1.1. Problems Encountered in the Design of an 


Interferometer 


The principles of operation of fringe counting 
interferometers have been extensively described in 
the literature. However, problems still exist in the 
application of the principles to a practical instrument. 

The basic problem in measuring length with an 
interferometer is to determine the fringe order 
accurately as the position of one mirror changes with 
respect to the other. The number of fringes counted 
between two positions of the moving mirror is the 
number of half wavelengths of light in the displace- 
ment. If the displacement corresponds to a length 
to be measured, the value of theintegral and fractional 
fringe count then represents a measurement of the 
length in half wavelengths. The fringes must be 
detected and analyzed with sufficient accuracy to 
allow inference to be made of the lengths repre- 
sented by them. Any factors that influence the 
accuracy of measurement must be either under 
control during operation or be measured and recorded 
at the time. 

Among the more serious problems is that of the 
mechanical precision required in the mechanism for 
moving the mirror. In some applications, corner 
reflectors have been used to eliminate the effect of 
angular rotation of the moving carriage, but in this 
application it was considered that their use did not 
constitute a suitable solution of the problem. If the 
carriage rotates in an unknown manner, the relation- 
ship between the center of reflection and the item 
being measured is in doubt. It is believed that more 
confidence can be placed in the accuracy of calibra- 
tions if plane mirrors are used and their parallelism 
maintained by electronic and mechanical means. 

Methods must also be provided for overcoming 
errors which may arise with variations in intensity 
of the light source and for maintaining necessary 
adjustments correct throughout the period required 
for a calibration. 

1 See biblic 


graphy at the end of this paper. 
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The method of interpolation must be applicable to 
operation without manipulation by an operator and 
yet maintain a precision to the specified tolerance. 
A sine-cosine method of interpolation requires that 
the amplitudes and phase angles of the quadrature 
fringe signals be held always within close limits. 

Solutions to these problems have been incorporated 
in an experimental instrument to be described in 
this paper. This was built to prove the design for a 
larger full scale instrument to be used in the length 
calibration laboratory. 


2. Optical Considerations 
Sa 


Interferometer 


A diagram of the optical arrangement is shown in 
figure 1. The Michelson type of interferometer is 
especially suited to the fringe counting application 
because the fringe intensity varies sinusoidally with 
mirror separation. Also, the range of travel of the 
moving mirror extends continuously on both sides of 
the zero-order fringe, allowing a total range double 
the maximum mirror separation at which fringes can 
be obtained. 

Circular fringes are used in order to obtain the 
maximum possible light flux through the interfer- 
ometer and to establish a mechanical symmetry in 
the system. Efficient use of light is important 
because of the low source intensity and the small 
apertures required at large mirror separations, the 
intensity being a limiting factor on the counting 
distance and the counting rate. 


Two aperture plates are shown; the first at the 
light source determines the maximum angle, @nax, 
at which light may pass through the interferometer. 
The second, placed in the image plane, has larger 
openings and is used to limit any stray light that 
might reach the phototubes and also to act in con- 
junction with a dispersing prism to form a mono- 
chromator. 

The limit (@,.) of the solid angle of rays passing 
through the aperture and the interferometer depends 
on the maximum mirror separation at which reliable 
counting is required. At greater separations, more 
than one order of interference occurs in the light 
reaching the phototube, and the fringe contrast is 
seriously reduced. The radius of the maximum 
allowable aperture may be taken as approximately 


teu te 
: J Jy | 


where Dmax iS twice the maximum mirror separation, 


f is the foeal length of the collimating or objective 


lens, and J is the wavelength of the light. This rela- 
tionship is derived in appendix 2 in addition to a cor- 
rection factor for an error in the measurement of 
length due to the use of a finite aperture. Accurate 
centering of the aperture is also required. 

The beam is divided into four parts, each quadrant 
producing a separate set of fringes which are detected 
by a separate phototube. This divides the inter- 
ferometer into essentially four separate interfer- 
ometers using portions of the same mirrors. The 
signals from the four phototubes are combined in 
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different ways for use in parallelism control, stabili- 
zation, counting, and interpolation. 

In order to take advantage of the 4-phototube 
arrangement and to obtain the sine and cosine 
voltages necessary for reversible counting and 
interpolation, the optical path difference in each 
qui adrant is adjusted so that a displacement of one- 
quarter fringe exists between adjacent quadrants. 
Beginning with one quadrant as a reference, the 
fringe phases are respectively, in rotation, 0°, 90° 
180°, and 270 

The adjustment of phase is made by means of a 
retarding plate placed in each quadrant, the effective 
path length varying with the angle of tilt of the plate. 
Very close adjustment of the relative phase of the 
fringes can:tbe made in this manner. 

In the form of Michelson interferometer 
fringe contrast 
order because of an extra reflecting surface in one 
beam. The contrast is decreased in proportion to 
the ratio of the intensities of the two beams and 
this depends on the square of the coefficient of 
reflection of the 45-degree mirror. For this reason, 
the efficiency of this mirror is more important than 
that of any other single part of the optical system. 
A coefficient of 0.80 for an aluminized 45-degree 
mirror gives a maximum contrast of 0.64. 

The effect of less-than-unity-fringe contrast on the 
fringe signal output of the phototubes is to decrease 
the signal-to-noise ratio by approximately the same 
factor as the contrast. A dielectric coating on the 
45-degree mirror with a coefficient greater than 95 
percent is possible and would make the effect 
negligible. 

Although the intensity of one beam is reduced by 
an extra reflection, that of the other is also reduced 
by the retarding plates, tending balance the 
beam intensities and to increase the contrast and 
signal to noise ratio slightly. 


used, the 


is always less than unity at zero 


to 


2.2. Light Source 
A suitable light 


source for measurement by inter- 
ferometry is 


which has high intensity and a 
spectral line width as narrow as possible, since the 
length over which interference fringes can be 
obtained (the coherence length) varies in reciprocal 
relationship to the width of the spectral line. The 
light sources available today have a coherence length 
of somewhat less than a half meter, but improved 
sources expected to be available in the near 
future. 

Although the orange line of Kr 86 is the primary 
standard of length, the Meggers Hg 198 lamp is 
probably a more suitable source for high- speed fringe 
counting sinee the intensity of the lines is higher. 
Also it is more convenient to operate. 

The Hg 198 electrodeless discharge lamp is excited 
by a high-frequency generator at a frequency which 
should be higher than 200 Mc/s. A microwave source 
operating at 2,500 Me/s has been convenient. The 
lamp is used end-on and is cooled in order to reduce 
Doppler broadening and self-reversal the lines. 
Immersion in flowing water is a simple method 


one 


are 
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cooling, but absorption of the exciting energy makes 
coupling from the source inefficient. 


Immersion in a 
jacket with recirculated, 


cooled oil has been quite 
satisfactory. For experimental work, selection of the 
green or blue lines was made by means of pairs of 
Wratten filters, 58 and 77A for the green, and 47B 
and 34A for the blue, but considerably more light 
intensity can be obtained if a prism monochromator 
is used, 


2.3. Phototubes 


In order to make use of the interference fringes 
generated in the interferometer, it necessary to 
convert the light flux into an equivalent electrical 
signal by means of a phototube. Since phototubes 
are available with operating characteristics covering 
a wide range, a few words might be said about the 
factors involved in the selection of a tube for the 
application. 

There are two principal sources of noise within the 
phototube, noise associated with the presence of light 
at the photocathode and noise generated by the so- 
called dark current. A typical photomultiplier tube, 
in total darkness, can have a dark current flow due 
to thermionic and other types of electron emission 
that is equivalent to the current that would result 
from a radiant flux of about 8107" lumens when 
operated at room temperature. Measurement with 
the Michelson interferometer used in the described 
equipment show that the light flux at the phototubes 
is approximately 1.5><10~* lumens when an aperture 
is selected for a measurement range of about 15 cm. 
Since the dark current is at least one order of magni- 
tude below the current due to this value of light flux, 
it can be neglected, as the signal noise (Ns) and the 
dark current noise (Np) add as yN+N,?. If the 
light source were less intense, or the aperture reduced 
in order to obtain a longer maximum interference 
path length, the dark current would become impor- 
tant. 

The noise associated with the presence of light flux 
is similar to shot noise in an electron tube and is due 
to the finite number of electrons per unit time in the 
current. Although the average phototube current is 
constant for a given incident light flux, the instanta- 
neous current has a random variation about the 
average due to the noise. It can be shown that this 
noise current varies as the square root of the electron 
current,? whereas the electron current varies directly 
with the light flux. Therefore the signal-to-noise 
ratio is proportional to the square root of the light 
flux. 

The magnitude of the noise that can be tolerated in 
the fringe signal must be viewed from a consideration 
of the statistical behavior of the instantaneous noise 
pulses. In the length of one meter, the number of 
fringes to be counted is in the vicinity of 4,000,000, 
depending on the wavelength. If the noise compo- 
nent amplitude is assumed to have a normal distri- 
bution, a calculation will show that the signal-to- 
noise ratio must be greater than 6:1 if the probability 


is 


Schottky’s theorem, W. Schottky, Ann Physik Vol. 57, 541, (1918). 








of a noise pulse exceeding the signal in amplitude is to 
be less than one in 510° fringes when the counter is 
operating at the maximum rate. The input circuit 
used with the counter described later reduces the 
probability of miscount under these conditions. 

The quantum efficiency of a photocathode is the 
ratio of the number of electrons released to the num- 
ber of photons incident on its surface. The efficiency 
of selected S11 type cathode materials may be of the 
order of 10 to 15 percent. Recently a multi-alkali 
photocathode has been made available with an effi- 
ciency near 20 percent at 4200 A. Furthermore, this 
new material offers a broader response curve over the 
visual spectrum plus an improvement in efficiency of 
3 to 6 times over conventional photocathodes in the 
red region. Since the orange line of krypton lies in 
this region, this tube will prove valuable for use with 
the standard light source. 

Measurements with the present equipment have 
indicated a potential due to the light flux of about 1 v 
root mean square at the photomultiplier anode. 
Since a level of about 5 v is needed at the trigger cir- 

cuits of the counters, it may be necessary to prov ide 
a d—c amplifier at the output of each phototube. Th 
10-stage multiplier tube used in the experimental 
equipment has a current amplification of about 600,- 
000. <A 14-stage photomultiplier operating at a 
higher supply voltage would have an amplification 
approaching 910°, eliminating the need for a d 
amplifier, but at a considerable increase in the tube 
cost. 


A small shift in the d—c level at the phototube out- 


























conditions, withno change of light level2* This can 
be reduced by operating the tube at reduced anode cur- 
rent. In the described system, the phototube anode 
load resistor is of the order of 5 megohms, and with 
the assigned voltage source, the anode current is less 
than 10 wa. An attempt to further lower the current 
by increasing the load resistor would result in a 
reduction of the response bandwidth. It was found 
to be essential to protect the phototubes from expo- 
sure to light levels that are in excess of their normal 
operating incidence flux. Fatigue effects of this sort 
usually persist for minutes or hours depending on 
the level of exposure. Because of the low level of 
flux present in the fringe pattern, the difference from 
normal operation to complete darkness that prevails 
with the equipment off does not result in an appre- 
ciable fatigue effect. 


3. Electrical Circuits 


A block diagram of the electronic circuits is shown 
in figure 2. The functions of the circuits include 
combining the outputs of the phototubes to produce 
stabilized sine and cosine signals, counting, inter- 
polating, and digital readout of the fringe count 
data. Also, the information necessary for correction 
of parallelism is derived and applied to actuators 
which maintain the moving carriage always in correct 
orientation. Remote and monitoring of 
operations are also required of the electronic system, 
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3.1. 


The fringes in the interferometer are sensed as 
fluctuations in the light flux reaching the phototubes. 
The flux has the form of a fixed value around which 
the intensity varies as a sine function of the optical 
path difference. The amplitude of variation depends 
on optical factors which determine the fringe con- 
trast, such as the flatness of the optical surfaces, 
reflection ratio of the beam splitter, the optical path 
difference and the spectral line width of the source. 
Since the fringe information lies in the variation only, 
it is necessary to isolate the signal in the form of a 
sinusoidal function with a baseline of zero. This is 
done conveniently by making use of the 4-phototube 
arrangement. 

By the effective thickness of the four 
plates in the interferometer beam, the phase retarda- 
tion Is adjusted so that there is an advance of one- 
quarter fringe between adjacent quadrants in rota- 
tion. Diagonally opposite quadrants are thus in 
phase opposition. Taking the electrical difference 
between diagonal phototubes results in a signi al 
which has twice the amplitude of each and a baseline 


Stabilization of the Electrical Fringe Signal 


Cc hoice of 


of zero. The baseline remains zero with variations 
in the intensity of the light source, although the 
amplitude of the resulting sinusoidal voltage is 


ie sadageo al to the intensity. A further discussion 
is given in appendix 1. Two outputs are available, 
one sae each pair of diagonal phototubes, 
representing the sine function of 
other the cosine. These are used to operate the 
bidirectional counter and the interpolator. The 
signal-to-noise ratio in the output from each of the 
diagonal pairs is y2 times that for each phototube 
since the signals are added coherently and the noise 
noncoherently. Precise adjustments of phasing can 


one 
distance and the 


be made easily to within 1/400th fringe. The accu- 
racy of the sine and cosine voltages thus obtained is 
such that interpolation to better than 1/100th fringe 
is practical by means of sine-cosine methods. 

Cireuit Details. As shown in figure 3, the fringe 
signals generated at the four phototubes are directed 
to individual amplifiers via cathode followers. The 
combination of anode resistance and circuit capaci- 
tance establishes a bandwidth of about 1,200 ¢/s. 
The sums and differences are taken for the diagonal 
pairs and BD. A-C represents the sine fringe 
signal, and B-D represents the cosine fringe signal. 
The sums A and B+D are used for gain level 
adjustment. 

To provide for balancing and adjustment of the 
photomultiplier gains, the sums and differences of 
each pair of diagonal phototube outputs are indicated 
on zero-center meters. The four meters have, 
associated with each, a variable resistor whose opera- 
tion influences the supply voltage of the photo- 
multiplier tubes in a manner such that its related 
meter is nulled. This is done while the interferom- 
eter is in operation. 


3.2. Reversible Counter 


The fringes sensed by the phototubes between two 
positions of the moving interferometer mirror are 
counted by an electronic counter and the total indi- 
cated as an integral number. 

The counter must be reversible in order to insure 
that the fringe count is true at the end of a measure- 
ment, at which point the motion of the interferom- 
eter may be reversed to correct an overshoot. 
This may occur, since in the intended application it is 
convenient if the interferometer is able to stop at 
each measurement. Other requirements of a suit- 
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able counter include a sufficiently high counting 
rate, convenient count indicators, and provision for 
automatically recording the count. 

These requirements have been met in a counter 
using neon glow-transfer tubes. A counter was 
built for the experimental instrument because a 
suitable one was not commercially available. The 
counter has 6 reversible decade stages, allowing a 
count of 10° without repeating. This exceeds the 
number of fringes in the range of the interferometer. 

Counter Circuit Details. Figure 4 shows a block 
diagram of the counter. The sinusoidal phototube 
output signals are first transformed into rectangular 
waves by means of voltage sensitive trigger circuits. 
The pulses which operate the counter in the forward 
or reverse direction are generated by differentiating 
the rectangular waves and are distributed to the app 
or SUBTRACT counter inputs by the use of gating 
circuits. One of the fringe signals operates a pair 
of gates which allow transfer of pulses derived from 
the other fringe signal to the respective counter 
inputs depending on the polarity of the derivative 
pulse at the time the gates are open. As shown in 
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| figure 3, the upper input becomes the gate signal and 
| the lower one the count signal for the waveforms 
| shown; the positive count pulse occurs during a 
positive gate signal. This common polarity operates 
| the forward aANp gate which in turn provides a pulse 
| to the amplifiers for registering an App count at the 
| glow-transfer tubes. The waveforms for reverse 
| operation are such that the positive gate signal 
| coincides with the negative count pulse. <A polarity 
| inverter reverses the negative pulse to positive and 
| delivers it to the reverse AND gate. A pulse then goes 
| to the amplifiers to register a SUBTRACT count. 

| ‘The second stage shown in figure 4 (as well as all 
| subsequent stages), possesses the reversibility feature. 
| In order to provide these stages with the proper 
| signals for reversible counting, 4 of the possible 
10-state positions from the previous counter stage 
are used as input For counting in the 
forward direction, the sequence of the four states is 
8, 9, 0, 1. Proceeding in a forward direction, the 
activation of state 8 at the first stage counter tube 
causes an output pulse to operate the lower flip flop 
into the position that renders the reverse AND gate 
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inoperative. On the 9th state the upper binary 
switch functions and places the forward ANp gate in 
a receptive position. The 10th or 0 state pulse 
passes through the forward ANpb gate and is registered 
at the counter tube in the manner previously de- 
scribed. The same 0-state pulse reverses the lower 
flip flop, placing the reverse gate in its receptive 
position. If the first counter tube should happen 
to reverse back to the 9th state, a pulse from this 


state would pass through the reverse AND gate. The 
counter would reverse or subtract out one count. 


The complete arrangement is compatible to the 
extent that all counters can be made to oscillate in 
synchronism between 9 and 10 by the application 
of a single forward or reverse pulse at this unique 
combination of counter positions. 

It may be worth while to say a few words about 
the glow-transfer counter tube. This electronic 
device has the ability to count input pulses and 
also visually indicate the totalized count. Its 
capacity is 1 decade, and it has provision for bringing 
out each of the 10 states for data transfer purposes 
and recording. In ~ for the tube to count, the 
input pulse is applied to 1 of 2 transfer electrodes. 
Connected across the unaaden electrodes is an induc- 
tance-capacitance delay network whose function is 
to direct the same input pulse to the second transfer 
electrode, but with a time delay of several micro- 
seconds. It is the application of this time-staggered 
pair of pulses that causes the ionized-gas glow to be 
transferred from one state to the next. The transfer 
electrode requires a fast negative pulse with an 
amplitude of about 100 vy. It is for this reason that 
a thermionic tube is used in the cireuit in lieu of a 
transistor, with its limited voltage swing. 

It was pointed out in the section on phototubes 
that, from a purely statistical viewpoint, a minimum 
value of design signal-to-noise ratio is established 
by the need to avoid miscounts. The probability of 
a noise pulse exceeding the fringe signal amplitude 
is one in 5108 counts at the maximum rate for a 
signal-to-noise ratio of 6:1. With higher signal-to- 
noise ratios the probability of noise pulses exceeding 
the signal is negligible. It was further mentioned 
that the probability can be lowered by the use of a 
special trigger circuit. The relative noise immunity 
is achieved by a modification of the Schmitt * trigger 
circuit increasing the hysteresis value. The hys- 
teresis value is the voltage difference between the 
level at which the trigger circuit changes state in one 
direction and that for a change in the opposite direc- 
tion. If the triggering levels are adjusted to be 
equally positive and negative from zero, and these 
are also at about half the peak voltage amplitude 
of the fringe signals, an optimum discrimination 
against noise results. In this case a rapidly repeated 
chang e of state which would be beyond the e: wp: ability 
of the reversible counter ean not occur unless a noise 
pulse exceeds the peak amplitude of the signal. 
If a small value of hysteresis were used, many rapid 
changes of state would occur with a small value of 
noise. 

Voltage sensitive 
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er circuit 





3.3. 


If the precision of the fringe count is to be greater 
than +43 fringe, a means is required for determining 
the fractional fringe beyond the last integral count. 
In determining the count, some reference point on 
the sinusoidal function representing the fringe must 
be chosen. This is usually taken to be one of the 
points at which the function equals zero. In the past 
a number of methods of interpolation have been used, 
including a means of changing the effective path 
length of the interferometer until the fringe fiducial 
point is reached. The fraction is then read on a 
scale which is calibrated by changing the path by 
one whole fringe. In another method, sine and 
cosine voltages from the interferometer are displayed 
on a cathode ray tube, the two voltages being put on 
separate deflection axes. The display spot then 
follows a circle with the advance of the fringe, and 
the fringe fraction is equal to the polar angle measur- 
ed around this circle, since there is one revolution per 
fringe. In order that the fraction may be accurately 
determined by this method, the voltages must 
accurately represent the sine and cosine values, the 
quadrature being accurate to a few degrees. Both of 
these methods are difficult to apply to automatic 
operation. 

Since automatic data recording is a requirement 
for the application, development of an electronic 
means for determining the fringe angle was under- 
taken. The result is a relatively simple analog 
computer which indicates the angle in digital form 
to 1/100th fringe in 1 millisecond. The decimal 
fraction is indicated on counters of the same type 
as are used for the integral count, so that the method 
of data recording is the same for both. 

The electronic operations in the conversion process 
may be described as taking place in three steps. Let 
¢ represent the fringe angle; ¢=360° for a displace- 
ment of one fringe. First, the two input voltages 
Asing and Acos@ from the phototube circuits are 
used to determine the amplitude of quadrature com- 
ponents of an alternating voltage with an arbitrary 
angular frequency w,. The result is Asingsinwt and 
Acos¢cosw;t. Second, the difference between these 
voltages is taken. The result can be expressed by a 
trigonometrical identity: 


Interpolation 


Asingsinw,t — Acos¢cosw,t = Acos(w,t+¢). 

The third step is to measure the phase shift between 
this sum and sinwf, since the phase difference is 
equal to the fringe angle ¢ 

Interpolator Circuits. A block diagram of the inter- 
polator is shown in figure 5. 

Electronic balanced modulation or multiplication 
is generally difficult and subject to errors in linearity. 
A technique which is believed to be new greatly 
simplifies the generation of the products required in 
the first step. The method makes use of linear gates 
in the resonance paths of LC networks. The 
principle of generation is shown in figure 6. 

Two resonant LC circuits, one for each of the 
quadrature signals, are excited to oscillation simul- 
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pared with zero by means of a voltage sensitive 
trigger circuit, and the time interval to the first 
positive zero crossing measured. The time intervals 
measured by the use of an alternating voltage of 
angular frequency 100 w; and 2 decade cycle counters 
which indicate the fringe angle in decimal form. 

Operation of the interpolator is not necessary 
during the process of counting, but is ne eded only to 
determine the fractional fringe at either end of a 
measurement. One initiating pulse at each position 
is sufficient to read the fraction, although it is con- 
venient to apply repeated pulses at perhaps 15 or 
20/see when watching the count visually. The ap- 
pearance that of a continuous reversible fringe 
count reading to 0.01 fringe. 


Is 


3.4. Ambiguity in the Fringe Count 


Since there is hysteresis in the trigger circuits of 
the integral counter, the purpose of which is to 
eliminate multiple triggering at zero crossings of the 
fringe signals, the indicated count does not change 
exactly at the fringe fiducial point used as a reference 
for the interpolator. This results in an ambiguity 
of a unit count in the region near the zero point of 
each fringe. 

In order to insure that the same fringe zero refer- 
ence point is used for both the integral counter and 
interpolator, it is necessary to make use of information 
developed by the interpolator. In forward counting, 
if the interpolator indicates that the zero point is 
passed, the integral counter must be foreed to indicate 
the next higher count. In reverse counting, the 
reverse must be true when the interpolator indicates 
less than zero. This necessary only when an 
interpolation is made, since only the final count is 
recorded. 

The antiambiguity circuit consists of a set of pulse 
generators responding to the 0.8, 0.9, 0.0, and 0.1 
states of the second glow-transfer tube of the inter- 
poli itor counter. At the end of the processing time 
of the interpolator, a pulse is generated which 
negative if the interpolation indicates 0.8 or 0.9, and 
positive if the indication is 0.0 or 0.1) One such 
pulse is applied to the counter trigger circuit and 
superimposed on the fringe signal. The pulse forces 
the signal voltage to exceed the triggering threshold in 
the direction of the count, eliminating the ambiguity. 


Is 
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3.5. Interpolation Bandwidth 


The signal-to-noise ratio allowable in the input of 
the interpolator depends on the desired accuracy of 
interpolation. The value of the sine or cosine input 
voltage may not vary by more than +6 percent of 
the maximum amplitude if the reading is to be 
accurate to +1/100th fringe. Allowing for the sta- 
tistical character of the noise and the characteristics 
of the counter and interpolator, a ratio of approxi- 
mately 1000:1 is required between the bandwidth 
for counting reliably and the bandwidth for interpo- 
lating to 1/100th fringe. The maximum fringe rate 
when counting or interpolating has this same ratio. 
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3.6. Parallelism 


Measurement of length with a Michelson interfer- 
ometer is accomplished by moving one of the plane 
mirrors which is mounted on a carriage, along a set 
of precise ways. Angular errors of motion of the 
carriage affect the geometry of the interferometer, 
and these may amount to as much as a micron per 
centimeter. If the error greater than one-half 
wavelength across the mirror, the fringes are entirely 
lost and counting is impossible. 


is 


Cube corners used as reflectors maintain the fringe 
pattern in spite of angular errors, but for the appli- 
cation it was considered more desirable to maintain 
sui ‘lism by correcting the error. If this is done, 
plane mirrors can be used. Also it is made certain 
that the position of a reference point on the moving 
carriage on which the mirror mounted always 
maintains a fixed relationship to the effective plane 
of reflection so that length measurements may be 
made accurately. 


is 


In the instrument described here the angular errors 
are eliminated by means of a servomechanism which 
includes first, a means for detecting the error, and 
second, a set of actuators to provide the necessary 
correction. A voltage indicating the error is obtained 
from a comparison of the phase of the fringes detected 
at diagonally opposite corners of the mirror. For 
this purpose the interferometer beam is divided into 
4 quadrants with a separate phototube for each, 
thus dividing the instrument into essentially 4 sep- 
arate interferometers working simultaneously. If 
the mirror is in perfect alinement, there is no phase 
error detected, but if there is an angular error and 
one mirror quadrant leads its diagonally opposite 
quadrant, a phase difference occurs between the 
respective phototube signals. An error of 1/100th 
fringe is easily detected. 

The actuators alter the relationship between the 
mirror and the carriage by an amount nearly equal 
to the angular error of the carriage so that the mirror 
does not depart from parallelism by more than is 
necessary to produce an error signal. Barium 
titanate is used for the actuators in the experimental 
equipment, but is limited to a few microns in its 
range of motion. Small hydraulic acutators have 
been used successfully for this purpose. In 
addition to their simplicity, both types of actuators 
have an advantage in that there is a minimum of 
energy dissipated in their operation, and 
required to hold a fixed position. 


also 


hone is 

Since in this equipment the phase is determined 
from the relative timing of the axis crossings of the 
phototube signals, the mirror must be moving at 
the time of correction. It necessary, therefore, 
to incorporate some means of maintaining the cor- 
rection voltage or pressure if the parallelism is to be 
held for long periods at rest. A suitable voltage 
memory is a motor-driven potentiometer controlled 
by the error signal. In the case of the hydraulic 
actuators, flexible tubes containing a fluid are raised 
»r lowered to provide the required pressure. 


is 





Circuit Description. In figure 7, A and C represent 
one diagonal pair of phototube outputs, and B and 
D the other. The discussion for one pair will hold 
equally well for the other. Each of the amplified 
phototube outputs is applied to a voltage sensitive 
trigger tube in order to generate square waves, the 
transitions occurring very close to the zero crossings 
of the sinusoidal fringe signal. The networks fol- 
lowing the amplifiers form the sum of the square 
waves, A’ and C”’, to their opposite derivatives, 
dC’dt and dA’dt. In figure 7, the waveforms shown 
are for the condition in which the mirrors are out of 
parallelism in a direction such that signal C lags 
signal A by slightly more than the normal 180 deg. 
Note that for this condition, the C’+ (dA’dt) output 
has a sharp positive pulse on top of the pedestal. 
If the mirror were to be tilted out of parallelism in a 
direction opposite to the previously assumed direc- 
tion, the upper waveform, A’+ (dC’dt), would have 
the positive pulse on its pedestal. 

For a mirror travel restricted to one direction only, 
the circuits described would be sufficient to provide 











a signal for correction purposes. 
bidirectional travel, the correction sense must be 
reversed if the direction is reversed. This is done 
by means of a circuit in which gates are actuated 
by signals from the B and PD signal channels. As 
shown in figure 7, the gate outputs are combined by 
diodes to form two separate output lines. Signals 
appearing at each line are further amplified (in addi- 
tion to a polarity inverter in one line), and are 
recombined at the grid of a cathode follower. The 
positive or negative output of the latter (depending 
upon mirror tilt and travel direction) is used as the 
error control signal in the servo loop for maintaining 
the mirror position. 

As mentioned previously, the servo loop includes 
a motor-driven voltage control for the purpose of 
providing position memory of mirror parallelism 
while the interferometer is at rest. 


However, for 


The d-e output 


| of the control in turn is amplified to a potential of 
several hundred volts in order to drive the barium 


| 
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titanate transducer. 
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Mirror paralle lism control circu 


CONTROL TRANSDUCER 


The phase detector circuit is capable of resolving 
a phase difference of 0.01 fringe. With a mirror 
3-cm square and with green light, this is equivalent 
10-7 radians, or about 0.02 sec of are. 


4. Accuracy of Measurements 


Since a length measured by the interferometer is a 
comparison with respect to the wavelength of the 
light source, any variation in the wavelength in the 
interferometer affects the accuracy of measurement. 
The stability of the wavelength of the krypton lamp 
may be 1 part in 10°, whereas the mercury 198 
lamp is affected by the presence of a carrier gas by 
a few parts in 10°. The Hg 198 lamp is calibrated, 
however, and may be relied on to the extent that the 
argon pressure is known. The wavelength is given 
to one part in 10°. 


A greater variation in the wavelength in the inter- 


ferometer is due to the refractive index of the 
ambient air. The refractive index error is pro- 
portional to the density of the medium and is 


therefore affected by the temperature, the pressure, 
and the proportion of other constituents such 
CO, and water vapor. If the wavelength is to be 
known to one part in 10°, the atmospheric pressure 
must be known to approximately 0.0025 em of Hg, 
the air temperature to 0.01 °C, water vapor content 
to 0.25 mm partial pressure, and CO, to 1/100.4 


as 


5. Aperture Error Correction 


Another factor influencing the accuracy of measure- 
ment is characteristic of the interferometer. In an 
interferometer in which a phototube is used to detect 
the fringes, the number of fringes counted not 
exactly equal to the distance of mirror travel divided 
by one-half the wavelength of light, since the photo- 
tube receives light from more than the central fringe. 
The ratio of the count obtained to the true count is a 
constant factor which depends on the degree of 
collimation of the light in the interferometer. The 
difference is significant since the numerical value of 
the necessary correction is generally greater than 
one part in 10° and may total more than one fringe 
in a given measurement. A correction factor for 
the error can be applied which is fortunately a con- 
stant for a given angular aperture and ap- 
plies equally throughout the range of measurement. 
It is also the same for a fixed aperture for any wave- 
length of light. 

If the source were a point on the axis of the colli- 
mator, the fringes would truly represent the separa- 
tion of the interferometer plates, but since the bright- 
ness of the source is finite, the light flux obtained from 
a point on the axis would be zero. In order that the 
light flux through the interferometer be finite, the 
source must be broad, and as a result, essentially all 
of the rays that are used are nonparallel to the axis. 
The average path of rays that are seen by the photo- 
tube is thus at an angle to the axis and the average 
optical path difference is not equal to that normal 
to the axis. 


‘W. F. Meggers, C. G 
Am. 43, 339 (1953 
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Peters, Bull. BS 14, 697 (1918); B. Edelen, J. Opt. Soe. 
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Integration of the light flux received by a photo- 
tube in an interferometer has been described in a 
number of published papers.®° An expression for the 
fringe contrast versus optical path difference and 
aperture is derived in appendix 2. The magnitude 
of the length antaieen is given. 


6. Conclusion 


The interferometer described here is applicable to 
uses in which the fringe count and interpolation is 
required to be continuously available, as in systems 
where fully automated control is desired. 

Further improvement in the measuring ability of 
the instrument will be in the direction of extended 
range and in the speed of counting. Such improve- 
ment depends primarily on the light source, the 
important factors being the monochromaticity of the 
radiation and the flux per Angstrom of spectrum per 
unit solid angle from a source of unit area. Both of 
these factors are characteristic of the coherent sources 
now being developed elsewhere and which should be 
available in the not too distant future. 

With present light sources, there is an approxi- 
mately reciprocal relationship between the maximum 
range of fringe counting and the speed of counting, 
the acceptance angle of the collimator aperture being 
the variable parameter. This true within the 
limitations of (1), the length over which interference 
can be obtained, and (2), the maximum rate of the 
counter. 

In the application of the interferometer to scale 
calibration, it is used in conjunction with a photo- 
electric microscope which has comparable precision. 
Measurement of a seale which is longer than the 
maximum range of the interferometer is done by 
measuring sections of the scale, allowing the measure- 
ments to overlap with several lines in common. The 
error in adding the sections to derive the total is 
thereby minimized by the redundancy in reading the 
lines. 
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7. Appendix 1. Relation of Phototube 
Signals 
If the source intensity is J and the fringe contrast 


the currents in the phototubes A, B, D are 
respectively : 


Mm, 


k,T(1-+mz, sin ¢) 
kpl(1+mz, cos ¢) 
k-I(1—me sin ¢) 
kpI(1—mp cos ¢) 


When balance adjustments are made, E.. ke, ke, Ks 
become equal. 


The difference between A and C is then: 


A—C=kI(m, 


1m) sin @ 
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Am. 43, 505 (1953 Joseph H. Jaffe; J. Opt. Soc. 
Huntoon, A. Weiss, J. Opt. Soc, Am, 44,4 1954). G.R. Harrison, G. W 
J. Opt Am. 45, 112 (1955 "i R. Peck, J. Opt. Am. 45, 931. ( 4 
Bruce Acta 4, 127 (1957). M.V.R.K Seuuie Teen ties Am. 50,7 (1960). 
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Also 
B—D=kI(mg+mp) cos ¢ 


The result is a sine function and a cosine function 
of @ with baselines equal to zero and amplitudes 
proportional to J. 


8. Appendix 2. Calculation of the 
Correction Factor 


A ray from a point within the aperture leaving the 
collimating lens at an angle @ with respect to the 
axis of the system passes through the interferometer 
at thissame angle. Since the optical path difference 
in the interferometer varies as the cosine of this 
angle, the average optical path difference for the 
complete set of rays is less than that for a ray which 
is parallel to the axis. The total light flux reaching 
the phototube through the aperture and interfer- 
ometer is found by integrating the interference func- 
tion over the angular aperture. 

The point intensity of a ray parallel to the axis is 
proportional to 

2xD 


rIRCcos 


1 (1) 


where D is twice the mirror separation, \ is the wave- 
length, and m is a factor representing the fringe 
contrast for a perfectly collimated beam at an optical 
path difference equal to D. For any ray the ex- 
pression becomes: 


IrD 
1+ meos( 2 cos 6 


. 


) 


Integrating over a sphere from zero to 6,, the limit 
of the angular aperture, 





. = 2xD se : 
Flux= E +m cos ( ar cos 0)| sin 6 dé 
J 6 a. 
. 2oD . (2nD 
‘ sin ss —sin (; r —- cos 6; ) 
—cos 1m — i : 
CC 1 n nD 
(3) 
By trigonometric identity this is equal to 
. {2xD 1—cos 6 
sin( ——> a ‘) a 
—_— ee \ 2 eee 
COS Oj + <m™M nD cos — . > +) 
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This is a cosine function of D the optical path 
difference in which the argument is similar to that 
of the intensity function for an axial ray except for 
the factor (1+ cos 6,)/2. This is the correction 
factor by which the fringe count must be multiplied 
in order to obtain a true measurement of length in 
half wavelengths of light. It is approximately equal 
to 

6. GF 


4 48 


] 





The expression also shows the variation of the 
amplitude of the cosine function with respect to D, 
assuming no reduction in contrast for other reasons. 
The amplitude becomes zero when D=)/(1—cos6,) & 
2\/@. The maximum useful counting distance 
Dogx 18 considerably less than this, perhaps }/6,? 
where the amplitude is about half the maximum. 
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Selected Abstracts 


Heat treatment and properties of iron and steel, T. G. 
Digges and 8. J. Rosenberg. NBS Mono. 18 (1960) 35 cents. 
(Supersedes C495). 

This Monograph has been prepared to give an understanding 
of heat treatment principally to those unacquainted with this 
subject. To this end, the basie theoretical and practical 
principles involved in the heat treatment of iron and steel are 
presented in simplified form. 


Household weights and measure, NBS Misc. Publ. 234 
(1960) 5 cents. (Supersedes M39) 

The purpose of this card is to present in convenient form the 
weights and measures tables most useful for household pur- 
poses, together with associated weights and measures infor- 
mation of general interest. 


Transistorized building blocks for data instrumenta- 
tion, J. A. Cunningham and R. Hill, NBS TN68 
(PB161569) (1960) $2.00. 

The National Bureau of Standards has developed a number of 
modular transistorized digital circuits that have been used in 
automatizing many data recording and preliminary processing 
tasks encountered in its scientific operations. These versa- 
tile building blocks can be connected together systematically 
to form digital circuits that accept raw data from experimen- 
tal equipment and transpose these data into a form suitable 
for input to a high-speed electronic computer. These blocks 
can be used where: (1) data are produced in large volumes; 
(2) data taking is extremely fast, extremely slow, or extremely 
precise; (3) a need exists to minimize human error, tedium, 
and eyestrain; or (4) computation is extensive. This report 
describes 5 major packages and 8 special purpose packages. 


Thermal voltage converters for accurate voltage meas- 
urements to 30 megacycycles per second, F. L. Hermach 
and E. 8. Williams, Commun. Electron. AVEE, No. 49, 200- 
206 (July 1960). 

Thermal voltage converters, each consisting of a resistor in 
series with a thermoelement in a coaxial line, have been de- 
veloped for measurements of rms voltages of 1 to 200 volts at 
frequencies from 3 eps (eyeles per second) to 30 me (mega- 
cycles). An accuracy of 0.1 percent or better may be ob- 
tained by a.c.-d.c. transfer techniques up to at least 10 me 
and 0.2 percent at 30 me. 


The strength of ten structural adhesives at tempera- 
tures down to minus 424 °F., W. Frost, Proc. 1959 
Cryogenic Engr. ( onf. 5, 375-384 (Se pt. 2-4, 1959). 
The shear strengths of lap-joints bonded with commercially 
manufactured structural adhesives were determined at 
107 °F, 323 °F, and $24°F. The effect of rapid cool- 
ing bond strength was studied at 323 °F. Epoxy- 
phenolic adhesives on glass cloth generally showed the highest 
strengths over the range of test temperatures. 


on 


Maser frequency stability, R. C. Mocklerand J. A. Barnes, 
Proc. 13 Annual Frequency Control Symp., 583-595 (May 
12-14, 1959 

The effort of the Bureau of Standards in maser work has been 
directed toward the attainment of high frequency stability 
for extended periods of time. A maser stabilized frequency 
multiplier chain has been controlled to 1 10-" for periods 
of 6 hours. These preliminary measurements and tests in- 
dicate that the maser can be held within narrow frequency 
limits for periods of weeks with simple modifications. The 
masers have been employed to study frequency multiplier and 





crystal oscillator stability. The data have enabed us to 
make some predictions of the behavior of extremely long 
atomic beam machines. 


Some mechanical properties of magnesium alloys at 
low temperatures, R. F. Reed, R P. Mikesell, and R. L. 
Greeson, Proc. 1959, Cryogenic Eng. Conf. 5, 897-405 (Sept. 
2-4, 1959). 

Unnotched tensile tests and Charpy V-notch impact tests 
were conducted on six magnesium alloys at four temperatures: 
300 °K, 195°K, 76°K and 20°K. The magnesium alloys 
tested were wrought AZ31B-0, ZE10A—H11, HM21A—T6 and 
HK31A-0; extruded HM31A-F and sand cast HK31A-T6 


It was found that the magnesium alloys generally increase in 
tensile strength, decrease in tensile elongation and decrease in 
impact strength with a decrease in temperature. Photo- 
micrographs which illustrate the principal modes of fracture 
were taken of the fracture surfaces. 


The behavior in tension and impact of magnesium alloys at 
low temperatures is compared with that of other common 
cryogenic fabricating materials, considering weight as an 
important criterion. 


Report on the standardization of pH and related termi- 
nology, R. G. Bates and E. A. Guggenheim, Intern. Union 
Pure and Appl. Chem. 1, No. 1, 163-168 (1950). 

The standardization of pH has been achieved in Great Britain, 
in the United States, and recently also in Japan. In their 
essential elements, the three standards are in agreement; only 
differences of detail exist. It is the object of this report to 
summarize, for the information of all national groups in- 
terested in pH, the extent of standardization already achieved 
and to indicate how this area of agreement could perhaps be 
extended to other conceptions related to pH. 


Regulated power supply for instruments, W. V. Loenben- 
stein, Electronics 33, No. 48, 132 (Nov. 1960). 

This paper describes a simple, low-cost device which will 
maintain output current constant to within Mo percent 
when the input is subject to normal fluctuations in line 
voltage. The device was used to supply d.c. heater current 
for a thermocouple-type vacuum gauge. 


Experimental investigation of creep deflection of ex- 
truded and riveted I-beams, L. Mordfin and N. Halsey, 
NASA Tech. Note D-662 (Dec. 1960). 

Five beams of 2024-T4 aluminum alloy were tested in creep 
bending at 450°F. These were extruded I-beams, two of 
which were reinforced with riveted flanges. The beams were 
simply supported at the ends and loaded with weights at the 
third-points. 


The creep deflections were compared with theoretical values 
calculated by the methods presented in another report. 
Agreement is found to be very good, thereby providing a 
measure of substantiation to the theory. 


Examples are used to demonstrate various aspects of the 
creep bending problem, such as the effect of rivet creep on 
stress distribution and deflection. 


The measurement of thermal conductivity, D. G. Gin- 
nings, Book, Thermoelectric Materials and Devices, edited by 
I. B. Cadoff and E. Miller, Chapter 8, 113-132 (Reinhold 
Publ. Corp., New York, N.Y., 1960). 

This chapter covers material given in a talk in a Thermo- 
electric Materials and Devices course at New York Univer- 
sity in June 1959. The talk was intended for novices in the 
field of heat measurements, with emphasis on measurements 
at high temperatures. This paper is one of twenty presented. 








Analytical study of creep deflection of structural beams, 
L. Mordfin, NASA Tech. Note D-661 (Dec. 1960). 


A study was made of existing methods of analyzing the creep 
deflections of beams. It was found that none of these 
methods are directly applicable to structural beams of the 
types used in airframe construction. 


Using Popov’s elastic followup technique, which is based on 
the hypothesis that plane sections before bending remain 
plane after bending, the stress and strain distributions are 
determined for the creep bending of beams of nonrectangular 
cross section having unequal creep properties in tension and 
compression. It is shown how the deflections are obtained, 
both for pure bending and for nonuniform bending. The 
corrections required by transverse shear and by creep in the 
riveted connections are then developed. It is shown that 
the analyses are adaptable to beams fabricated from two or 
more materials. Finally, the effects of shear lag on creep 
deflections are discussed. 


Evaluation of ball bearing separator materials operat- 
ing submerged in liquid nitrogen, W. A. Wilson, K. B. 
Martin, J. Brennan, and B. W. Birmingham, ASLE. 
ASME Lubrication Conf., October 17-19, 1960 (Boston, Mass.), 
Am. Soc. Lubrication Engrs. Preprint No. 60 LC-4 (1960). 


The demand for the reliable operation of rotating machinery 
in eryogenic applications has resulted in a continuing pro- 
gram directed toward determining the life and durability of 
ball bearings operating at cryogenic temperatures. On the 
basis of past work in this field, i. was determined that the ball 
separator, when made from conventional materials, usually 
was the cause of failure. A program of screening separator 
materials with bearings subjected to actual operating condi- 
tions was undertaken. A number of different separator 
materials and arrangements were tested to determine their 
relative durability. The apparatus used provided a method 
of loading the bearing axially and included a means of meas- 
uring bearing torque, wear and load during operation. All of 
the tests of this screening program were conducted using liquid 
nitrogen. The results are presented in the form of bearing 
life, progressive bearing torque and progressive bearing wear 
for the various configurations and types of separator materials 
tested. Photographs showing the test apparatus and some 
test bearings are included in addition to a graphic presenta- 
tion of the test results. This test program proved that a ball 
separator made from filled PTFE material is far superior to 
the other separator materials tested in liquid nitrogen. 


Aircraft storage batteries, W. J. Hamer, AJEE Trans. 79, 
Pt. II, 1-11 (Sept. 1960). 


The characteristics and performance of lead-acid, nickel- 
cadmium, and silver-zine storage batteries of aircraft design 
are discussed. A standard test program is outlined. In this 
program the batteries are subjected to a cycle of tests with 
constant-potential charging interspersed at intervals; part of 
the program is conducted manually, another part automat- 
ically. Performance data obtained at 80°F and 0°F are 
given for the various types of batteries. Discharge charac- 
teristics of the three types of batteries at various rates and 
temperatures are presented. Charge characteristics of the 
three types of batteries on constant-potential charging at 
120 °F, 80 °F, 0 °F, and 20 °F are considered. Finally an 
outline of physical tests and areas for future research are 
presented. 


Stress-strain relationships in yarns subjected to rapid 
impact loading. Part VJ: Velocities of strain waves 
resulting from impact, J. C. Smith, J. M. Blanford, and H. 
F. Schiefer, Textile Research J. 30, No. 10, 752-760 (Oct. 1960). 
When a textile yarn is impacted transversely, strain waves 
and transverse waves are initiated, which travel outwards 
away from the point of impact. Each strain-wave front 
travels at constant velocity. A transverse-wave front, how- 
ever, changes its velocity each time it meets a reflected strain- 
wave front. This effect forms the basis of a new method used 
to measure strain-wave propagation velocities in a selection 
of textile yarns. Velocities obtained in this way are in good 
agreement with velocities obtained by two other methods. 





These velocites range in values from 1400 m/see for undrawn 
nylon to 5000 m/see for high tenacity rayon and glass fiber. 


A comparison of atomic beam frequency standards, R. 
E. Beehler, R. C. Mockler, and C. 8. Snider, Nature 187, No. 
41738, 681-682 (Aug. 20, 1960). 

A general qualitative description of atomic beam devices is 
given with particular emphasis on cesium atomic frequency 
standards. 

The various uncertainties in frequency measurements are 
discussed including the effeets of the uniform C field, the effect 
of a phase difference between the two oscillating fields that 
excite the atomic transition, and the effect of exciting an 
atomic transition with signals that are not monochromatic. 
The two cesium beam devices constructed at the National 
Bureau of Standards are described in some detail. The 
various tests are described for determining the uncertainties 
in the frequency measurements. The results of these tests are 
given. 

The accuracy of these two quite dissimilar machines is con- 
sidered to be 1.5 10-" and the precision of measurement 
is 210°". The latter figure represents the standard devia- 
tion of the mean for measurement times of a few hours or 
about 100 separate measurements. 
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An optical study of the boundary layer transition processes 
in a supersonic air system, W. Spangenberg and W. R. 
Rowland, Phys. of Fluids 3, No. 5, 667-684 (Sept.—Oct. 
1960). 

Thermodynamic structure of the outer solar atmosphere. 
VI. Effect of departures from the Saha equation on in- 
ferred properties of low chromosphere, 8S. R. Pottasch and 
R. N. Thomas, Astrophys. J. 132, 195-201 (July 1960). 

Polymer decomposition: Thermodynamics, mechanisms, and 
energetics, L. A. Wall, Soc. Plastic Engrs. Pt. I, 810-814 
(Aug. 1960); Pt. II, 1031-1035 (Sept. 1960). 

The role of surface tension in determining certain clay-water 
properties, W. C. Ormsby, Bull. Am. Ceramic Soc. 39, No. 
8, 408-432 (Aug. 1960). 
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Procedure for the determination of the noble metal content 
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Rapid frequency analysis of fading radio signals, J. M. Watts 
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On Stokes flow about a torus, W. H. Pell and L. E. Payne, 
Mathematika 7, 78—92 (1960). 
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3318, 1-9 (Aug. 1, 1958). 

Modern theories of materials, C. Truesdell, Trans. Soc. Rhe- 
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Pressure dependence of rotationally perturbed lines in the 
ultraviolet band spectrum of CN, H. P. Broida and S. 
Golden, Can. J. Chem. 38, 1666-1677 (1960). 

Neighbor interactions and internal: rotations in polymer 
molecules. IV. Solvent effect on internal rotations, S. 
Lifson and I. Openheim, J. Chem. Phys. 33, No. 1, 109- 
115 (July 1960). 

Improvements in radio propagation prediction service, W. B. 
Chadwick, Elec. Engr. 1-4 (Sept. 1960). 

Carbon resistance thermometry with mixed de and rf currents, 
J. J. Gniewek and R. J. Corruccini, Rev. Sci. Inst. 31, 
No. 8, 899-900 (Aug. 1960). 

Some SEAC computations of subsonic flows, P. Davis and 
P. Rabinowitz. Book, Bergman’s Linear Integral Operator 
Method in the Theory of Compressible Fluid Flow, by 
M. Z. Krzywoblocki, p. 148-172 (Wien, Springer-Verlag, 
Berlin, Germany, 1960). 

Spatial distribution of energy dissipated by fallout beta rays, 
A. E. Boyd and E. E. Morris, Health Phys. 2, 321-325 
(1960). 

Comment on a paper of Mori on time-correlation expressions 
for transport properties, M. 8. Green, Phys. Rev. 119, 
No. 3, 829-830 (Aug. 1960). 
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M. M. Davis and M. Paabo, J. Am. Chem. Soc. 82, 5081 
5084 (1960). 

Absolute isotopic abundance ratio and the atomic weight 
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H. S. Parker, and M. D. Burdick, J. Am. Ceram. Soc. 43, 
No. 12, 641 647 (Dec. 1960). 

Green and purple sulfur: Electron-spin resonance studies, 
H. E. Radford and F, O. Rice, J. Chem. Phys. 33, No. 3, 
774-776 (Sept. 1960). 
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C. W. Gartlein, J. Geophys. Research 65, No. 9, 2759 
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